AD-A231  588 


■  ‘  •  M(:  0  3  3S  2-  3-  - 

? ;  5  . 5  . 

L  r 

EIGENSPACE  DESIGN  OF  HELICOPTER  FLIGHT 

CONTROL  SYSTEMS 

FINAL  REPORT 

WILLIAM  L.  GARRARD 
EICHER  LOW 

NOVEMBER,  1990 

U.  S.  ARMY  RESEARCH  OFFICE 

DAAL03-86-K-0056 


DEPARTMENT  OF  AEROSPACE  ENGINEERING 

AND  MECHANICS 
UNIVERSITY  OF  MINNESOTA 
MINNEAPOLIS,  MINNESOTA  55455 


APPROVED  FOR  PUBLIC  RELEASE; 
DISTRIBUTION  UNLIMITED. 


THE  VIEW,  OPINIONS,  AND/OR  FINDINGS  CONTAINED  IN  THIS 
REPORT  ARE  THOSE  OF  THE  AUTHOR(S)  AND  SHOULD  NOT  BE 
CONSTRUED  AS  AN  OFFICIAL  DEPARTMENT  OF  THE  ARMY 
POSITION,  POLICY,  OR  DECISION,  UNLESS  SO  DESIGNED  BY 
OTHER  DOCUMENTATION. 


MASTER  COPY 


KEEP  THIS  COPY  FOR  REPRODUCTION  PURPOSES 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No  0704-0188 


5.  FUNDING  NUMBERS 


'C' \  C^CTC"  z1  c"  •»  •;  *.<*Mqe  ‘  'Ouf  3f '  ''Coding  tr-e  ti^e  'O'  ir^if^cr.ons.  vr*rc">irig  e*<it<"3  data  vow'cev 

;jfhef>r'q  jirtji«<nq  aata  ne-aed  Aoa  jon'o*e,.*','q  j"d  'fv^^c  'm‘t  ccUecucd  0*  •“  .-'T^iion  Send  r*z*fc>rq  it*.*  Cuf3*n  wti^ate  o*  *«v  ot^ff  a%c*a  ot  t*n 

*qu#rtion  cf  int?c«ijtion.  n<’ufl;nq  ^gq^Mt'O^s  *o<  recucipq  ?"•*  ogre* ~  *:  Aa^h.^qtc"  •"■eaaa  ware's  5e«-v«c*v  Cir-aD'afe  *c'  r'v^auc"  ODe'atiop's  a-o  ae:>orts.  ' 2  '5  ;etfeY\ori 
Oav%  «.qr^a#.  Swte  *204  -r.rqicn  ,i  ;^02-4]0«'  ♦'■■a  to  t*e  0*«<e  :*  *.«  a^age****^  if'C  3ucqer  *>ac>er«*or>i  Seagate"  P'Cjer:  3>C4-0  *S9>  Aasn.nqicn  ;c  *350 


1.  AGENCY  USE  ONLY  (Leave  blank)  2.  REPORT  OATE  3.  REPORT  TYPE  AND  DATES  COVERED 

December  1990 


4  TITLE  AND  SUBTITLE 

Eigenspace  Design  of  Helicopter  Flight  Control 
Systems,  Final  Report 


6.  AUTHOR(S) 

William  L.  Garrard  and  Eicher  Low 


4)  ft  ft  L  0  3-#(,  -k~oOS 


7.  PERFORMING  ORGANIZATION  NAME(S)  ANO  ADDRESS(ES) 

Department  of  Aerospace  Engineering  and  Mechanics 
University  of  Minnesota 
107  Akerman  Hall,  110  Union  St.  SM 
Minneapolis,  MN  55455- 


9.  SPONSORING /MONITORING  AGENCY  NAME(S)  ANO  AOORESS(ES) 

U.  S.  Army  Research  Office 
P.  0.  Box  12211 

Research  Triangle  Park,  NC  27709-2211 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


10.  SPONSORING/ MONITORING 
AGENCY  REPORT  NUMBER 


()U  <2  3v<U  fc_<- 


11.  SUPPLEMENTARY  NOTES 

The  view,  opinions  and/or  findings  contained  in  this  report  are  those  of  the 
author (s)  and  should  not  be  construed  as  an  official  Department  of  the  Army 
position,  policy,  or  decision,  unless  so  designated  by  other  documentation. 


IZa.  DISTRIBUTION /AVAILABILITY  STATEMENT  12b.  DISTRIBUTION  CODE 

Approved  for  public  release;  distribution  unlimited. 


13.  ABSTRACT  (Maximum  200  words) 

An  eigenstructure  based  design  methodology  for  helicopter  flight 
control  systems  is  developed  and  evaluated  in  this  report.  The 
report  provides  a  detailed  review  of  the  application  of  multivariable 
design  techniques  to  helicopter  flight  control  systems,  a  review  of 
applicable  handling  quality  specifications,  a  description  of  the 
mathematical  models  used,  a  presentation  of  the  theory  of  eigen¬ 
structure  design,  application  of  the  theory  to  a  helicopter  design 
problem,  and  evaluation  of  the  performance  and  robustness  properties 
of  the  resulting  control  laws. 


14.  SUBJECT  TERMS 


Helicopter,  Flight  Control  Handling  Qualities, 
Eigenspace  Assignment 


15.  NUMBER  OF  PAGES 
181 


16.  PRICE  COOE 


17  SECURITY  CLASSIFICATION  18  SECURITY  CLASSIFICATION  119.  SECURITY  CLASSIFICATION  20.  LIMIT  A  TlON  OF  ABSTRACT 
OF  REPORT  OF  THIS  PAGE  OF  ABSTRACT 


UNCLASSIFIED 


\5N  7540-0!  290-5500 


UNCLASSIFIED 


UNCLASSIFIED 


89 


- 11  - 


TABLE  OF  CONTENTS 

Page 


Title  .  i 

Table  of  Contents  .  ’i 

List  of  Figures  .  v 

List  of  Tables  .  x 

Summary  .  xii 

Publications  .  xiii 

Participants  .  xv 

Chapter  1.  INTRODUCTION  .  1 

I.  Motivation  .  1 

II.  Previous  Work  .  7 

Model  Following  Concepts  .  9 

Quadratic  Optimal  Cooperative  -  Control  Synthesis  .  12 

MIMO  Proportional  Plus  Integral  (P  +  I)  Control  Law  .  13 

Linear  Quadratic  Regulator  (  LQR  )  /  Linear  Quadratic  Gaussian 

(  LQG  )  Methods  .  14 

H  Infinity  (  H«> )  Design  Techniques  .  15 

Eigenstructure  Assignment  .  20 

III.  Organization  of  the  Thesis  .  22 

Chapter  2.  HELICOPTER  HANDLING  QUALITIES  .  33 

I.  Introduction  .  33 

II.  Specification  Overview  .  34 

Mission  Task  Element  (  MTE  )  34 

Usable  Cue  Environment  (  UCE  )  35 

Pilot  Attention  .  35 

Speed  Range  .  35 

III.  Specification  Criteria  .  35 

A.  Bandwidth  Criterion  .  35 

B.  Response  -  Types  .  39 


-  iii  - 

C.  Dynamic  Requirements,  Hover  and  Low  Speed  .  41 

1.  Small  Amplitude  Attitude  Changes  (Pitch  and  Roll)  .  41 

2.  Moderate  Amplitude  Attitude  Changes  (Pitch  and  Roll)  ....  44 

3.  Large  Amplitude  Attitude  Changes  (Pitch  and  Roll)  .  44 

4.  Inter-Axis  Coupling  .  45 

5.  Response  to  Disturbance  .  46 

Chapter  3.  THEORY  OF  EIGENSTRUCTURE  ASSIGNMENT  .  48 

I.  Introduction  .  48 

II.  Full  State  Feedback  .  49 

III.  Full  State  Estimator  for  LTR  .  53 

IV.  Asymptotic  Linear  Quadratic  Regulator  Theory  .  56 

A.  Asymptotically  Finite  Modes  .  59 

B.  Asymptotically  Infinite  Modes  .  59 

Chapter  4.  MATHEMATICAL  MODEL  .  66 

I.  Introduction  .  66 

II.  Equations  of  Motion  .  66 

A.  Rigid  Body  Dynamics  .  66 

B.  Rotor  Dynamics  .  73 

III.  Design  Model  .  77 

IV.  Error  Model  .  84 

V.  Simulation  Model  .  85 

Chapter  5.  NOMINAL  DESIGN  .  88 

I.  Introduction  . .  88 

II.  Design  Philosophy  .  88 

A.  Inner  Loop  .  90 

B.  Outer  Loop  .  93 

III.  Inner  Loops  via  Eigenstructure  Assignment  .  94 

A.  Full  State  Feedback  .  94 

1 .  Nominal  Design  .  94 

B.  Full  State  Estimator  .  108 

IV.  Outer  Loop  .  113 

A.  ACAH  113 


-  iv  - 

B.  TRCPH  .  129 

Chapter  6.  ROBUSTNESS  ANALYSIS  .  134 

I.  Introduction  .  134 

II.  Unstructured  Singular  Value  Analysis  .  135 

A.  Stability  Robustness  .  135 

B.  Nominal  Performance  .  136 

III.  Structured  Singular  Value  Analysis  .  UtO 

A.  Theory  .  140 

B.  Mathematical  Formulation  .  143 

C.  Robust  Stability  .  146 

D.  Nominal  Performance  .  147 

E.  Robust  Performance  .  148 

F.  Weighting  Functions  .  149 

G.  Results  .  152 

IV.  Simulations  . 155 

A.  Analysis  of  Worst  Case  Perturbation  .  155 

B.  Time  Histories  with  High  Order  Model  .  164 

Appendices 

A.  Dimensional  12th  and  8th  Order  Matrices  .  182 

B.  Closed  Inner  Loop  Transfer  Function  Matrix  .  186 


Aooesslon  tot 


HTIS  GRAM 
DTIC  TAB 
Unannounced 
Justification. 


i n 


□ 

□ 


By - - 

Distribution/ 


Availability^  Codea_ 
Avail  and/or 
Dlst  1  Special 


-  V  - 


LIST  OF  FIGURES 

Figure  Page 

1.1  ADOCS  Modeling  Following  Control  System .  9 

1.2  Cooperative  Control  Synthesis  Design  Structure .  12 

1.3  Typical  H Control  Structure .  16 

2.1  Definition  of  Gain  Margin  Limited  Bandwidth .  36 

2.2  Definition  of  Phase  Margin  Limited  Bandwidth .  37 

2.3  Requirements  for  Small  Amplitude  Attitude  Changes .  42 

2.4  Limits  on  Mid-Term  Eigenvalues.  .  43 

2.5  Requirements  for  Moderate  Amplitude  Roll  Attitude  Changes .  44 

2.6  Yaw  Cross  Coupling  Criteria .  45 

4.1  Definition  of  the  Rigid  Body  Degrees  of  Freedom  of  a  Helicopter, 

and  the  Forces  and  Moments  due  to  the  Main  Rotor,  Tail  Rotor  and 
Gravity .  67 

4.2  Blade  Flapping  Harmonics.  .  74 

4.3  Rotor  Force  Diagram .  75 

4.4  Singular  Value  Plot  of  12th  Order  Open  Loop  Helicopter  at  Hover _  78 

4.5  Singular  Value  Plot  of  8th  Order  Residualized  Helicopter  Model.  ...  82 

4.6  Block  Diagram  of  the  Simulation  Model .  85 

4.7  Actuator  Dynamics.  86 

5.1  Singular  Value  Plot  of  Open  Loop  Helicopter  at  Hover .  90 

5.2  Integral  Characteristics  of  the  Desired  Transfer  Functions .  93 

5.3  Full  State  Open  and  Closed  Inner  Loop  Time  Responses  to  a  Step 

Roll  Rate  Command .  101 


-  vi  - 

5.4  Full  State  Open  and  Closed  Inner  Loop  Time  Responses  to  a  Step 

Roll  Rate  Command .  102 

5.5  Bode  Magnitude  of  Roll  Angle  to  Lateral  Command .  103 

5.6  Singular  Values  of  [  I  -  Gclosed  Inner  Loop  ] .  104 

5.7  5.7  Singular  Values  of  a  [  GcIosed  Inner  ,oop  ].  105 

5  8  Full  State  Closed  Inner  Loop  Response  of  Accelerations  to  Vertical 

Gust .  106 

5.9  Block  Diagram  for  Inner  and  Outer  Loops  Feedback  Configuration....  107 

5.10  o  [  I  +  (  K(s)G(s)  )'’  ]  Full  State  Inner  Loop  and  Estimator  in  the 

Feedback  Loop .  110 

5.11  Singular  Value  Plot  of  Compensator  (  Regulator/ Estimator ) .  113 

5.12  Definition  of  Phase  Margin  Limited  Bandwidth .  114 

5.13  Full  State  Closed  Inner  and  Outer  Loop  Time  Responses  to  a  Step 

Roll  Command .  117 

5.14  Full  State  Closed  Inner  and  Outer  Loop  Time  Responses  to  a  Step 

Roll  Command .  118 

5.15  Full  State  Closed  Outer  Loop  Response  to  a  5  degrees  Step  Roll 

Command  of  Duration  4  seconds.  .  119 

5.16  Control  Actuator  Deflections.  .  121 

5.17  Full  State  Closed  Inner  and  Outer  Loop  Time  Responses  to  a  5  ft./s 

Heave  Command.  122 

5.18  Full  State  Closed  Inner  and  Outer  Loop  Time  Responses  to  a  5  ft./s 

Heave  Command.  123 

5.19  Full  State  Closed  Inner  and  Outer  Loop  Time  Responses  to  a 

5  deg./s  Yaw  Rate  Command .  124 


-  vii  - 

5.20  Full  State  Closed  Inner  and  Outer  Loop  Time  Responses  to  a 

5  deg./s  Yaw  Rate  Command .  125 

5.21  Phase  Delay  of  ACAH  System  plus  Effective  Time  Delay .  127 

5.22  Level  1  Small-Amplitude  Requirements  for  Minimum  and  Maximum 

Gain  Tasks.  .  127 

5.23  Closed  Outer  Loop  Full  State  Roll  Attitude  Bode  Magnitude  and 

Phase  Plots .  128 

5.24  Block  Diagram  for  Forward  Velocity  TRCPH .  130 

5.25  Block  Diagram  for  Position  Control .  130 

6.1  <2.  [  I  +  (  K(s)G(s)  )''  ]  Full  State  Inner  Loop  and  Estimator  in  the 

Feedback  Loop .  137 

6.2  Minimum  Singular  Value  for  Full  State  Closed  Outer  Loop  and 

Maximum  Singular  Value  of  Error  Matrix  vs.  Frequency .  138 

6.2  Perfcrni"-"?  Test  for  Cloned  Inner  Loop  System .  139 

6.4  — -  ofGCjami  +  KijwiGUo)))-1  ].  .  140 

I  r(ja>)  I 

6.5  General  Interconnection  Structure.  .  141 

6.6  Robust  Performance  Block  Diagram  Description .  142 

6.7  Perturbed  Pilot  Stick  Inputs  to  Errors .  143 

6.8  Matrix-Vector  Multiplication.  .  143 

6.9  Block  Diagram  for  the  Standard  Feedback  Loop .  144 

6.10  Partitioned  Matrix-Vector  Multiplication .  144 

6.11  Linear  Fractional  Transformation .  145 

6.12  Linear  Fractional  Transformation  for  Nominal  Performance .  148 

6.13  a(  Mj j  ],  Maximum  Singular  Value  for  Robust  Stability.  ..  153 


-  viii  - 

6.14  o  [  Mo  ],  Maximum  Singular  Value  for  Nominal  Performance .  154 

6.15  n  [  M!  ],  Structured  Singular  Value  for  Robust  Performance .  155 

6. 16  Full  State  Closed  Inner  and  Outer  Loop  Time  Responses  to  a 

5  Degrees  Step  Roll  Command  with  Parameter  Variation  of  -  20%  in 
A42.  158 

6.17  Full  State  Closed  Inner  and  Outer  Loop  Time  Responses  to  a 

5  Degrees  Step  Roll  Command  w'ith  Parameter  Variation  of  -  20%  in 

A.  :.  159 

6.18  Full  State  Closed  Inner  and  Outer  Loop  Time  Responses  to  a 

5  Degrees  Step  Roll  Command  with  Parameter  Variation  of  +  20%  in 

B42.  160 

6.19  Full  State  Closed  Inner  and  Outer  Loop  Time  Responses  to  a 

5  Degrees  Step  Roll  Command  with  Parameter  Variation  of  +  20%  in 
B42.  161 

6.20  Full  State  Closed  Inner  and  Outer  Loop  Time  Responses  to  a 

5  Degrees  Step  Roll  Command  with  Parameter  Variation  of  -  20%  in 

A4  2  and  +  20%  in  B4  2.  162 

6.21  Full  State  Closed  Inner  and  Outer  Loop  Time  Responses  to  a 

5  Degrees  Step  Roll  Command  with  Parameter  Variation  of  -  20%  in 

A4  j  ar*d  +  20%  in  B4  2.  163 

6.22  12th  Order  and  8th  Order  Closed  Loop  Responses  to  a  5  ft./s  Step 

Vertical  Velocity  Command .  165 

6.23  12th  Order  and  8th  Order  Closed  Loop  Responses  to  a  5  ft./s  Step 

Vertical  Velocity  Command .  166 


-  ix  - 

6.24  12th  Order  and  8th  Order  Closed  Loop  Responses  to  a  5  ft./s  Step 

Vertical  Velocity  Command .  167 

6.25  12th  Order  and  8th  Order  Closed  Loop  Responses  to  a  5  Degrees 

Step  Roll  Attitude  Command.  .  168 

6.26  12th  Order  and  8th  Order  Closed  Loop  Responses  to  a  5  Degrees 

Step  Roll  Attitude  Command.  .  169 

6.27  12th  Order  and  8th  Order  Closed  Loop  Responses  to  a  5  Degrees 

Step  Roll  Attitude  Command.  .  170 

6.28  12th  Order  and  8th  Order  Closed  Loop  Responses  to  a  5  Degrees 

Step  Pitch  Attitude  Command.  .  172 

6.29  12th  Order  and  8th  Order  Closed  Loop  Responses  to  a  5  Degrees 

Step  Pitch  Attitude  Command.  .  173 

6.30  12th  Order  and  8th  Order  Closed  Loop  Responses  to  a  5  Degrees 

Step  Pitch  Attitude  Command.  .  174 

6.31  12th  Order  and  8th  Order  Closed  Loop  Responses  to  a  5  deg./s  Step 

Yaw  Rate  Command .  175 

6.32  12th  Order  and  8th  Order  Closed  Loop  Responses  to  a  5  deg./s  Step 

Yaw  Rate  Command .  176 

6.33  12th  Order  and  8th  Order  Closed  Loop  Responses  to  a  5  deg./s  Step 

Yaw  Rate  Command .  177 


-  X  - 


LIST  OF  TABLES 

Table  Page 

2.1  Requirements  for  Large  Amplitude  Attitude  Changes .  45 

4.1  Open  Loop  Eigenvalues  and  Corresponding  Modes  of  12th  Order 

Model .  79 

4.2  Comparison  of  Open  Loop  Eigenvalues  of  12th  Order  Model  with  8th 

Order  Truncated  and  Residualized  ModcL .  81 

4.3  Open  Loop  Eigenvalues  and  Associated  Non-dimensional 

Eigenvectors.  83 

5.1a  Desired  Closed  Loop  Eigenvalues.  .  96 

5.1b  Desired  Closed  Loop  Eigenvectors.  .  97 

5.2a  Achievable  Non-Dimensional  Regulator  Eigenvectors .  98 

5.2b  Non-Dimensional  Regulator  Feedback  Gain  Matrix .  98 

5.3a  Non-Dimensional  Feedforward  Gain  Matrix .  99 

5.3b  Achievable  Form  of  Non-Dimensional  Bd .  99 

5.4  Desired  Estimator  Eigenvalues .  Ill 

5.5  Desired  Estimator  Eigenstructure.  .  Ill 

5.6  Non-Dimensional  Estimator  Gain  Matrix .  112 

6.1  Summary  of  Time  Delay  Contributors .  136 

6.2  Comparison  of  Closed  Inner  Loop  Eigenvalues  of  the  Nominal 

Design  with  Parameter  Variations .  156 

6.3  Comparison  of  Closed  Outer  Loop  Eigenvalues  of  the  Nominal 

Design  with  Parameter  Variations .  157 


-  XI  - 


6.4  Comparison  of  Closed  Inner  and  Outer  Loop  Eigenvalues  of  the  12th 
Order  and  8th  Order  Models  Using  the  Same  Full  State  Regulator — 


164 


-  XU  - 


EIGENSPACE  DESIGN  OF  HELICOPTER  FLIGHT  CONTROL  SYSTEMS 


SUMMARY 

An  eigenstructure  based  design  methodology  for  helicopter  flight  control  systems  is 
developed  and  evaluated  in  this  report.  The  report  provides  a  detailed  review  of  the 
application  of  multivariable  design  techniques  to  helicopter  flight  control  systems,  a  review 
of  applicable  handling  quality  specifications,  a  description  of  the  mathematical  models 
used,  a  presentation  of  the  theory  of  eigenstructure  design,  application  of  the  theory  to  a 
helicopter  design  problem,  and  evaluation  of  the  performance  and  robustness  properties  of 
the  resulting  control  laws. 

The  methodology  detailed  in  this  report  allows  the  designer  to  synthesize  control  laws 
which  result  in  desirable  response  types  such  as  attitude  command  attitude  hold,  rate 
command,  and  translational  rate  command  position  hold.  Eigenstructure  assignment  is 
used  to  design  inner  loop  control  laws  which  provide  decoupled  first  order  rate  command 
response  in  each  control  axis.  Outer  loop  control  laws  are  then  designed  to  achieve  desired 
response  types  using  standard  single  input  single  output  frequency  response  design 
techniques. 

Control  laws  were  designed  for  a  typical  high  performance  helicopter  operating  in  the 
hover  flight  condition.  The  performance  of  these  control  laws  was  evaluated  with  the 
mathematical  models  used  for  design  and  the  control  laws  were  found  to  yield  excellent 
dynamic  response  characteristics  in  both  the  time  and  frequency  domains.  The  control  laws 
yielded  good  tracking  of  pilot  commands  and  attenuation  of  atmospheric  disturbances. 
They  also  provided  Level  1  handling  qualities  when  evaluated  using  the  design  model. 

The  stability  robustness  of  the  control  laws  in  the  presence  of  modeling  uncertainties 
was  evaluated  using  unstructured  singular  values.  A  simple  time  delay  model  of  the  errors 
due  to  neglected  high  frequency  dynamics  was  used  in  this  analysis  and  the  stability 
robustness  of  the  system  was  shown  to  be  acceptable.  Much  more  detailed  robustness 
analyses  were  accomplished  using  structured  singular  value  techniques  and  computer 
simulations.  Actuator  dynamics,  rotor  dynamics,  time  delays  due  to  digital  implementation 
of  flight  control  laws,  and  errors  in  aerodynamic  stability  and  control  derivatives  were 
included  in  these  robustness  analyses.  The  stability  and  performance  robustness  of  the 
feedback  control  laws  were  shown  to  be  acceptable  and  in  no  case  were  handling  qualities 
less  than  Level  2  obtained. 

Eigenstructure  assignment  appears  to  provide  a  good  method  for  design  of  helicopter 
flight  control  systems  which  result  improved  handling  qualities.  Eigenstructure  assignment 
provides  a  straightforward  methodology  which  results  in  control  laws  which  are  relatively 
simple  to  implement  and  which  do  not  require  high  order  dynamic  compensators  in  the 
feedback  loops..  The  eigenstructure  control  laws  provide  good  robustness  properties.  It  is 
suggested  that  implementation  of  these  control  laws  on  a  piloted  simulator  be  investigated 
with  the  goal  of  eventual  flight  tests. 
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Chapter  1 


Introduction 


I.  Motivation 

The  objective  of  the  research  described  in  this  thesis  is  to  develop  and  evaluate  a 
methodology,  based  on  eigenstructure  assignment,  for  the  design  of  advanced  flight  control 
laws  suitable  for  modem  high  performance  helicopters. 

When  operating  at  hover  or  low  speed  (  up  to  approximately  45  knots  )  helicopters 
exhibit  poor  handling  qualities  and  are  difficult  to  fly.  Interactions  between  the  longitudinal 
and  lateral  axes,  coupling  of  the  control  inputs,  and  inherent  low  frequency  instabilities  are 
undesirable  characteristics  typical  of  helicopters  when  flying  at  low  speeds.  As  a  result, 
helicopters  can  benefit  substantially  from  the  use  of  automatic  flight  control  systems  to 
enhance  Km. ’ling  qualities.  Until  recently,  the  development  of  rotorcraft  flight  controls  has 
lagged  behind  more  advanced  automatic  flight  controls  which  are  available  for  fixed  wing 
aircraft  [l]-[3]. 

Recent  advances  in  rotorcraft  technology  have  greatly  expanded  the  role  of  both 
military  and  civilian  helicopters.  This  expanded  role  has  resulted  in  a  requirement  for 
improved  handling  qualities.  Until  the  niddle  1980's,  helicopter  handling  quality 
specifications  had  been  poorly  defined  compared  with  those  for  fixed  wing  aircraft.  As 
part  of  the  U.S.  Army  LH  (  Light  Helicopter  )  program,  an  extensive  effort  to  revise 
military  handling  quality  specifications  was  under  taken.  This  program  is  essentially 
completed  and  has  resulted  in  a  comprehensive  set  of  rotorcraft  handbag  qualities.  These 
new  handling  quality  specifications  are  mission  oriented  and  require  that  advanced 

rotorcraft  exhibit  good  handling  qualities  over  a  diverse  spectrum  of  operational  missions. 
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Precision  flight  path  control  and  attitude  stability  are  required  for  single  pilot  operation  in 
hovering  and  flight  close  to  the  earth,  especially  in  poor  weather  conditions  and/  or  at  night. 
For  rotorcraft  which  are  intended  for  air-to-air  combat,  agility  and  high  maneuverability  are 
required.  These  new  handling  qualities  specifications  are  being  used  for  civilian  as  well  as 
military  helicopters. 

The  revised  handling  quality  specifications  are  structured  in  terms  of  specific  Mission- 
Task-Elements  (  MTEs  )  and  Operational  Flight  Envelopes  (  OFEs  ).  The  specifications  are 
expressed  in  terms  of  dynamic  response  to  control  inputs  or  response  types  [4],  These 
response  types  are  formulated  as  desired  frequency  response  characteristics  between  pilot 
inputs  and  vehicle  outputs.  These  outputs  may  include  attitude  angles,  attitude  rates,  and/or 
translational  velocities.  Important  characteristics  of  these  frequency  response 
characteristics  are  their  bandwidth,  magnitude,  and  phase.  If  these  desired  frequency 
response  characteristics  are  modeled  as  either  first  order  or  second  order  transfer  functions, 
it  is  simple  to  relate  the  desired  frequency  response  characteristics  to  the  natural  frequencies 
and  damping  factors  or  the  time  constants  of  the  transfer  function  models. 

The  revised  military  specifications  provide  a  unified  framework  for  design  and  for 
evaluation  of  handling  qualities.  Requirements  and  limits  are  set  at  levels  which  if  not  met, 
will  result  in  poor  flying  qualities.  Thus  the  requirements  are  necessary,  but  may  not  be 
sufficient,  to  guarantee  Level  1  handling  qualities.  Minimum  Level  1  control  response 
bandwidths  for  rate  and  attitude  response  types  ir  the  pitch  and  roll  axes  have  been  defined 
and  incorporated  into  the  new  military  specifications.  These  bandwidths  have  been 
adjusted  to  reflect  the  requirements  imposed  by  various  mission-task-elements. 

Three  primary  response  types  [4]-[5]  have  been  proposed.  These  are  Attitude 
Command/  Attitude  Hold  (  ACAH  ),  Rate  Command  (  RC  )  or  Rate  Command/  Attitude 
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Hold  (  RCAH  ),  and  Translation  Rate  Command/  Position  Hold  (  TRCPH  ).  The 
characteristics  for  a  step  control  input  for  each  of  these  response  types  are  as  follows: 

1.  The  ACAH  response  type  results  in  a  pitch  ( roll )  attitude  angle  proportional  to  the 
longitudinal  ( lateral )  stick  input.  When  the  input  is  removed,  the  helicopter  returns 
to  the  trim  attitude  before  the  input  was  applied.  This  response  type  typically  does 
not  have  the  agility  of  the  RCAH  type  but  is  highly  desirable  in  low  visibility 
conditions,  especially  in  hover,  since  releasing  the  control  input  will  return  the 
rotorcraft  to  its  reference  trim  attitude.  A  trim  control  of  some  type  is  required  in 
order  to  change  the  trim  attitude  reference. 

2.  For  the  RCAH  response  type,  angular  velocities  about  the  vehicle  roll,  pitch,  or 
yaw  axis  are  proportional  to  the  corresponding  cockpit  controller  force  and/or 
deflection.  The  RCAH  response  type  is  characterized  as  having  a  high  level  of 
agility  with  capability  for  extended  periods  of  hands-off  operation. 

3.  For  the  TRCPH  response  type,  a  constant  controller  force  input  must  result  in 
constant  translational  rate.  The  rotorcraft  must  hold  position  if  the  force  on  the 
cockpit  controller  is  zero.  The  TRCPH  response  type  is  useful  for  precision 
hovering  tasks  and  provides  an  important  tactile  cue  for  low  speed  maneuvering  in 
poor  visual  cueing. 

Unaugmented  helicopters  provide  control  responses  which  are  of  the  rate  command 
type.  Typical  open  loop  rigid  body  helicopter  dynamics  include 

1.  extensive  coupling  between  the  lateral  and  longitudinal  responses, 

2.  low  frequency  unstable  response  in  forward  and  side  velocities. 
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3.  insufficient  bandwidth  for  Level  1  flying  qualities  in  all  four  axes. 


Most  unaugmented  rotorcraft  will  not  meet  the  new  handling  quality  specifications  and 
feedback  control  systems  are  necessary  to  improve  handling  qualities  so  that  safe  operation 
close  to  the  earth  in  poor  weather  conditions  and/or  at  night  is  possible.  Control  systems 
should  result  in  the  following  characteristics: 

1.  accurate  tracking  of  pilot  commands, 

2.  high  bandwidth  and  low  time  delay  response  characteristics, 

3.  attenuation  of  gusts  and  other  disturbances, 

4.  achievement  of  specified  response  types, 

5.  decreased  pilot  workload, 

6.  modal  decoupling, 

7.  acceptable  stability  and  performance  robustness. 

Typical  helicopters  have  four  control  inputs.  For  a  single  rotor  vehicle,  these  are  as 
follows: 

1.  Collective  pitch  which  controls  the  lift  produced  by  the  main  rotor  and  is  the 
primary  control  for  vertical  velocity. 

2.  Lateral  cyclic  pitch  which  controls  the  rolling  moment  produced  by  the  main  rotor 
and  is  the  primary  control  for  roll  rate. 
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3.  Longitudinal  cyclic  pitch  which  controls  the  pitching  moment  produced  by  the  main 
rotor  and  is  the  primary  control  for  pitch  rate. 

4.  Tail  rotor  collective  pitch  which  controls  the  lift  produced  by  the  tail  rotor  and 
consequently  the  yawing  moment  on  the  helicopter.  This  is  the  primary  control  for 
yaw  rate. 

Typical  flight  control  sensors  measure  normal  acceleration  or  velocity,  roll,  pitch,  and 
yaw  rates,  and  side  and  forward  velocity.  Other  measurements  can  include  pitch,  roll  and 
heading  angles  and  altitude.  There  is  considerable  coupling  between  inputs  and  outputs. 
For  example,  a  collective  input  will  create  yawing  moment  as  well  as  lift  and  a  longitudinal 
cyclic  input  will  produce  a  rolling  moment  as  well  as  a  pitching  moment.  Part  of  the  pilot's 
task  is  to  compensate  for  these  couplings.  This  results  in  high  pilot  work  load,  particularly 
during  hover.  The  principal  source  of  coupling  is  due  to  the  main  rotor  which  is  required 
to  produce  both  lift,  forward  thrust,  and  pitching  and  rolling  moments  for  control. 
Production  of  thrust  and  control  moments  requires  changes  in  the  rotor  tip  path  plane 
orientation.  These  changes  in  tip  path  plane  orientation  produce  gyroscopic  effects  which 
influence  all  of  the  forces  and  moments  acting  on  the  helicopter. 

Due  to  the  fact  that  responses  to  control  inputs  are  highly  coupled,  helicopter  dynamics 
are  characterized  by  multi-input-multi-output  (  MIMO  )  mathematical  models  and  design  of 
flight  control  systems  for  such  vehicles  is  a  true  multivariable  synthesis  problem.  Modem 
MIMO  design  techniques  are  well  suited  for  the  design  of  rotorcraft  control  laws. 

Typical  mathematical  models  used  for  helicopter  control  system  design  are  based  on 
rigid  body  models  of  the  fuselage  dynamics  with  a  quasi-static  representations  of  the  rotor 
dynamics  [6]-[7].  This  representation  is  appropriate  as  the  state  variables  representing  the 
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rigid  body  motions  of  the  helicopter  are  the  quantities  which  are  to  be  controlled;  however 
effects  such  as  rotor,  structural,  actuator,  sensor,  and  engine  dynamics  can  have  a 
significant  effect  on  dynamics  response  [8]. 

Rotor  dynamics  can  produce  non-minimum  phase  zeros  which  can  result  in 
instabilities  if  control  gains  are  too  large.  At  present  it  appears  that  information  on  rotor 
states  cannot  be  used  in  control  laws  due  to  noise  and  other  difficulties  associated  in  the 
measurement  of  these  states. 

Actuators  produce  phase  lags  which  can  result  in  instabilities  and  also  have  both  rate 
and  deflection  limits.  Since  saturation  of  actuators  is  undesirable,  inputs  to  the  actuators 
should  be  kept  small  enough  to  avoid  saturation  and  this  may  also  limit  feedback  gains. 
Also,  since  the  actuators  themselves  have  a  limited  bandwidth  (  30  to  60  rad/s  ),  the  control 
laws  should  not  require  high  frequency  actuator  responses. 

Interaction  between  the  rotor  and  the  airframe  generates  what  are  called  "n/rev" 
vibrations  (  n  equals  the  number  of  blades  ).  The  sensors  used  to  measure  the  vehicle 
outputs  will  pick  up  these  vibrations  and  the  control  system  will  be  excited  by  vibrations  at 
this  frequency.  For  example,  a  typical  helicopter  rotor  operates  at  an  angular  velocity  of 
about  30  rad/s,  so  for  a  four-bladed  helicopter,  the  sensors  would  be  subjected  to  excitation 
of  120  rad/s.  Thus  sensor  outputs  must  be  filtered  before  they  are  fedback  into  the  control 
system.  Typically  a  Bessel  filter  is  used  for  this  purpose,  as  this  type  of  filter  gives  the 
least  phase  shift  of  any  type  of  linear  filter  in  the  band  pass  region.  This  filter  does 
introduce  significant  phase  shift  at  frequencies  higher  than  the  design  filter  bandwidth. 
Thus  the  bandwidth  of  the  feedback  compensator  is  limited  to  some  fraction  of  n  times  the 
rotor  angular  velocity.  Typically,  compensators  should  be  designed  with  a  bandwidth  of 
no  more  than  10  rad/s. 
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Structural  dynamics  of  the  fuselage  may  also  affect  stability  and  performance  and  lack 
of  accurate  models  for  these  modes  makes  it  difficult  to  include  them  directly  in  the  control 
laws.  Gas  turbine  engines  used  for  propulsion  introduce  both  lags  and  coupling  between 
yawing  motion  and  collective  pitch  changes.  Also,  it  is  important  to  avoid  large  changes  in 
rotor  angular  velocity  (  rotor  droop  )  when  control  inputs  are  applied. 

Rotor,  actuator  and  structural  dynamics  affect  the  upper  pan  range  of  frequencies  for 
flight  control  and  the  closed  loop  response  characteristics  must  be  insensitive  to  these 
effects  as  well  as  to  delays  due  to  filtering  and  digital  implementation  of  the  flight  control 
laws.  Uncertainties  in  aerodynamic  stability  and  control  coefficients  can  also  be  significant 
and  will  introduce  low  frequency  errors  in  the  mathematical  model  used  for  design. 

The  key  design  criteria  for  the  flight  control  systems  are 

1.  decouple  both  the  internal  dynamics  and  the  response  to  controller  inputs, 

2.  stabilize  the  helicopter, 

3.  increase  the  bandwidth, 

4.  provide  the  desired  response  type  or  types, 

5.  reject  atmospheric  disturbances  and  sensor  noise, 

6.  provide  acceptable  response  characteristic  even  if  the  dynamics  of  the  helicopter 
differ  significantly  from  those  of  the  mathematical  model  used  for  control  system 
design  (  stability  and  performance  robustness  ). 

II.  Previous  Work 
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Methods  used  for  design  of  helicopter  flight  control  systems  [9]-[  1 3]  include  classical 
single  input  single  output  (  SISO  )  approaches,  model  following  methods  [  14]-[  18],  Linear 
Quadratic  Regulator  (  LQR  )  and  Linear  Quadratic  Gaussian/  Loop  Transfer  Recovery 
(  LQG/LTR  )  procedures  [19]-[23],  dynamic  inversion  [24]-[26],  H-infinity  (  H ) 

techniques  [27]-[30],  and  eigenstructure  assignment  techniques  [31]-[34],  An  excellent 
review  of  research  in  the  area  of  multivariable  design  techniques  applied  to  helicopter  flight 
control  system  design  is  given  in  [35]. 

Classical  SISO  techniques  applied  to  each  control  loop  individually  were  used  to 
design  control  systems  for  all  current  production  rotorcraft.  These  systems  were  designed 
for  stability  augmentation  and  are  largely  mechanical  in  nature.  Classical  SISO  design 
techniques  concentrate  on  achieving  a  maximal  loop  gain  over  a  specified  frequency  range 
consistent  with  specified  gain  and  phase  margins.  Although  this  has  proven  to  be  a  very 
reliable  and  successful  technique,  classical  SISO  techniques  treat  each  control  loop 
separately  and  do  not  explicitly  include  coupling  between  loops.  Since  this  coupling  is 
very  important  in  rotorcraft,  a  variety  of  modem  techniques  are  being  studied  for  use  in 
advanced  rotorcraft  flight  control  systems. 

The  success  of  digital  flight  control  technology  for  improving  handling  qualities  of 
fixed  wing  aircraft  has  motivated  similar  effects  for  rotorcraft  [36]  and  at  least  six  vehicles 
have  been  flight  tested  using  advanced  digital  control  systems.  These  are  the  Boeing 
Advanced  Digital  Optical  Control  System  (  ADOCS  )  [37],  implemented  on  a  Sikorsky 
UH  60B,  the  McDonnell  Douglas  Apache  AV05  YAH-64,  the  German  DLR  MBB 
Bo  105,  the  NASA  Ames  Research  Center  CH-47  and  UH-IH,  and  the  Aerospatiale 
Dauphin.  The  ADOCS  and  DLR  Bo  105  used  explicit  model  following.  A  variety  of 
techniques  were  used  on  the  NASA  CH-47.  The  YAH-64  used  a  MIMO  proportional  plus 
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integral  (  P  +  I )  controller.  The  NASA  UH-1H  used  a  dynamic  inversion  control  and  the 
Dauphin  control  system  was  designed  using  eigenstructure  assignment  but  no  details  of  this 
design  are  available.  The  remainder  of  this  section  will  review  work  which  had  been  done 
so  far  in  the  application  of  multivariable  techniques  to  the  design  of  helicopter  flight  control 
systems. 

Model  Following  Concepts 

A  variety  of  model  following  concepts  have  been  proposed  for  use  in  flight  control 
systems.  The  ADOCS  design  is  representative  of  model  following  designs  and  will  be 
described  in  some  detail  below.  The  ADOCS  model  following  concept  uses  feedforward 
and  inverse  plant  dynamics  to  cancel  the  inherent  rotorcraft  dynamics  and  replace  them  with 
the  desired  command  responses  from  a  command  model.  Feedback  is  used  for  stability 
enhancement  and  reduction  of  sensitivity  to  variations  in  open  loop  dynamics. 
Feedforward  is  used  to  obtained  desired  response  to  commands.  Using  this  approach,  the 
command  and  stabilization  response  characteristics  can  be  set  independently.  A  block 
diagram  of  the  ADOCS  system  is  shown  in  Fig.  1.1. 


Fig.  1.1  ADOCS  Modeling  Following  Control  System. 

In  Fig.  1.1,  the  following  nomenclature  applies: 
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=  pilot  command  (  collective  pitch,  lateral  cyclic  pitch,  longitudinal  cyclic 
pitch,  or  tail  rotor  collective  pitch  ), 

=  shaped  pilot  command  to  achieve  desired  control  sensitivity, 

=  command  from  pre-compensator, 

=  command  from  feedback  loop, 

=  actuator  command, 

=  output  command, 

=  output, 

=  error  between  output  command, 

=  transfer  function  model  of  stable  short  period  helicopter  dynamics, 

=  desired  transfer  function  between  pilot  command  and  output, 

=  model  of  actuator  and  other  stable  system  dynamics  not  included  in  P(s), 

=  the  actuator  and  other  system  dynamics  not  included  in  P(s). 
block  diagram  manipulations  yield 


ecm  =  (  1  +  P(s)  A(s)  H(s)  )-*  (  A(s)  -  Am(s) )  D(s)  5s(s) 


If  A^s)  — »  A(s),  ecm  — »  0  and  the  output  approaches  the  output  command.  Furthermore, 
H(s)  can  be  selected  to  ensure  that  1  +  P(s)  A(s)  H(s)  has  desirable  stability  margins.  If 
H(s)  is  large  at  low  frequencies  where  command  following  is  desired,  then 


and  the  desired  transfer  function  between  input  and  output  will  be  achieved. 


For  the  ADOCS  system,  the  transfer  functions  for  controlling  yaw  rate  are  as  follows: 


1.  The  desired  response  is  First  order  with  a  bandwidth  coc 
K 


D(s)  = 


S  +  CD, 


2.  The  yaw  rate/  tail  rotor  collective  pitch  short  period  transfer  function  is  modeled  as 

p<s>  -  ;ttn7 


3.  The  actuator  and  other  stable  dynamics  not  included  in  P(s)  are  modeled  by 
a 


a„(s)  =  (r?-5)c 


-Ts 


where  a  is  the  actuator  roll  off  frequency  and  x  represents  a  time  delay 
approximation  of  rotor  dynamics,  filters,  computational  lags  etc. 

4.  The  feedback  control  laws  is  proportional  plus  integral 

H(s)  = 

The  ADOCS  model  following  approach  is  SISO  and  ignores  inter  axial  coupling.  Also 
since  P(s)  must  be  inverted,  it  must  be  minimum  phase  and  must  not  contain  unstable  poles 


-  11  - 


(  when  inverted  these  will  become  non-minimum  phase  zeros  ),  stable  poles  near  the 
imaginary  axis  (  when  inverted  these  will  become  near  non-minimum  phase  zeros  ),  non¬ 
minimum  phase  zeros  (  when  inverted  these  will  become  unstable  poles  ),  or  near  non¬ 
minimum  phase  zeros  (  when  inverted  these  will  become  near  unstable  poles  ).  Helicopters 
have  intrinsically  unstable  low  frequency  poles  which  must  be  neglected  in  the  inversion 
model,  thus  this  model  will  be  inaccurate  at  low  frequencies.  Despite  these  shortcomings, 
the  ADOCS  system  performed  adequately  in  flight  test. 


A  model  following  approach  similar  to  ADOCS  was  used  in  the  successful  flight 
testing  of  a  MBB  Bo  105  [38].  It  was  found  that  the  accurate  modeling  of  the  rotor 
dynamics  was  critical  in  the  success  of  this  design  [39]-[40],  This  model  following  design 
technique  was  also  used  to  design  control  laws  for  UH-1H  and  was  tested  on  the  NASA 
Ames  flight  simulator  [15],  [18].  The  UH-1H  and  Bo  105  have  completely  different  rotor 
designs  which  result  in  different  control  characteristics  for  the  two  helicopters. 


Quadratic  Optimal  Cooperative  -  Control  Synthesis 


Inputs 


Outputs 


SCAS  :  Stability  and  Control  Augmentation  System 


Fig.  1.2  Cooperative  Control  Synthesis  Design  Structure. 
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An  extension  of  LQR/LQG  methods,  the  quadratic  optimal  control  synthesis  (  CCS  ) 
method  developed  by  Innocenti  and  Schmidt  [9],  was  applied  by  Townsend  to  the  NASA 
CH-47  tandem  rotor  helicopter  [11].  The  basic  structure  of  the  CCS  control  system  is 
shown  in  Fig.  1.2.  The  CCS  method  incorporates  an  optimal  pilot  model  as  well  as  an 
aircraft  model  and  minimizes  the  weighted  sum  of  the  squares  of  the  aircraft  states,  the 
control  inputs,  and  the  states  of  the  pilot  model.  Townsend  claims  the  CCS  method  offers 
the  following  two  advantages  over  other  methods 

1 .  The  CCS  method  uses  output  feedback  rather  than  full  state  feedback,  thereby 
making  it  simpler  for  implementation. 

2.  The  CCS  method  requires,  in  principle,  no  detailed  a  priori  design  criteria,  because 
an  assumed  analytic  pilot-model  structure  is  an  inherent  feature  of  the  approach  as 
shown  in  Fig.  1.2. 

The  second  feature  can  be  an  advantage  in  cases  where  no  design  criteria  exist .  This 
is  not  the  case  with  helicopters  because  of  the  new  handling  qualities  specifications. 
Townsend  obtained  two  CCS  designs;  one  using  the  original  approach  of  Innocenti  and 
Schmidt  and  a  modified  approach  developed  by  specifically  for  helicopters.  Results  from 
fixed  base  piloted  simulations  indicated  that  the  modified  CCS  design  compared  favorably 
with  a  classical  frequency  response  design. 

MIMO  Proportional  plus  Integral  (  P  +  I  )  Control  Laws 

An  extremely  successful  approach  in  classical  SISO  design  is  the  used  of  a  P+I 

controller.  This  type  of  controller  gives  high  loop  gain  at  low  frequencies  for  good 

command  tracking  and  provides  good  high  frequency  roll  off  for  stability  robustness. 

MIMO  generalizations  of  the  P+I  technique  have  been  developed  and  applied  to  helicopters. 
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Enns  [10]  used  a  multivariable  low  frequency  approximation  of  the  YAH-64  dynamics  and 
designed  a  MIMO  P+  I  control  law  for  this  helicopter.  A  pre-compensator  was  used  to 
approximately  cancel  the  short  period  dynamics  and  achieve  desired  command  response 
characteristics  and  a  P+I  feedback  control  law  was  used  to  obtain  desirable  feedback 
properties  of  command  tracking,  stability  and  insensitivity  to  variation  in  forward  loop 
dynamics.  These  control  laws  were  successfully  flight  tested  [41].  A  similar  approach  is 
being  studied  for  use  in  the  Agusta  EH  101  helicopter  [42]. 

Linear  Quadratic  Regulator  (  LQR  )  /  Linear  Quadratic  Gaussian 

( LQG  )  Methods 

LQR  theory  results  in  a  linear  control  law  which  minimizes  the  integral  over  an  infinite 
time  interval  of  the  weighted  sum  of  the  squares  of  the  elements  of  the  system  state  and 
control  vectors.  The  LQR  control  law  requires  that  all  system  states  be  fedback  (  full  state 
feedback  )  but  results  in  excellent  multivariable  stability  margins  [20]-[21],  [43]-[45].  In 
most  applications,  all  of  the  system  states  cannot  be  measured  and  a  state  estimator 
(  Kalman  Filter )  is  required  in  the  feedback  loop.  This  estimator  can  seriously  degrade  the 
stability  margins  of  the  closed  loop  systems;  however,  in  the  late  1970's  Doyle  and  Stein 
developed  a  method  for  using  LQG  theory  to  recover  the  closed  loop  feedback  properties  of 
the  LQR.  This  procedure,  called  loop  transfer  recovery  (  LTR  ),  can  be  used  in 
conjunction  with  LQG  theory  to  shape  the  frequency  response  characteristics  of  the 
feedback  loops.  The  LQG/LTR  methodology  has  provided  an  important  tool  for  control 
law  design  and  is  commonly  used  in  design  of  flight  control  laws  for  fixed  wing  aircraft 
[46]. 

Several  early  papers  have  described  the  use  of  LQR  methods  for  helicopters  [47]-[49]. 
More  recent  results  have  used  the  LQG  [50]  or  LQG/LTR  [51]  methodologies  and  at  least 
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one  design  has  been  flight  tested  [50],  Application  of  LQG  loop  shaping  to  helicopter 
flight  control  was  studied  by  Holdridge,  et  al.  [50]  and  Rodriguez  [51].  Gribble  and 
Murray-Smith  [52]  have  recently  used  an  LQR  version  of  model  following  to  design  a 
flight  control  system  for  the  Westland  Lynx  helicopter  operating  at  a  forward  speed  of  80 
knots.  Their  design  was  based  on  a  rigid  body  model  with  the  actuator  and  blade  flapping 
dynamics  were  represented  as  pure  delays  in  order  to  motivate  the  choice  of  the  LQR  input 
weighting  matrix.  Prasad  et  al.  [53]  present  a  design  approach  based  on  output  feedback 
with  gains  selected  to  minimize  a  quadratic  performance  index.  This  approach  uses  a  lower 
order  state  estimator  in  the  feedback  loop  than  is  required  by  the  LQG  approach,  but  it  does 
not  provide  a  direct  method  for  achieving  desired  response  types  nor  does  it  result  in 
guaranteed  stability  margins. 

Although  LQR  and  LQG/LTR  theories  guarantee  stability  robustness,  these  theories  do 
not  provide  a  direct  means  to  obtain  the  desired  response  types.  Response  types  are 
obtained  indirectly  through  choice  of  state  and  control  weighting  matrices.  The  LQR 
method  requires  full  state  feedback  which  is  usually  not  feasible.  The  LQG/LTR  maintains 
stability  margins  without  full  state  feedback  but  requires  an  estimator  in  the  feedback  loop 
and  LQG  approach  to  LTR  can  result  in  closed  loop  poles  far  into  the  left  half  plane, 
resulting  in  high  compensator  bandwidth. 

H  Infinity  ( //<*>  )  Design  Techniques 

Recently  there  has  been  interest  in  the  application  of  Hoo  design  techniques  to 
helicopter  flight  control  systems.  The  H ^  approach  [54]  can  be  regarded  as  a  MIMO 
generalization  of  classical  SISO  frequency  response  techniques.  Since  the  desired  response 
types  for  helicopters  are  given  in  terms  of  frequency  response  relationships  between  pilot 
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inputs  and  helicopter  outputs,  the  H approach  appears  to  provide  a  natural  method  for 
design  of  helicopter  flight  control  laws. 


disturbances,  d 


The  Hoo  design  procedure  produces  a  stabilizing  control  law,  K(s),  which  fits  the 
feedback  structure  shown  in  Fig.  1.3.  The  //«,  norm  is  defined  for  a  transfer  matrix  as  the 

supremum  of  the  maximum  singular  value  of  the  matrix  taken  over  all  frequencies.  The 
Hoo  design  approach  is  based  on  the  the  selection  of  a  control  law  which  minimizes  the 

norm  of  transfer  function  matrices  associated  with  system  performance.  Thus  the  matrices 
selected  for  //«,  optimization  should  have  outputs  which  are  to  be  small.  In  the  case  of 

helicopters,  tracking  errors,  e,  actuator  inputs,  u,  and  responses  to  disturbances,  y,  are 
quantities  which  should  be  kept  as  small  as  possible  and  consequently  are  candidates  for 
Hoo  optimization.  From  Fig.  1.3,  the  relationship  between  the  command  input,  yc,  and  the 

error,  e,  is 

e  =  ( I  +  GK  )_1  yc 


also  the  relationship  between  the  output  and  disturbances  at  the  output  is 
y  =  ( I  +  GK  )->  d 
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The  matrix  (  I  +  GK  )-1  is  commonly  called  the  sensitivity  matrix.  The  transfer 
function  matrix  between  the  command  input,  yc,  and  the  actuator  input,  u,  is 

u  =  K  ( I  +  GK  V1  yc 

The  largest  singular  value  of  the  sensitivity  matrix  also  determines  the  size  of  the 
smallest  unstructured  output  inverse  multiplicative  perturbation  that  could  destabilize  the 
system  and  for  good  stability  margins  at  the  plant  output,  the  norm  of  the  sensitivity 
should  be  less  than  unity.  The  H norm  of  ( I  +  KG  )_1  determines  the  smallest  input 
inverse  perturbation  which  could  destabilize  the  system.  The  H ^  norm  of  the 

complementary  sensitivity  ( I  +  KG  l'1  KG  determines  the  smallest  unstructured  input 
multiplicative  perturbation  which  could  destabilize  the  system  and  the  //«>  norm  of  K  ( I  + 

GK  )_1  determines  the  largest  additive  unstructured  perturbation  which  could  destabilize  the 
system.  Thus  minimization  of  the  //<*,  norm  of  the  sensitivity  and  other  related  matrices  is 
important  for  purposes  of  stability  as  well  as  performance.  H M  design  procedures  tend  to 
give  nearly  constant  values  for  the  maximum  singular  values  of  the  matrices  which  are 
optimized.  Thus  in  practice  it  is  necessary  to  weight  these  matrices  with  appropriate 
frequency  dependent  weighting  matrices  in  order  to  obtain  desired  frequency  response 
characteristics.  Selection  of  these  weighting  matrices  is  one  of  the  most  important  tasks  for 
the  designer  using  the  H approach. 

Yue  and  Postlethwaite  [55]-[56]  have  described  the  application  of  //,*>  techniques  to 

the  determination  of  feedback  control  laws  for  a  mathematical  model  of  the  Westland  Lynx 

helicopter  at  hover.  They  were  interested  in  designing  control  laws  which  resulted  in 

ACA’i  responses  in  the  pitch  and  roll  axes,  and  first  order  yaw  rate  and  vertical  velocity 

command  responses.  They  had  to  include  as  outputs  pitch  and  roll  rates  ;i>  well  as  the 

commanded  variables,  pitch  and  roll  angles,  yaw  rate  and  vertical  velocin  m  order  to 
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achieve  acceptable  performance.  These  extra  two  outputs  increased  the  order  of  the  control 
laws.  Yue  and  Postlethwaite  chose  to  minimize  the  H ^  norm  of 

W ,  (  I  +  GK  )_1  W  3 
W2  K  (  I  +  GK  )-’  W3 

Since  (  I  +  C-K  )-1  is  related  to  tracking  errors  and  disturbance  rejection,  Wj  was 
selected  to  be  large  at  lower  frequencies  to  improve  tracking  and  disturbance  rejection. 
That  is  Wj  had  integral  like  characteristics  and  a  bandwidth  of  about  8  rad/s.  Since  ( I  + 
GK  )4  K  is  related  to  actuator  inputs,  W2  was  selected  to  be  small  at  low  frequencies  and 
large  at  higher  frequencies  to  penalize  high  frequency  inputs  to  the  actuators.  W2  was 
comprised  of  first  order  high  pass  filters  with  a  cut  off  frequency  of  10  rad/s  to  limit  system 
bandwidth.  W3  was  a  constant  matrix  selected  to  improve  tracking  performance. 

The  flight  control  laws  designed  by  Yue  and  Postlethwaite  using  H techniques  were 
tested  on  a  piloted  moving  base  simulator  at  the  Royal  Aircraft  Establishment  in  Bedford  in 
the  U.  K.  Hover  and  forward  flight  at  80  knots  were  simulated.  The  flight  characteristics 
at  these  two  conditions  are  very  dissimilar.  At  the  hover  design  condition,  pilot  ratings 
were  good.  At  80  knots  the  ratings  were  less  good  but  still  better  than  the  unaugmented 
helicopter. 

Walker  and  Postlethwaite  [57]  studied  two  //«,  techniques  for  the  design  of  flight 
control  laws  for  the  Lynx.  In  order  to  reduce  the  order  of  these  compensators,  they 
designed  fixed  gain,  inner  loop  controllers  for  stability  enhancement.  By  this  method,  they 
were  able  to  use  only  the  four  commanded  variables  as  outputs.  Their  first  design  used 
essentially  the  same  H0 0  norm  as  in  Ref  [28].  The  other  method  used  a  loop  shaping/ 

normalized  left  coprime  factor  robust  design  procedure  developed  by  McFarlane  and 
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Glover  [58]-[59].  The  use  of  only  four  outputs  gave  poorer  performance  than  obtained 
from  the  six  output  design  of  Yue  and  Postlethwaite. 

//oo  techniques  are  relatively  new.  Until  recently,  algorithms  to  calculate  control 

laws  were  primitive,  however,  new  computational  reliable  algorithms  have  been  developed 
[54].  The  underlying  mathematical  theory  of  //«>  design  is  sophisticated.  Selection  of 

weighting  functions  is  critical  to  the  success  of  the  design  and  is  still  an  art  which  may 
require  repeated  design  iterations.  H0 „  techniques  can  result  in  high  order  controllers,  and 

practical  implementation  often  requires  reduction  of  the  order  of  the  controller.  Since  fixed 
gain  controllers  probably  will  not  provide  adequate  handling  qualities  through  out  the  flight 
envelope  of  typical  high  performance  helicopters,  some  type  of  control  law  scheduling  will 
be  required  regardless  of  the  design  technique  used.  The  problem  of  scheduling  H x 

controllers  is  not  easy  as  there  is  no  restriction  on  the  control  law  structure  for  different 
flight  conditions,  i.e.  a  controller  designed  for  one  flight  condition  may  be  completely 
different  in  structure  from  that  designed  for  another  flight  condition.  Although  more 
powerful  computers  for  on  board  flight  control  may  allow  different  control  structures  to  be 
used  for  different  flight  conditions,  the  problem  of  smooth  transitions  between  modes  has 
not  yet  been  studied.  On  the  other  hand,  //<*, is  a  powerful  design  procedure  which  when 

coupled  with  structured  singular  value  ideas  allows  the  explicit  consideration  of  both 
nominal  performance  and  stability  and  performance  robustness  in  the  design  process.  The 
question  of  whether  this  power  is  required  in  the  design  of  helicopter  flight  control 
systems,  has  not  yet  been  resolved. 

Nonlinear  Dynamic  Inversion 

Nonlinear  dynamic  inversion  or  feedback  linearization  is  a  simple  idea  for  control  of 
nonlinear  systems.  It  consists  of  using  feedback  to  cancel  the  nonlinear  terms  and  then  the 
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design  of  the  linearized  system  can  be  accomplished  using  standard  linear  controller 
synthesis  procedures.  There  are  several  important  issues  which  can  affect  the  use  of  this 
method.  First,  exact  linearization  of  the  system  depends  on  exact  knowledge  of  the 
nonlinear  terms.  There  will  always  be  errors  in  the  mathematical  representation  of  the 
helicopter  and  the  control  system  must  be  robust  to  the  effects  of  these  errors.  Second, 
linearization  can  not  be  obtained  if  any  of  the  control  actuators  are  saturated.  Third,  if 
aerodynamic  data  is  available  from  look  up  tables  as  is  usually  the  case,  considerable 
computational  time  may  be  required  to  calculate  the  terms  in  the  nonlinear  inversion  control 
laws.  Despite  these  difficulties,  there  is  considerable  interest  in  this  method  and  a 
controller  using  a  nonlinear  inversion  controller  has  been  flown  by  NASA  on  UH-1H 
helicopter  [24], 

Eigenstructure  Assignment 

Eigenstructure  assignment  is  a  technique  for  synthesis  of  feedback  control  laws  which 
allows  the  designer  to  place  closed  loop  eigenvalues  in  specified  configurations  and  shape 
closed  loop  eigenvectors.  Eigenstructure  assignment  can  be  a  very  useful  technique  when 
performance  specifications  can  be  expressed  in  terms  of  closed  loop  eigenvalues  and 
eigenvectors  as  is  the  case  with  aircraft  handling  quality  specifications  and  many 
researchers  have  applied  this  technique  to  fixed  wing  aircraft  [60]-[67].  In  fact,  the  flight 
control  system  for  the  Aerospatiale  A320  Airbus  was  designed  using  eigenstructure 
assignment  [67],  Eigenstructure  techniques  provide  a  direct,  simple  method  for  achieving 
desired  nominal  performance  but  do  not  provide  any  guarantee  of  either  stability  or 
performance  robustness. 

Eigenstructure  assignment  is  also  useful  for  systems  which  are  not  completely 
controllable.  In  such  systems,  the  uncontrollable  eigenvalues  cannot  be  vi'.m-ed  by 
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feedback  but  the  associated  eigenvectors  can  be  altered  [68]  to  adjust  the  magnitude  of  the 
uncontrollable  mode  in  each  of  the  states. 

The  use  of  eigenstructure  assignment  for  helicopter  control  systems  has  been 
considered  by  several  authors.  Early  work  in  the  area  was  reported  by  Parry  and  Murray- 
Smith  [69],  Garrard  and  Liebst  [31],  and  Ekblad  [70],  These  papers  used  eigenstructure 
assignmenr  for  decoupling  and  stabilization  but  not  the  direct  achievement  of  response 
types.  Apkarian  [71]  also  applied  eigenstructure  assignment  to  a  helicopter  to  illustrate  his 
approach  to  robustness  enhancement  of  control  laws  designed  by  eigenstructure 
assignment.  Innocenti  and  Stanziola  [33]  have  described  performance/robustness  trade¬ 
offs  for  eigenstructure  control  laws  for  rotorcraft.  In  a  recent  paper,  Hughes,  Maness,  and 
Murray-Smith  [34]  used  eigenstructure  assignment  to  design  inner  loop  control  laws  which 
stabilize,  simplify  and  decouple  the  responses  to  control  inputs.  Simplification  is 
accomplished  by  selecting  eigenvector/eigenvalue  pairs  to  cancel  stable  invariant  zeros  and 
the  associated  left  zero  vectors.  This  is  analogous  to  pole-zero  cancellation  for  SISO 
systems.  Decoupling  and  stabilization  are  accomplished  by  suitable  choices  of  the 
remaining  eigenvalue/eigenvector  pairs.  Command  following  is  then  achieved  by  a  pre¬ 
compensator  designed  using  the  Broussard  Command  Generator  Approach.  An  ACAH 
system  is  synthesized.  Samablancat,  Apkarian,  and  Patton  [72]  use  a  combined 
eigenstructure  and  H ^  approach.  Inner  loop  controllers  which  decouple  the  response  into 

longitudinal,  lateral,  and  directional  modes  are  designed  using  eigenstructure  assignment. 
An  ACAH  response  type  is  then  achieved  using  the  H technique  in  a  manner  similar  to 

that  of  Yue  and  Postlelwaite.  The  objective  of  the  decoupling  is  to  reduce  the  order  of  the 
H0 o  compensators  without  loss  of  performance. 

Garrard,  Low,  and  Prouty  [32]  use  eigenstructure  assignment  to  design  a  feedback 

controller  which  achieves  desired  response  types  directly.  This  is  accomplished  by  using 

-21  - 


the  desired  response  types  to  generate  a  state  space  models  which  in  turn  generates  desired 
eigenvalue/eigenvector  sets.  Eigenstructure  assignment  is  then  used  to  design  feedback 
control  laws  which  give  the  best  least  squares  approximation  of  the  desired  eigenstructure. 
This  control  law  requires  full  state  feedback.  Since  only  angular  body  rates  and  vertical 
velocity  were  assumed  to  be  measured,  a  state  estimator  was  required  in  the  feedback  loop. 
LTR  was  achieved  by  an  eigenstructure  assignment  technique  due  to  Kazerooni  and  Houpt 
[73].  This  control  structure  gave  good  decoupling,  stability  robustness  and  the  desired 
time  and  frequency  responses  but  was  some  what  complex. 

This  study  represents  a  further  extension  of  the  work  of  Garrard,  Low,  and  Prouty. 
Eigenvalue  placement  is  used  to  design  an  inner  loop  control  law  which  provides 
decoupled,  rate  command  response  in  each  control  loop.  Eigenvalue  placement  is  used  to 
stabilize  unstable  modes  and  provide  appropriate  bandwidth,  and  eigenvector  shaping  is 
used  for  modal  decoupling.  With  the  inner  loops  closed,  the  helicopter  response  to  pilot 
commands  appears  as  four  decoupled,  first  order  systems.  It  is  then  possible  to  use 
classical  SISO  frequency  response  techniques  to  design  outer  loop  controllers  which  give  a 
variety  of  desired  response  types.  It  is  believed  that  the  integration  of  eigenstructure  and 
classical  techniques  presented  in  this  thesis  will  encourage  wider  use  of  eigenstructure 
design  techniques  in  helicopter  and  aircraft  flight  control  law  design. 

III.  Organization  of  the  Thesis 

The  remainder  of  this  thesis  is  organized  into  seven  chapters.  Chapter  II  reviews 
portions  of  the  new  specifications  on  handling  qualities  which  are  relevant  to  control 
system  design  for  hover  and  low  speed  flight.  Chapter  III  describes  the  theory  of 
eigenstructure  assignment.  Chapter  IV  discusses  the  mathematical  models  used  in  control 
system  design  and  evaluadon.  Chapter  V  covers  the  nominal  design  lor  hover  and  low 
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speed  flight-  Chapter  VI  addresses  the  topics  of  stability  and  performance  robustness  using 
a  variety  of  methods.  Chapter  VII  describes  the  conclusions  and  recommendations  for 
further  research. 
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Chapter  2 


Helicopter  Handling  Qualities 


I.  Introduction 

Recent  advances  in  rotorcraft  technology  have  expanded  the  role  of  the  modem 
helicopter.  It  is  desirable  that  the  handling  qualities  of  modern  helicopters  reflect  the 
demands  of  particular  mission  environments.  The  distinctive  characteristics  of  particular 
mission  environments  make  handling  qualities  a  challenge  for  the  control  designers.  The 
current  helicopter  handling  qualities  specification,  MIL-H-8501A  [1],  written  in  1952  and 
mainly  unchanged  since  then,  does  not  address  specific  handling  qualities  issues  relevant  to 
modem  helicopters.  The  U.S.  Army  is  currently  engaged  in  revising  and  updating  this 
document  to  reflect  advances  in  rotorcraft  technology  since  1952  and  link  it  to  the  design 
mission  environments  of  helicopters.  The  latest  version  [2]  has  been  adopted  by  AVSCOM 
as  an  Airworthiness  Design  Standard  ADS  33  (  May  1988  )  and  used  for  designing  the 
U.S.  Army  Light  Helicopter  Family  (  LH  )  in  June  1988.  This  new  handling  quality 
specification  also  appears  to  be  gaining  acceptance  for  civilian  use. 

Selected  features  and  criteria  are  proposed  for  design  guidance  and  the  evaluation  of 
handling  qualities.  The  primary  feature  of  the  proposed  method  is  the  classification  of 
response  types  required  to  complete  the  specified  missions  in  the  presence  of  specified 
visual  environments  and  displays.  In  addition,  separate  criteria  are  suggested  for  small, 
moderate,  and  large  amplitude  maneuvering.  Frequency  domain  criteria  are  suggested  for 
small  amplitude,  precision,  closed  loop  control,  and  time  domain  criteria  are  provided  for 
low  frequency,  and /  or  large  amplitude  maneuvering.  The  new  handling  quality 
specifications  [2]  are  very  lengthy  and  only  those  specifications  which  direcil;.  i elate  to 
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flight  control  system  design  for  hovering  and  low  speed  flight  are  outlined  in  this  Chapter. 
For  higher  speed  flieht,  the  handling  qualities  specifications  are  similar  to  fixed  wing 
aircraft. 

II.  Specification  Overview 

For  application  in  the  specification,  the  requirements  have  to  be  translated  into  specific 
Mission  Task  Elements  (  MTE  )  and  an  Operational  Flight  Envelope  (  OFE  ).  The  MTEs 
define  which  requirements  need  to  be  satisfied  and  at  what  level  of  performance.  The  OFE 
defines  the  envelope  of  speed,  normal  load  factor,  altitude,  rate  of  climb,  sideslip  angle  and 
other  parameters  that  are  of  operational  significance,  within  which  the  rotorcraft  is  capable 
of  successfully  performing  its  mission.  Within  this  OFE,  the  rotorcraft  must  achieve  Level 
1  handling  qualities. 

The  Usable  Cue  Environment  (  UCE  )  provided  to  a  pilot  is  determined  'v  combining 
information  from  the  operational  environment  with  vision  aids  and  displays.  The 
applicable  MTE  and  UCE  are  then  combined  to  determine  the  required  response  types  and 
the  character  of  these  responses. 

The  requirements  which  have  to  be  met  by  a  helicopter  are  dependent  on  several 
parameters:  Mission  Task  Element  ( MTE  ),  Usable  Cue  Environment  (  UCE  ),  pilot 
attention  level,  and  speed  range.  These  are  be  discussed  below. 

Mission  Task  Element  (  MTE  ):  Rotorcraft  performing  different  tasks  may 
require  different  handling  qualities.  Therefore,  the  rotorcraft  specification  defines  a  set  of 
maneuvers  which  represent  the  essential  ingredients  of  most  missions;  these  are  termed 
Mission  Task  Elements.  MTEs  can  further  be  divided  into  stationary  and  translating  MTEs 
for  hover  and  low  speed  flight. 
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Usable  Cue  Environment  (  UCE  ):  The  visual  cues  available  to  a  pilot  are 
compiled  to  produce  a  Usable  Cue  Environment  (UCE).  These  visual  cues  include 
whatever  the  pilot  can  see  out  of  the  window  with  or  without  vision  aids.  Three  levels  of 
UCE  are  defined,  with  level  1  being  the  best  and  providing  cues  for  aggressive  and/or 
precise  maneuvering,  and  level  3  being  the  worst  and  providing  cues  only  for  small  and/or 
gentle  maneuvers. 

Pilot  Attention  :  The  two  levels  of  response  to  controls  and  disturbances  which  are 
defined  in  the  specification  are  the  two  levels  of  pilot  attention.  Fully  attended  tasks  are 
performed  when  pilots  apply  full  attention  and  are  allowed  more  lenient  standards  then 
partially  attended  tasks,  which  required  the  pilot  to  engage  in  other  activities  in  addition  to 
flight  path  management  tasks. 

Speed  Range  :  There  is  a  major  break  in  the  requirements  at  the  speed  of  45  knots. 
For  speeds  less  than  45  knots,  (  hover  and  low  speed  ),  maneuvers  are  ground  referenced, 
and  these  rotorcraft  dynamics  are  assumed  to  be  typical  of  a  hovering  helicopter  with  pitch 
and  roll  having  similar  dynamics,  and  yaw  and  vertical  velocity  ideally  having  uncoupled 
first  order  responses.  Above  45  knots,  (  forward  flight ),  rotorcraft  dynamics  are  assumed 
to  approach  those  of  fixed-wing  aircraft. 

III.  Specification  Criteria 

Basic  stability  and  control  requirements  are  divided  into  two  categories:  hover  and  low 
speed,  and  forward  flight.  Due  to  the  symmetry  of  the  xz  plane  of  a  helicopter,  pitch  and 
roll  have  similar  requirements.  However,  yaw  and  vertical  velocity  responses  have 
separate  requirements. 

A .  Bandwidth  Criterion 
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Bandwidth  describes  the  ability  to  follow  a  range  of  input  frequencies.  It  is  defined  as 
the  frequency  at  which  the  output  magnitude  is  3  dB  less  than  the  input  magnitude  (  output 
magnitude  is  0.707  of  the  input  magnitude  ).  A  good  system  normally  has  a  high 
bandwidth  relative  to  the  maximum  input  frequency  that  it  is  designed  to  track.  Bandwidth 
limits  are  set  by  system  stability  considerations. 

The  bandwidth  is  referenced  to  the  aircraft  as  the  maximum  crossover  frequency  with 
all  augmentation  loops  closed.  Additional  outer  loops,  ( including  the  pilot )  can  be  closed 
later.  The  frequency  response  data  required  to  measure  the  bandwidth  parameters  include 
all  the  elements  of  the  analog  or  digital  flight  control  system,  (  e.g.  anti-aliasing,  bending 
mode  and  stick  filters,  zero  order  hold  effects,  actuators,  sensors  and  computational 
delays  ). 


FREQUENCE  (rad/s) 

Fig.  2. 1  Definition  of  Gain  Margin  Limited  Bandwidth. 


As  shown  in  Figures  2.1  and  2.2,  there  are  two  definitions  of  bandwidth  frequency, 
namely  gain  margin,  Wgwgain’  phase  margin,  coswphase-  These  two  definitions  of 
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bandwidth  frequency  are  defined  for  a  second  ordcr  system.  These  figures  describe  the 
margin  above  the  augmented  vehicle's  response  in  which  the  pilot  can  increase  the  gain  or 
add  a  time  delay  or  phase  lag  without  causing  any  instability  to  the  design  model. 


FREQUENCY  <rad/s> 


Fig.  2.2  Definition  of  Phase  Margin  Limited  Bandwidth. 

The  gain  margin  bandwidth  is  defined  as  the  frequency  for  6  dB  of  gain  margin.  To 
ascertain  the  gain  margin  bandwidth,  the  gain  corresponding  to  0)]80  ( the  frequency  at 
which  the  phase  shift  is  -180° )  must  be  determined.  The  frequency  at  which  the  magnitude 
of  the  gain  is  6  dB  above  the  magnitude  at  co180  equals  the  gain  margin  bandwidth.  Low 
gain  margin  bandwidths  tend  to  result  in  system  configurations  which  are  PIO  (  pilot 
induced  oscillation,  ~1  Hz  )  prone  because  small  changes  in  the  pilot  gain  result  in  a  rapid 
reduction  in  phase  margin.  For  phase  margin  bandwidth,  the  frequency  for  co180  can  be 
determined  from  the  Bode  phase  plot  (  Figure  2.2  ).  The  bandwidth  for  phase  margin 
systems  is  the  frequency  at  which  the  phase  margin  is  45  degrees.  The  definition  of 
bandwidth  for  ACAH  Response  -  Types  does  not  include  a  gain  margin  limit.  Also,  gain 

margin  limited  ACAH  Response  -  Types  lead  to  PIO  for  super  precFir>n  tasks. 
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Nevertheless,  gain  bandwidth  can  be  used  as  a  requirement  for  Rate  or  RCAH  Response  - 
Types,  where  the  pilot  attitude  loop  closure  is  necessary  to  maintain  stability  at  hover. 

For  first  order  systems,  the  bandwidth  and  the  eigenvalue  Xn,  are  identical.  The 
minimum  bandwidth  [3]  that  will  satisfy  Level  1  small  amplitude  yaw  rate  or  vertical 
velocity  requirements  for  minimum  gain  tasks  is  2  rad/s. 

For  a  second  order  system,  Eq.  2. 1  from  Ref.  4  can  be  used  to  correlate  the  bandwidth 
requ.  'ments  given  by  the  handling  quality  specification  to  an  appropriate  set  of  time 
domain  parameters  which  includes  the  natural  frequency,  con  and  damping  ratio  £,  of  the 
system 

®bw  =  CDnU  +  V  C2  +  1  )  (2.1) 

The  typical  values  for  coBW,  0)n,  and  £  that  will  satisfy  Level  1  handling  quality 
requirements  for  small  amplitude  pitch  (  roll  )  attitude  changes  are  coBW  =  4  rad/s, 
con  =  2.071  rad/s  and  £  =  0.707. 

Pilots  are  sensitive  to  whether  the  phase  curve  drops  off  rapidly  or  slowly  at 
frequencies  beyond  the  bandwidth  frequency.  This  is  characterized  by  the  phase  delay 
parameter.  Phase  delay  is  a  measure  of  the  steepness  at  which  the  phase  drops  off  after 
-180°  and  is  an  indication  of  the  behavior  of  the  rotorcraft  as  the  pilot  increases  the 
crossover  frequency.  Large  phase  delay  means  that  there  is  a  small  margin  (  range  of 
frequencies )  between  normal  tracking  at  45  degrees  of  phase  margin  and  instability.  In 
addition,  rotorcraft  with  large  phase  delay  are  PIO  prone.  Phase  delay  is  defined  as  a  two- 
point  approximation  to  the  phase  curve  defined  at  co180  and  2co180  as  in  Eq.  2.2. 

T  _  _ -_^2wl8Q _  (j  2) 

p  57.3  (2colgo) 
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where  A<t>2coi80  *s  the  phase  difference  between  -180°  and  phase  at  2g)180  ( the  frequency 
that  is  two  times  the  frequency  when  the  phase  angle  is  -180°  ).  If  phase  is  nonlinear 
between  colgo  and  2co180,  xp  is  determined  from  a  linear  least  squares  fit  to  the  phase  curve 
between  co180  and  2co180.  Hoh  [5]  has  shown  that  for  two  systems  which  have  essentially 
the  same  bandwidth,  the  system  that  has  a  significantly  higher  phase  delay  has  significantly 
degraded  pilot  ratings.  This  is  due  to  the  fact  that  phase  delay  is  a  measure  of  the  shape  of 
the  phase  curve  around  the  instability  frequency.  In  order  to  achieve  Level  1  handling 
quality,  it  is  desired  to  have  a  value  of  0.2  seconds  or  less  for  Xp. 

B  .  Response  -  Types 

As  discussed  briefly  in  Chapter  I,  three  Response  -  Types  have  been  developed  to 
quantify  mission  oriented  rotorcraft  handling  quality  specifications  [5]-[6].  These  are  as 
follows  : 

1 .  Attitude  Command  with  Attitude  Hold  (  ACAH  ) 

2 .  Rate  Command  ( RC )  or  Rate  Command  with  Attitude  Hold  (  RCAH ) 

3 .  Translational  Rate  Command  With  Position  Hold  ( TRCPH ) 

In  ACAH,  a  constant  control  input  must  result  in  a  proportional  attitude  displacement 
in  the  presence  of  external  disturbances.  For  a  pulse  control  input,  the  attitude  must  return 
to  within  10%  of  its  peak  value  within  20  seconds.  This  response  type  is  required  for 
hover  and  low  speed  operations  in  conditions  of  degraded  visual  cueing  and  for  divided 
attention  tasks. 

In  RCAH  systems,  attitude  must  diverge  away  from  trim  for  at  least  4  seconds 

following  a  step  command.  This  response  type  is  preferred  for  most  MTEs  for  fully 
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attended  operations  in  conditions  of  good  visual  cueing.  For  TRCPH,  a  constant  controller 
force  input  must  produce  a  constant  translational  rate  and  the  rotorcraft  must  hold  position 
if  the  force  on  the  cockpit  controller  is  zero.  TRCPH  is  necessary  to  achieve  Level  1 
handling  quality  in  nap  of  the  earth  maneuvers  in  fair  to  poor  visual  cue  environments. 


The  desired  roll  attitude  and  roll  rate  transfer  functions  for  ACAH  and  RCAH 
respectively  are  [3],  [5],  [7] 

C02a 


1 

4>c 


s2  +  2CCO0S  +  co2< 

2.  _  (  TS  +  1  )  (D2^ 

Pc  S2  +  2£<D6S  +  co  2d 


(2.3) 


(2.4) 


Transfer  functions  similar  to  Eqs.  2.3  and  2.4  are  required  for  the  pitch  axis. 


The  desired  vertical  velocity  and  yaw  rate  transfer  functions  for  RC  are 

_w  _  *w 
wc  s  +  Xw 

L  = 

rc  s  +  A,r 


(2.5) 

(2.6) 


The  vertical  velocity  and  yaw  rate  bandwidths  are  given  by  and  \  ,  respectively.  In 
order  to  satisfy  ultra  high  gain  tasks,  the  values  for  these  eigenvalues  are  typically  4  rad/s. 
For  ordinary  flight  path  management,  can  be  much  less. 


TRCPH  can  be  achieved  by  second  order  transfer  functions  of  the  form 

_u.  _  _ Wju _ 

uc  s2  +  2£(ous  +  (02u 


(2.7) 
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V 


(2.8) 


o £v _ 

vc  s2  +  2£covs  +  co  2  v 

In  general,  bandwidths  of  the  order  of  1  rad/s  are  required  for  TRCPH.  This  is 
considerably  less  than  the  4  rad/s  required  for  RC  and  ACAH. 

C .  Dynamic  Requirements,  Hover  and  Low  Speed 

Dynamic  response  to  control  requirements  are  divided  into  small,  medium  and  large 
amplitude  responses.  The  small  amplitude  response  is  characterized  by  attitude  and 
velocity  stability,  or  low  pilot  workload  in  the  primary  control  axis.  This  mainly  involves 
closed  loop  tracking  tasks  and  can  also  be  characterized  by  bandwidth  frequency. 

1 .  Small  Amplitude  Attitude  Changes  (  Pitch  and  Roll  ) 

The  small  amplitude  response  criteria  are  divided  into  high  frequency  and  mid-to-low 
frequency  spectra.  Small  amplitude  criteria  can  be  evaluated  using  linear  models.  The  high 
frequency  criteria  are  specified  in  terms  of  the  maximum  closed  loop  bandwidth  parameter. 
The  bandwidth  is  defined  as  a  function  of  phase  or  gain  margin  of  the  open  loop  Bode  plot 
as  shown  in  Figures  2.1  and  2.2. 

It  can  be  seen  from  Fig.  2.3  that  the  required  bandwidth  increases  as  the  phase  delay 
parameter  increases  above  approximately  0.17  seconds.  Fig.  2.3  also  shows  that  the 
bandwidth  is  very  sensitive  to  lags  and  delays  in  the  flight  control  system.  Phase  delay  is  a 
measure  of  how  quickly  the  phase  lag  increases  at  frequencies  beyond  that  for  which  phase 
shift  is  - 1 80  degrees.  A  rotorcraft  with  rapid  phase  drop  off  will  be  more  sensitive  than 
one  with  gradual  phase  drop  off.  Therefore  a  higher  bandwidth  is  required  to  maintain 
Level  1  handling  qualities. 
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u>BW,  ra  d/s 

fll I  Other  MTEs  -fill  UCE 

Fig.  2.3  Requirements  for  Small  Amplitude  Attitude  Changes. 

Helicopters  exhibit  relatively  low  frequency  responses  in  side  and  forward  velocity 
which  are  somewhat  analogous  to  the  phugoid  response  of  fixed  wing  aircraft.  These  are 
called  the  mid-term  responses.  The  handling  qualities  for  these  mid-term  responses  are 
specified  in  terms  of  the  frequency  and  damping  ratio.  For  Level  1  handling  qualities,  a 
damping  ratio  of  at  least  0.35  is  required  except  when  frequencies  fall  below  0.5  rad/s.  For 
frequencies  below  0.5  rad/s,  unstable  responses  are  acceptable  provided  the  time  to  double 
is  not  too  large. 

The  mid-term  response  is  specified  in  terms  of  frequency  and  damping  of  characteristic 
oscillations  as  shown  in  Fig.  2.4.  Mid-term  is  defined  as  a  frequency  range  below 
bandwidth  but  above  steady  state.  Mid-term  response  characteristics  represent  the  handling 
qualities  most  affected  by  mission  requirements  for  divided  attention.  Simply  stated,  the 
pilot  should  be  able  to  relinquish  control  of  the  helicopter  for  short  periods  without 
encountering  significant  excursions.  Fig.  2.4  represents  a  more  stringent  Level  1  boundary 
and  less  stringent  Level  2  boundary  [8]  than  previously  used.  Since  the  old  Mil  H-8501A 
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[1]  boundary  allowed  essentially  zero  damping  at  ton  =  1  rad/s,  the  revised  Level  1 
boundary  is  more  consistent  with  the  proposed  bandwidth  requirement  for  cogw  >  1  rad/s . 


Fig.  2.4  Limits  on  Mid-Term  Eigenvalues. 


For  Level  1  and  divided  attention  tasks,  a  damping  ratio  of  at  least  0.35  is  required  for 
all  mid  and  low  frequency  oscillations.  If  full  attention  is  possible,  low  frequency 
instabilities  at  frequencies  below  0.5  rad/s  are  allowed  for  Level  1.  A  minimum  damping 
ratio  of  -0.2  is  required  for  these  low  frequency  unstable  eigenvalues.  In  terms  of  doubling 
time,  Level  1  requires  tjou^  >  6.76  seconds.  As  can  be  seen  from  Fig.  2.4,  Level  2 
requires  low  frequency  unstable  oscillatory  modes  to  have  a  maximum  bandwidth  of  0. 1 39 
rad/s.  This  corresponds  to  a  minimum  allowable  t,jouble  of  4.96  seconds. 
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2 .  Moderate  Amplitude  Attitude  Changes  (  Pitch  and  Roll  ) 


nil  Other  MTEs  -  fill  UCE  <  roll  ) 

Fig.  2.5  Requirements  for  Moderate  Amplitude  Roll  Attitude  Changes. 

Evaluation  of  moderate  amplitude  attitude  changes  require  nonlinear  analyses  and 
simulations.  The  requirement  for  moderate  attitude  changes  are  shown  in  Fig.  2.5.  The 
parameter  ppk/Ad>pk  is  the  ratio  of  the  maximum  roll  rate  achieved  while  changing  from  one 
steady  roll  attitude  to  another,  divided  by  the  roll  attitude  change.  Moderate  amplitude 
specifications  given  in  terms  of  peak  angular  rate  per  attitude  change  are  the  same  for  hover/ 
low  speed  and  forward  flight. 

3 .  Large  Amplitude  Attitude  Changes  (  Pitch  and  Roll  ) 

The  large  amplitude  requirements  are  specified  as  the  minimum  achievable  angular  rate 
for  rate  response  types  and  as  the  minimum  achievable  angle  for  attitude  response  types. 
As  for  moderate  amplitude  changes,  evaluation  of  these  criteria  require  nonlinear  analyses 
and  simulations.  The  criteria  are  grouped  into  three  levels  of  aggressiveness  corresponding 
to  the  needs  of  the  Mission  Task  Elements,  as  shown  in  Table  2.1. 
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Rate  Response  Types 


Attitude  Response  Types 


Mission  Task  Element 

Minimum  Achievable 
Angular  Rate,  deg/sec 

Minimum  Achievable  Angle, 
deg 

Level  1 

Level  1 

P 

q 

r 

<!> 

e 

Limited  Maneuvering 

+  21 

±  6 

±  9.5 

±  13 

±  13 

Moderate  Maneuvering 

±50 

±  13 

±22 

±60 

+  20 
-30 

Aggressive  Maneuvering 

±50 

±  30 

±60 

±60 

±30 

Table  2.1  Requirements  for  Large- Amplitude  Attitude  Changes. 


4 .  Inter-axis  Coupling 


Vaiu  due  to  Col  I ec  t i ue  Coup  I  i  ng  Limits 


Fig.  2.6  Yaw  Cross  Coupling  Criteria. 


Cross-coupling  can  be  a  significant  problem  during  aggressive  maneuvering.  Fig.  2.6 
illustrates  the  requirement  that  collective  control  is  to  achieve  vertical  velocity  response  and 
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shall  not  cause  objectional  response  in  yaw  or  other  axes.  Limits  on  this  form  of  cross¬ 
coupling  are  specified  in  terms  of  yaw  rate  oscillation  parameters  per  rate  of  climb. 


From  Fig.  2.6,  the  vertical  axis  defines  the  absolute  ratio  of  yaw  rate  response  at  1 
second  to  vertical  velocity  response  at  3  seconds  due  to  a  step  collective  command. 
Likewise  the  horizontal  axis  defines  the  absolute  ratio  of  yaw  rate  response  at  3  seconds  to 


vertical  velocity  response  at  3  seconds  due  to  a  step  collective  command.  For  Level  1 


handling  qualities, 


must  be  smaller  than  0.65  and 


must  be  between  0.15  and  0.2. 


In  addition,  specific  requirements  are  provided  for  roll  due  to  pitch,  and  pitch  due  to 
roll.  These  requirements  apply  for  aggressive  maneuvering  MTEs.  Limits  are  placed  on 
the  ratio  of  off-axis  response  to  the  on-axis  response.  For  at  least  5  seconds  of  time 
response,  the  ratio  of  pitch  rate  to  roll  rate  for  a  lateral  cyclic  step  input  ( or  roll  rate  to  pitch 
rate  for  a  longitudinal  cyclic  step  input  )  shall  be  less  than  0.25  for  Level  1  handling 
qualities. 

5 .  Response  to  Disturbances 

To  ensure  adequate  disturbance  rejection,  the  specification  requires  short-term 
response  to  wind  gusts  to  have  the  same  bandwidth  as  the  response  to  control.  The 
specification  further  requires  inputs  to  be  made  directly  at  the  control  actuator  (  the 
swashplate  on  the  helicopter  )  so  that  the  resulting  response  can  meet  the  same  bandwidth 
criterion  as  response  to  control.  The  disturbance  rejection  requirements  only  apply  to  Level 
1  handling  qualities. 
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Chapter  3 


Theory  of  Eigenstructure  Assignment 

I.  Introduction 

For  a  controllable  system,  eigenstructure  assignment  techniques  use  feedback  control 
to  place  closed  loop  eigenvalues  and  shape  closed  loop  eigenvectors  to  achieve  performance 
specifications.  Eigenstructure  assignment  can  be  achieved  directly  by  least  squares 
techniques  [l]-[7],  or  asymptotically  by  linear  quadratic  regulator  (  LQR  )  methods  [8]-[9] 
Both  methods  have  advantages  and  disadvantages.  For  a  system  with  full  state  LQR 
controllers,  the  minimum  singular  value  of  the  return  difference  matrix  is  always  greater 
than  or  equal  to  one.  This  guarantees  gain  margins  of  at  least  ±6  dB  with  no  variations  in 
phase,  and  phase  margins  of  ±60  degrees  with  no  variations  in  gain  [  1 0]-[  1 1];  however, 
since  a  state  estimator  is  almost  always  required  in  the  feedback  loop,  the  stability  margins 
may  be  lost  in  practice  (  as  shown  by  Doyle  in  [121  )  unless  some  type  of  loop  transfer 
recovery  (  LTR  )  [  1 3]-[  14]  procedure  is  used  to  partially  restore  them. 

The  direct  approach  does  not  result  in  guaranteed  stability  margins;  however  in 
practice,  stability  margins  are  often  good  and  Gilbert  [15]  has  shown  that  sensitivity  of 
eigenvalue  locations  to  parameter  variations  can  be  minimized  if  the  closed  loop 
eigenvectors  are  made  as  orthogonal  as  possible.  The  asymptotic  LQR  approach  requires  a 
controllable  system  while  the  direct  approach  does  not.  In  fact,  the  direct  approach  can  be 
used  to  shape  eigenvectors  associated  with  uncontrollable  eigenvalues.  The  direct  approach 
can  be  used  to  design  a  gain  matrix  which  uses  sensor  outputs  for  direct  output  feedback 
although  in  this  case  stability  is  not  guaranteed. 
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Finally,  in  the  limit,  the  asymptotic  LQR  approach  drives  some  of  the  closed  loop 
eigenvalues  to  infinity.  In  practice,  this  may  result  in  excessively  high  bandwidth.  The 
direct  approach  can  be  used  to  maintain  eigenvalues  in  their  open  loop  position  if  the 
response  associated  with  these  eigenvalues  is  satisfactory  and  it  is  not  desired  to  move 
them.  When  compared  with  eigenstructure  assignment  by  LQR  theory,  the  advantages  of 
the  direct  method  outweigh  its  disadvantages,  so  the  direct  method  will  be  used  in  this 
study.  In  addition,  an  eigenstructure  method  for  the  design  of  state  estimators  which 
exhibit  LTR  has  been  developed  [16]. 

II.  Full  State  Feedback 

Consider  a  linear  time  invariant  system  described  by  the  equations 

x  =  A  x  +  B  u  (3.1) 

y  =  Cx  (3.2) 

where  x  e  SHn,  u  e  9im,  ye  SRr  are  called  the  state,  control  input  and  output  respectively. 
A,  B  and  C  are  real  constant  matrices  of  appropriate  dimensions  with  the  assumptions  that 

rank  [  B  ]  =  m  *  0 

and  rank  [  C  ]  =  r  *  0 

If  the  system  described  in  Eq.  3.1  and  3.2  is  controllable  and  observable  and  a  linear 
control  law 

u  =  -  K  x  (3.3) 

is  assumed,  the  following  results  can  be  proven  [2],  [4] : 
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(1)  The  position  of  max(m,r)  closed  loop  eigenvalues  can  be  assigned  arbitrarily 
provided  that  {  Vj*  },  i  =  1,  2,  n  is  a  self-conjugate  set  of  scalars. 

(2)  The  shape  of  max(mj)  eigenvectors  can  be  partially  assigned  with  min(m,r)  entries  in 
each  vector  arbitrarily  chosen.  {  v?  },  i  =  1,  2,  n  is  a  self-conjugau  set  of  n- 

eigenvectors  (  linearly  independent  )  corresponding  to  the  self-conjugate  set  of 
eigenvalues  {  }. 

(3)  The  eigenvector  must  reside  in  the  m-dimensional  subspace  spanned  by  the  columns 
of  ( Aql -  A  j'1  B. 

The  sets  {  A.^  }  and  {  v  ^  }  correspond  to  the  desired  closed  loop  eigenvalues  and 
eigenvectors  pair  while  the  matrix  K  is  the  constant  feedback  gam  matrix. 

The  solution  of  the  control  problem  is  a  real  ( mxn  )  matrix  K  such  that  the  eigenvalues 
of  A  -  BK  are  precisely  those  of  V?  with  corresponding  eigenvector  V|.  Since  a  fixed 

eigenstructure  uniquely  defines  ( A  -  BK  ),  it  is  true  that  if  K  exist  and  rank  [  B  ]  =  m,  then 
K  is  unique.  The  number  of  control  variables  available  determines  the  dimensions  of  the 
subspace  which  the  achievable  eigenvectors  must  reside.  The  orientation  of  the  subspace  is 
determined  by  the  matrices  A,  B  and  the  desired  closed  loop  eigenvalues  A,^.  In  general, 

more  than  min(m,r)  entries  in  a  particular  vector  have  to  be  assigned  in  order  to  shape  the 
transient  response  characteristics  of  the  system.  This  may  result  in  a  desired  eigenvector 
that  does  not  reside  in  the  subspace  spanned  by  the  columns  of  (  Aql  -  A  )']  B . 

Consequently,  the  desired  closed  loop  eigenvector  that  does  not  lie  precisely  within  the 

prescribed  subspace  can  not  be  achieved  exactly  and  can  only  be  achieved  in  a  least  squares 
sense  ( the  orthogonal  projection  of  v  ^  onto  the  subspace  spanned  by  the  columns  of  (  Aql  - 

A  )_1  B  ). 
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Consider  the  closed  loop  system  given  in  Eq.  3.4 


x  =  (  A  -  BK)  x  (3.4) 

If  it  is  desired  to  move  as  eigenvalue,  the  design  procedure  consists  of  determining  the  gain 
matrix  K  such  that  for  all  desired  closed  loop  eigenvalue/eigenvector  pairs  ( and  v ; ), 


(  A  -  BK  )  Vj  =  A-i  v £  (3. 5. a) 

or  vj  =  (Xjl-  A^BKvj  (3.5.b) 

We  define  the  m-vector  m  ,  as 

mi  =  -  Kvj  (3.6) 

This  is  equivalent  to  finding  m  j  such  that 

(Xjl  - A)  vj  =  B m{  (3.7) 

Once  the  m  j's  have  been  found,  the  gain  matrix  is  calculated  as 
K  =  -  MV'1  (3.8) 


where  M  is  a  matrix  whose  columns  are  m  j  and  V  is  a  matrix  whose  columns  are  v  j.  This 
is  the  feedback  gain  when  the  desired  eigenvector  resides  in  the  prescribed  subspace  and 
can  be  attained  exacdy. 

In  general,  the  desired  eigenvectors  are  not  achievable,  therefore,  the  m(s  are  selected 
to  minimize  the  weighted  sum  of  the  squares  of  the  difference  between  the  elements  of  the 
desired  and  attainable  eigenvectors.  This  is  expressed  as  the  following  performance  index 

Ji  =  (  vf-vf)*Pj(  vf-vj)  (3.9) 
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where  is  the  achievable  closed  loop  eigenvector,  and  Pj  is  a  positive  definite  symmetric 

matrix  whose  elements  are  coefficients  that  weight  certain  components  of  the  corresponding 
desired  and  achievable  eigenvectors  more  heavily  than  others.  For  a  fixed  it  is  desired 
to  minimize  Jj  subject  to  Eq.  3.8.  Adjoining  Eq.  3.8  to  Jj  with  the  Lagrange  multipliers  jij 

Ji  =  (vf-vffPiCvf-vfj  +  H^KXil-A)  vf-BmJ  (3.10) 


Now,  Eq.  3.10  is  minimized  when 


3J 

dm  i 


=  -  (Ij*  B  =  0 


(3.11) 


3J_ 

a 


=  (  vj.  vf  )*  Pi  +  M-i*  (  >.il  -  A  )  =  0 


(3.12) 


In  matrix  form.  Equations  3.6,  3.9,  and  3.10  become 


— 1 

>* 

> 

-B 

0 

“  vr 

-  o  - 

0 

0 

-BT 

= 

0 

Pi 

0 

>> 

> 

* 

1 _ 

-  P.  - 

-  p,''id . 

r  va  - 

-  0 

i 

Ni 

m, 

= 

0 

-  P.  - 

.  PjV,d  . 

Garrard  and  Liebst  (17]  have  shown  that  if  the  system  is  controllable,  Nj  is  non¬ 
singular  even  if  ^  remains  at  its  open  loop  position.  The  achievable  v,  and  mx  are  then 
given  in  Eq.  3.13. 
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(3.13) 


It  is  obvious  that  by  setting  ml  =  0  assures  that  the  associated  Vj  and  remain  in  their 
open  loop  configuration. 

III.  Full  State  Estimator  for  LTR 

Once  the  full  state  feedback  gain  matrix,  K,  has  been  determined,  the  problem  of 
constructing  a  state  estimate,  x,  from  measurements,  y,  must  be  addressed.  A  common 
regulator/estimator  scheme  is  as  follows.  Considering  the  control  law  given  by 

u  =  -  Kx  (3.14) 

where  £  is  the  estimated  state  vector  given  by  the  state  estimator's  equation 

x  =  Ax  +  Bu  +  L(y-Cx)  (3.15) 

Using  the  estimates  of  the  state  in  the  control  law  results  in  the  compensator  transfer 
function 

K(s)  =  K  (  si  -  A  +  BK  +  LC  I'1  L  (s)  (3.16) 

The  objective  is  to  choose  the  estimator  gain  matrix,  L,  such  that  (  A  -  LC  )  is  stable 
and  the  performance  of  the  regulator/estimator  closely  resembles  that  of  the  full  state 
feedback  for  some  bounded. frequency  range. 

In  the  full  state  feedback  case,  the  return  difference  matrix  is 


I  +  K  (  si  -  A  )_1  B  (s) 
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(3.17) 


which  in  the  combined  regulator/  estimator  case  is 

I  +  K(s)  G(s)  (3.18) 

where  K(s)G(s)  is  the  properties  when  the  loop  is  broken  at  the  input  of  the  plant.  It  is 
obvious  that  in  the  SISO  case  if  one  wishes  to  recover  F(s),  then  L(s)  should  be  chosen 
such  that 

K(s)  -  (3.19) 

inverting  the  plant  from  the  left  by  moving  the  poles  of  K(s)  towards  the  zeroes  of  G(s). 
Hence,  it  is  important  that  the  plant  is  minimum  phase  (  no  right  half  plane  zero  )  to  ensure 
that  the  eigenvalues  of  the  estimator,  (  A  -  LC  ),  are  stable. 

A  technique  for  designing  the  filter  gain  matrix,  L  in  the  MIMO  case  to  recover  F(s)  by 
eigenstructure  assignment  was  shown  by  Kazerooni  and  Houpt  [16].  As  shown  by  Doyle 
and  Stein  in  [13]  and  [18],  L  is  chosen  such  that  the  filter  gain  has  the  following  asymptotic 
behavior  as  q  — >  oo 

»  BW  where  WWT  =  I  (3.20) 

4 

When  this  L  is  used  in  the  loop  transfer  expression,  as  q  becomes  large  (  assuming  G(s)  is 
minimum  phase  ),  the  filter  gain  behaves  in  such  a  way  that 

^limJK(s)  G(s)  — >  K  (  si  -  A  )_1  B(s)  (3.21) 

where  the  convergence  is  pointwise  in  s.  This  design  procedure  essentially  inverts  the 
plant  from  the  left  i.e. 

lim  K(s)  =  K  ( si  -  A  )  ’  B  [  C  ( si  -  A  )-'  B  ]  ‘(s)  (3  22) 

q— >oo 
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This  inversion  is  why  G(s)  is  not  allowed  to  have  zeros  in  the  right  half  plane.  However, 
the  recovery  procedure  is  effective  as  long  as  G(s)  has  no  right  half  plane  zeros  with 
magnitudes  in  frequency  ranges  where  high  loop  gain  is  required. 

Thus,  by  choosing 

L  =  qBW  (3.23) 

Eq.  3.20  will  be  satisfied  and  the  LTR  of  Eq.  3.21  will  be  achieved  as  q  -» <».  Kazerooni 
[  16]  shows  that,  for  a  minimum  phase  plant,  as  q  -4  » the  consequences  of  choosing  L  as 
in  Eq.  3.23  are: 

1 .  The  j  <  n  -  m  finite  closed  loop  eigenvalues,  Xp  of  (  A  -  LC  )  and  (  A  -  BK  -  LC  ) 
approach  the  finite  transmission  zeroes,  zi;  of  the  plant. 

2.  The  remaining  n  -  j  closed  loop  eigenvalues  approach  infinity  in  the  left  half  plane  with 
the  direction  of  their  associated  eigenvectors  being  arbitrary. 

3 .  The  left  closed  loop  eigenvectors,  v  j  of  (  A  -  LC  )  associated  with  the  finite  closed  loop 
estimator  eigenvalues,  X[,  approaches  the  left  zero  direction,  |ij,  of  the  transmission 
zero,  Zj,  which  satisfies 

Zjl  -  A  B 

[  c,T  1 

-C  0 

4.  As  the  infinite  closed  loop  poles  move  away  from  the  imaginary  axis,  an  equal  number 
of  closed  loop  finite  eigenvalues  of  the  compensator  will  cancel  out  the  left  finite 
transmission  zeros  of  the  full  state  model  and  hence  achieving  pointwise  loop  transfer 
recovery 
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K(s)  G(s)  ->  K  (  si  -  A  )*1  B(s) 


Therefore,  for  a  minimum  phase  system,  an  estimator  which  combines  any  stable  infinite 
eigenstructure  with  the  finite  eigenstructure  of  conditions  1  through  3  would  result  in  LTR. 
Using  the  procedure  for  an  unstable  system  results  in  a  loss  of  accuracy  in  the 
approximation  of  the  desired  loop  shape.  If  there  are  no  right  half  plane  zeros  with 
magnitudes  in  frequency  ranges  where  high  loop  gain  is  required,  then  placing  the 
estimator  poles  at  the  mirror  image  of  the  right  half  plane  zeros  achieves  reasonable  LTR. 

IV.  Asymptotic  Linear  Quadratic  Regulator  Theory 

In  LQR  methods,  the  designer  selects  the  weighting  matrix  which  includes  all  the 
performance  requirements  and  is  constrained  to  be  quadratic  in  states  and  controls.  This 
section  describes  a  scheme  for  the  selection  of  quadratic  weights  asymptotically  to  achieve  a 
desired  closed  loop  eigenstructure  as  the  weights  on  the  control  tend  to  zero.  Eigenvalue 
and  eigenvector  properties  are  developed  both  for  the  poles  which  remain  asymptotically 
finite  and  for  those  which  go  asymptotically  to  infinity.  These  properties  are  used  to 
provide  a  complete  unique  characterization  of  quadratic  weights  during  eigenstructure 
assignment  design. 

The  detailed  derivation  of  the  procedure  and  the  proof  of  the  related  theorems  and 
lemmas  are  provided  by  Harvey  and  Stein  [8].  This  section  also  includes  a  summary  of 
asymptotic  properties  of  multivariable  regulators.  A  constructive  procedure  to  define 
quadratic  weights  is  outlined.  This  procedure  is  based  on  modal  design  criteria  and  the 
asymptotic  properties  of  optimal  controllers. 

Consider  the  standard  multivariable  linear  time  invariant  feedback  system 
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x  =  A  x  +  B  u 


(3.24) 


u  =  -  Kx 

(3.25) 

r  =  Hx 

(3.26) 

where  x  is  the  n-dimensional  state  vector,  u  is  an  m-vector  of  controls,  and  H  is  an  m- 
vector  of  responses.  A,  B,  K  and  H  are  constant  matrices  of  appropriate  dimension.  It  is 
assumed  that  the  pair  (  A,B  )  is  controllable,  the  pair  (  A,H  )  is  observable,  and  that  the 
mxm  matrix  HB  has  full  rank  ( i.e.  rank  [  HB  ]  =  m  ). 

The  feedback  gains,  K  are  a  result  of  the  linear  quadratic  (  LQ  )  optimization  problem 
[19].  An  LQ  optimization  problem  for  the  system  in  Eq.  3.24  is  one  described  by  linear 
differential  equations  with  a  performance  index 

oo 

J  =  J  (  xtQx  +  puTRu  )  dt  (3.27) 

o 

represented  by  the  integral  of  a  quadratic  function  of  states  and  control  inputs  where 
Q  =  QT  >  0,  R  =  RT  >  0,  (  A,Q1/2  )  observable  and  scalar  p  >  0.  This  criterion  is 
minimized  by  the  control  function 

u  =  K*  x  (3.28) 

where 

K*  =  R'1  B  -  (3.29) 

P 

P  =  PT  >  0  is  the  unique  positive  definite  solution  of  the  algebraic  Riccati  equation 
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(3.30) 


0  =  P  A  +  AT  P  +  Q  -  P  B  R1  BT  - 

P 

The  behavior  of  both  eigenvalues  and  eigenvectors  are  used  to  provide  complete,  unique 
specification  of  the  weights  Q  and  R,  and  in  particular,  as  the  scalar  p  tends  to  zero. 

The  arbitrary  performance  index  given  in  Eq.  3.27  is  equivalent  to 

oo 

J'  =  J  (  rTr  +  puTRu  )  dt  (3.31) 

o 

with  r  defined  by  Eq.  3.26  form  some  mxn  matrix  H.  These  two  performance  indices 
produce  the  same  gain  matrix,  K*,  as  the  solution  of  their  respective  LQ  optimization 
problems. 

Consider  Q  matrices  of  the  form  Q  =  HTH  with  row  dimension  of  H  equal  to  the 
dimension  of  u.  If  m  is  unity,  then  a  dyadic  Q  suffices.  If  H  is  full  rank,  it  is  always 
possible  to  rearrange  states  and  normalize  H  as  follow 

H  =  WHo  (3.32) 


where 


Ho  =  [Hu  Im  ]  and  W  is  a  non-singular  mxm  matrix. 

Consider  the  linear  quadratic  regulator  from  Eq.  3.28  to  Eq.  3.30  with  performance 
index  given  in  Eq.  3.31  and  normalized  H  as  in  Eq.  3.32.  Assume  further  that 

(1)  rank  [  HB  ]  =  m 
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(2)  the  zeros  of  det  (  H  (  si  -  A  )'!  B  )  lie  on  the  left  half  plane,  are  distinct,  and  do 
not  belong  to  the  c(  A ),  the  spectrum  of  A. 

Then  the  optimally  controlled  system  has  the  following  properties. 

A .  Asymptotically  Finite  Modes 

As  p  — >  0,  there  are  (  n  -  m )  eigenvalues  of  the  form 
s- (p)  -4  s° 

I  s°  I  <  °°,  /=1,  2,  ...,  n  -  m  (3.33) 

with  associated  eigenvectors 

x;  (  p  )  -4  (  s°  -  A  )_1  B  v°  (3.34) 

where  s°  and  v?  are  defined  by 

H  (  s°  -  A  )■>  B  v°  =  0  (3.35) 

The  finite  modes  are  characterized  by  (  n  -  m  )  eigenvalues  s°  plus  (  m  -  1  )x(  n  -  m  ) 
parameters  for  their  associated  directions  v?.  These  provide  a  unique  definition  of  Tl0  via 
Eq.  3.35. 

B .  Asymptotically  Infinite  Modes 


As  p  — >  0,  there  are  m  eigenvalues  of  the  form 


with  associated  eigenvectors 


Xy  (  p  )  Bv“  (3.37) 

where  s “  and  v°J  are  defined  by 

R  =  N-t  S'2  N-1  (3.38) 

WT  W  =  (  H0B  )’t  N-t  N-1  (  HqB  )-’  (3.39) 

and 

N  =  [  v"  Vj  v“  ...  v~  ]  (3.40) 

S  =  diag  (  s~  s“  s”  ...  s~ )  (3.41) 


The  asymptotically  infinite  modes  are  characterized  by  (  m  -  1  )  parameters  for  the 
normalized  eigenvalues  s “  and  (  m  -  1  )x(  m  )  parameters  for  their  associated  eigenvectors 
VJ.  These  provide  unique  definitions  of  R  and  W  via  Eq.  3.38  to  Eq.  3.41. 

Note  that  the  scalar  p  is  the  tradeoff  parameter  for  control  energy  against  the  degree  to 
which  the  asymptotic  properties  are  achieved.  In  addition,  the  eigenvectors  of 
asymptotically  infinite  modes  are  in  the  range  space  of  B.  Hence,  if  the  state  model 
includes  actuator  dynamics  with  B  non-zero  only  for  actuator  rows,  the  asymptotically 
infinite  modes  are  actuator  modes.  The  procedures  for  the  construction  of  Hq  are  given 

below 


1 .  Assume  the  matrix  B  has  been  transformed  to  the  form 
T  0  1  n  -  m 

B21  J  m 


B  = 
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with  mxm  matrix  B2i  non-singular. 


2.  Assume  that  for  each  desired  asymptotically  finite  mode,  s°,  the  eigenvectors  x*  which 
describe  the  desired  modes  can  be  defined  as 

x  ( t )  =  x*  e(  si i  =  1,  2,  n  -  m  (3.42) 


Note  that  only  a  few  of  the  components  in  x*  are  actually  specified,  the  rest  can  be 
arbitrary. 


3 .  Once  the  desired  eigenvectors  are  specified,  reorder  and  partition  x*  as  follows  : 


(3.43) 


where  y*  is  the  specified  subvector,  \)'s  are  the  unspecified  components,  and  {  •••  )R> 
denotes  a  row  reordering  operation. 


4.  The  desired  mode  in  Eq.  3.43  may  not  belong  to  the  set  of  eigenvectors  achievable  by 
LQ  design  : 


V? 

u 


=  {  (s°I-A)-!B  }Ri  v° 


-  L 
.  M  J 


\P: 


(3.44) 


5 .  Using  the  least  square  fit  method  to  select  v?  so  as  to  best  approximate  y*  with  y°; 


v*  =  ((L)TLr'(L)TY* 


(3.45) 
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6.  Assume  that  {  v?  ;  i  =  1,2,  ...,  n  -  m  }  have  been  determined  via  Eq.  3.34.  Then  the 


rollowing  matrix  of  eigenvectors 

X  =  [  x®  Xj  x®  ...  xn°m  ]  (3.46) 

can  be  obtained  where  no  linear  combination  of  the  x°'s  fall  into  the  range  space  of  B. 

7 .  Define  the  projection  matrix 

P  =  In  -  X  [  (  X  )T  X  ]'*  (  X  )T  (3.47) 

as  follows  : 


n-m  n 


and  note  that 

PBp  *  0  if  p  *  0  (3.48) 

8.  The  non-singularity  of  B2i  and  the  condition  given  by  Eq.  3.48  imply  that  P22  is  non¬ 
singular.  Thus,  Hq  is  defined  as 

Ho  =  [  0  P-j'j  ]  P  =  [  P2'2  PT2  I„  J  (3.49) 

Once  the  normalized  matrix  H0  is  obtained  from  Eq.  3.49,  we  can  then  select  the 
eigenvectors  and  eigenvalue  ratios  for  the  asymptotically  infinite  modes.  These  were 
selected  to  achieve  asymptotic  decoupling  of  the  modes.  Matrices  W0  and  R0  were  then 
computed  via  Eq.  3.38  and  Eq.  3.39.  The  final  step  of  the  performance  index  selection 
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procedure  is  to  specify  p.  This  can  be  accomplished  by  computing  LQ  solution  for  a  few 
monotonically  decreasing  values. 
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Chapter  4 


Mathematical  Model 


I.  Introduction 

The  dynamic  stability  and  response  to  control  inputs  of  the  rotorcraft  rigid  body 
dynamics  define  its  handling  qualities.  The  six  rigid  body  degrees  of  freedom  are  the  basic 
motions  involved  in  the  analysis  of  flight  dynamics.  The  rotor  is  also  a  major  factor  in  the 
dynamic  response  characteristics  of  helicopters,  and  therefore  it  is  necessary  to  consider  the 
rotor  degrees  of  freedom  as  well,  particularly  the  tip  path  plane  flapping  motion.  The 
arrangement  of  the  rotor  or  rotors  on  a  helicopter  is  also  an  important  factor  in  its  behavior, 
notably  its  stability  and  control  characteristics.  Because  of  the  rotation  of  the  rotor,  there  is 
considerable  coupling  between  the  longitudinal  and  lateral  directional  motions. 

Since  the  handling  qualities  of  the  helicopter  have  different  characteristics  in  hover  and 
in  forward  fight,  only  the  hover  regime  will  be  analyzed  here.  The  hovering  analysis  is 
simpler  than  forward  flight  because  of  the  axisymmetry  of  the  rotor  aerodynamics  in  hover. 
The  equations  are  developed  in  this  Chapter  in  the  following  way.  The  helicopter  is  first 
regarded  as  a  single  rigid  body,  and  the  equations  of  motion  are  derived  with  respect  to  a 
set  of  axes  fixed  in  it.  The  derivation  will  be  developed  for  the  single  main  rotor  and  tail 
rotor  configuration.  This  configuration  of  helicopter  uses  the  tail  rotor  to  provide  the  anti¬ 
torque  in  order  to  balance  the  torque  produced  by  the  main  rotor  and  the  yaw  control. 

II.  Equations  of  Motion 
A.  Rigid  Body  Dynamics 
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Consider  the  rigid  body  motions  of  a  helicopter  with  a  single  main  rotor  and  a  tail 
rotor.  It  is  assumed  that  the  helicopter  has  complete  axisymmetry,  although  the 
longitudinal  -  lateral  dynamics  are  coupled,  the  vertical  and  the  longitudinal-lateral 
dynamics  are  completely  separated.  Such  separation  is  a  basic  feature  of  the  rotor  in  hover, 
and  generally  holds  for  the  hover  handling  qualities  even  though  the  entire  helicopter  is  not 
exactly  axisymmetric. 


T,  -0 


t  v,  ny 

l  X 


Ttr 


Fig.  4. 1  Definition  of  the  Rigid  Body  Degrees  of  Freedom  of  a  Helicopter,  and  the  Forces 
and  Moments  due  to  the  Main  Rotor,  Tail  Rotor  and  Gravity. 


A  body  fixed  orthogonal  coordinate  system  is  used,  with  the  origin  at  the  helicopter 


center  of  gravity  (  mass  ).  The  z-axis  is  vertical  with  the  center  of  gravity  directly  under  the 
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rotor  hub.  The  x-axis  is  positive  forward,  the  y-axis  is  positive  to  the  right,  and  the  z-axis 
is  positive  downward  and  is  lying  in  the  plane  of  symmetry.  It  is  assumed  that  these  axes 
are  the  principal  axes  of  the  helicopter,  so  inertial  cross  coupling  between  the  yaw  and  roll 
axes  is  neglected.  Figure  4.1  shows  the  free  body  diagram  of  the  helicopter. 

The  linear  displacement  of  Lite  center  of  gravity  has  the  components  xB,  yB,  and  zB;  the 
roll,  pitch,  and  yaw  angular  rotations  of  the  helicopter  about  the  center  of  gravity  are  given 
by  the  Euler  angles  0B,  0B,  and  \|/B.  The  helicopter  angular  velocity  has  components  in  the 
body  axes 

P  =  <t>B  '  Vb  sin  (0B)  (a) 

q  =  0Bcos(0B)  +  cos  (0B)  sin  (<>B)  (b)  (4.1) 

r  =  -  0B  sin  (0B)  +  \j/B  cos  (0B)  cos  (<j>B)  (c) 


inverting  these  equations  gives  the  angular  rates 

4>b  =  P  +  <7  ^n  (0B)  sin  (<J>B)  +  Man  (0B)  cos  (<}>B)  (d) 

0B  =  q  cos  (<|)B)  -  r  sin(0B)  (e)  (4.2) 


r( 


cos (0b) 
cos (0B) 


(0 


The  components  of  the  linear  velocity  of  the  center  of  gravity  V,  are 

u  =  xB  cos  (0B)  cos  (\(/B)  +  yB  [  sin  (0B)  sin  (0B)  cos  (vj/B)  -  cos  (0B)  sin  (\j/B)  )  + 
zB  [  cos  «>B)  sin  (0B)  cos  (yB)  +  sin  (0B)  sin  (V|/B)  ]  (a) 

v  =  xH  cos  (0B)  sin  (y)B)  +  yB  [  sin  (<>g)  sin  (0B)  sin  (\|/B)  +  cos  (0B)  cos  (\|/B)  ]  + 
zB  |  cos  (0B)  sin  (0B)  sin  (\|/B)  -  sin  (0B)  cos  (\jtB)  ]  (b)  (4.3) 
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W  =  -  Xq  sin  (0B)  +  yB  sin  (4>g)  cos  (0B)  +  zB  cos  (<J)B)  cos  (0B)  (c) 


q  y 

The  translational  acceleration  components  of  along  the  body-fixed  axes  will 
contain  not  only  the  components  C^p  and  but  also  centripetal  acceleration 

components  due  to  the  rotation  of  the  body-fixed  axis  *  sto.vi  relative  to  inertial  space.  The 
set  of  linear  accelerations  is  given  by 
d  u 

-jj-  =  u  +  qw  -  rv  tai 


dv 

=  v  +  ru-  pw 


(b)  (4.4) 


w  +  pv  -  qu 


By  Newton’s  2nd  Law  for  a  rigid  body,  force  is  equal  to  the  rate  of  change  of  linear 
momentum.  The  mass  of  the  helicopter,  m,  remains  essentially  constant  during  a 
maneuver.  The  force  F  is  resolved  into  components  Fx,  Fv,  and  F,  along  the  x,  y,  and  : 
body  axes,  and  the  three  scalar  equations  can  be  written  as 

Fx  =  m  =  m(  u  +  qw  -  rv  )  (a) 


Fy  =  m  =  m(  v  +  ru  -  pw  ) 


(b)  (4.5) 


r-  dw  .  •  . 

Fz  =  m  -jp  =  m  (w+pv-qu) 


The  forces  and  moments  that  are  generated  by  the  rotor  systems  and  all  other 
aerodynamically  generated  forces  and  moments  are  considered  to  be  the  external  forces/ 
moments  in  the  equations  of  motion.  All  forces  and  moments  are  expressed  relative  to  the 
body  center  of  gravity  for  use  in  the  six  degrees  of  freedom  rigid  body  equations  of 
motion.  The  external  aerodynamic  and  propulsive  forces  and  moments  are  denoted  by  X. 
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Y,  Z,  and  L,  M,  N  respectively.  Since  the  center  of  gravity  coincides  with  the  center  of 
mass,  the  gravity  force  produces  no  external  moments  about  that  point.  Gravity  can 
contribute  components  only  to  the  summation  of  the  external  forces.  The  components  of 
the  weight  in  the  directions  of  the  axes  are 


Xg  =  -  mg  sin  (0B) 

Yg  =  mg  cos  (0B)  sin  (<J>B) 
Zg  =  mg  cos  (0B)  cos  (<{>B) 


(a) 

(b)  (4.6) 

(c) 


with  the  aerodynamic  forces  ( including  propulsive  forces  ),  the  resultant  external  forces  are 


Fx  =  X  -  mg  sin  (0B) 


Fv  =  Y  +  mg  cos  (0B)  sin  (<J)B) 


Fz  =  Z  +  mg  cos  (0B)  cos  (<(>3) 


(a) 

(b)  (4.7) 

(c) 


The  resulting  equations  of  motion  for  the 
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translational  degrees  of  freedom  are 


f  -sin0B  ^ 

(X  } 

cos0B  sin<J)B 

1 

+  — 
m 

Y 

(4.8) 

^  COS0B  COS(J)B  J 

U ) 

Since  gravity  does  not  contribute  to  the  moments  about  the  mass  center,  then  the 
external  couples,  L,  M,  and  N,  acting  on  the  helicopter  are  entirely  aerodynamic.  For  all 
practical  purposes,  the  inertia  tensor  is  time  invariant.  The  body  fixed  x:  plane  of  the 
helicopter  is  further  assumed  to  be  the  plane  of  symmetry,  so  the  products  of  inertia  1^  and 
Iv,  are  zero.  Given  these  assumptions,  the  equations  for  the  external  moments  are 
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Equations  4.8  and  4.10  are  the  six  degree  of  freedom,  rigid  body  equations  of  motion. 
If  equations  4.8  and  4.10  are  combined  with  the  kinematical  equations  governing  the 
evolution  of  the  Euler  angles  (  Eqs.  4.1  or  4.2  ),  then  the  dynamics  of  the  helicopter  can  be 
described  by 
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X  =  f(x,u) 


(4.11) 


where  x  is  the  state  vector  and  u  is  the  input  vector.  Equations  4.1,  4.8,  and  4.10  contain 
products  of  the  dependent  variables,  some  of  which  appear  as  transcendental  function  and 
are  in  general  nonlinear. 

The  equation  of  motion  are  frequently  linearized  for  the  use  in  stability  and  control 
analyses.  It  is  assumed  that  the  motion  of  the  helicopter  consists  of  small  perturbations 
from  a  trim  condition  of  steady  flight.  Equation  4.1 1  can  then  be  linearized  about  a  trim 
flight  condition  to  yield  the  coefficients  of  a  linear,  first  order  set  of  differential  equations 
that  represents  the  rigid  body  dynamics  of  the  helicopter  for  small  perturbations.  The 
operating  point  equations  are  obtained  by  recognizing  that  zero  translational  and  rotational 
accelerations  are  implicit  in  the  concept  of  a  trim  condition.  Then  denoting  the  trim 
condition  with  zero  subscript,  the  perturbed  angular  velocity  has  components 

p  =  <>b  -  vb  sin  <e0) 

Q  =  0B  (4.12) 

r  =  yB  COS  (60) 

The  components  of  the  linear  velocity  perturbation  are 

u  =  (  jc0  +  )  cos  (0O)  -  Jt0  0g  sin  (0n)  +  zB  sin  (0O) 

V  =  ,t0\j/Bcos  (0O)  +  >'B  (4.13) 

w  =  -  ( ,x0  +  *B  )  sin  (0O)  -  Xq  0b  cos  (0O)  +  zB  cos  (0O) 

The  differential  equations  are  in  the  standard  state  space  representation 
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5x  =  A  8x  +  B  8u 


where  8x  represents  the  perturbation  from  the  trim  condition  of  the  state  variables  and  8u 
represents  the  deviations  of  the  control  displacements  from  trim.  Although  the  components 
of  A  and  B  matrices  are  function  of  forward  velocity,  they  are  often  approximated  by  an 
average  value  so  that  the  resulting  linear  model  is  velocity  and  time  invariant.  The 
linearized  six-degree-of  freedom  representation  of  the  helicopter  dynamics  is  valid  only  if 
the  initial  body  angular  rates  are  zero. 

15 .  Rotor  Dynamics 

The  rotor  of  a  helicopter  is  used  for  generating  lift,  forward  thrust  and  control 
moments.  Consequently,  modeling  the  dynamics  of  the  rotor  is  very  important  in 
accurately  representing  the  overall  dynamic  response  of  a  helicopter.  The  rotor  model  used 
in  this  study  is  the  tip  path  plane  representation  of  blade  flapping  [3].  This  model  is 
commonly  used  in  helicopter  dynamics  and  control  studies. 

The  blade  loads  and  motions  of  a  helicopter  blade  are  periodic  in  the  azimuth  angle  \\f, 
and  the  blade  flapping  angle  p  can  be  represented  as  a  Fourier  series  in  \\ / 

oo 

P(  t ,y )  =  Pn(  t )  +  I[  pnc(  t )  cos  (ny)  +  pns(  t )  sin  (ny)  ]  (4.14) 

n=l 

The  azimuth  angle  is 

y  =  fit  (4.15) 

where  Q  is  the  rotor  angular  speed.  Since  the  frequency  of  the  rotor,  approximately 

30  rad/s,  is  so  much  greater  than  that  of  the  rigid  body  motions,  less  than  4  rad/s,  only  the 

lowest  harmonics  in  the  Fourier  expansion  for  the  blade  flapping  angle  are  required  to 
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describe  the  motions  of  the  rotor  adequately  for  purposes  of  dynamic  stability  analysis  and 
control  law  synthesis.  In  fact  the  tip  path  plane  model  uses  only  the  first  harmonics  of 
blade  motion. 

P(  t,y  )  =  p0(  t )  +  plc(  t )  cos  (y)  +  pls(  t )  sin  (y)  (4.16) 


rotor 


viewed  from  the  left  viewed  from  aft 

Fig.  4.2  Blade  Flapping  Harmonics. 


A  graphical  interpretation  of  the  blade  cone  angle  P0(  t  )  and  first  two  flapping 
harmonics  pls(  t )  and  plc(  t )  are  shown  in  Fig.  4.2.  The  two  flapping  harmonics  define  a 
plane  known  as  the  tip  path  plane.  The  higher  harmonics  distort  this  plane  but  are  not 
needed  if  only  rigid  body  motions  are  to  be  modeled. 

A  diagram  of  the  forces  acting  on  the  rotor  is  shown  in  Fig.  4.3.  The  lift  on  the  blade 
is  denoted  as  L  and  the  components  of  the  moment  created  by  the  lift  about  the  c.g.  of  the 
helicopter  are 

McGx  =  *  h  L  sin  (Pj)  sin  (\|i)  =  -  h  L  p,  sin  (\j/)  (4.17) 

McGy  =  -  h  L  sin  (pj)  cos  (\}/)  ^  -  h  L  Pj  cos  (v|/)  (4. 1 8) 
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The  angle  p 
P,(t)  ‘ 


Plane  view,  view  from  forward 


Pig.  4.3  Rotor  Force  Diagram. 

is  represented  by  the  first  two  elements  in  its  Fourier  series  expansion, 
Pic(  t )  cos  (vj/)  +  pls(  t )  sin  (\j/)  (4.19) 
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Since  \\i  =  Q  t,  these  components  are  periodic  functions  of  time.  If  the  moment 
expressions  given  in  Eqs.  4.17  and  4.18  were  inserted  directly  into  the  rigid  body 
equations,  a  set  of  ordinary  differential  equations  with  time  variable  coefficients  would 
result.  Few  theories  exist  for  the  design  of  control  laws  for  systems  described  even  by 
linear  time  variable  ordinary  differential  equations.  Furthermore,  those  theories  which  are 
available  lead  to  time  variable  gains  which  are  not  desirable  because  of  difficulty  of 
implementation.  Fortunately,  the  frequency  of  the  rotor  is  so  much  greater  than  that  of  the 
rigid  body  motions  that  the  averaged  values  of  the  rotor  moments  can  be  used  in  the  rigid 
body  equations.  These  averaged  equations  are  obtained  by  integrating  Equations  4.17  and 
4. 1 8  with  respect  to  \\r  over  the  interval  2tt,  and  result  in 

Mcoxavo  -  'TPi>  <4-20) 

Mcgyavc  =  (4:1) 


These  expressions  can  be  substituted  into  the  rigid  body  equations  to  give  a  set  of  ordinary' 
differential  equations  which  give  a  good  representation  of  the  rigid  body  motions  of  the 
helicopter  fuselage. 


Since  the  blades  themselves  have  inertia,  their  motions  must  be  described  by 
differential  equations.  The  dynamics  of  the  coning  angle  have  little  effect  on  rigid  body 
dynamics  and  can  be  ignored;  however,  accurate  modeling  of  the  rigid  body  motions 
requires  that  the  dynamics  of  the  harmonic  coefficients  be  included.  The  harmonics  are 
described  by  a  set  of  coupled  second  order  linear  differential  equations 


■  1  0  if  P 

.0  ij[p 


1  S 

1c  - 


1 


+  c 


E 


-E 

1 


“  ■“ 1 

Pis 

+  k 

-  Pic- 

- 

lc  J 


-76- 


u 


-f  a  -g 
a  f  v 


v 

g 


_ 

V 

+  b 

i  -y 

p 

.  Y  1  . 

L  q  J 

e 


lateral  cyclic 


e 


longitudinal  cyclic  J 


(4.22) 


The  "damping"  and  "stiffness"  matrices  are  determined  by  aerodynamic  damping, 
centrifugal  and  structural  stiffness  of  the  blade,  and  gyroscopic  forces.  The  inputs  are  the 
latera'  and  longitudinal  cyclic  pitch.  The  above  flapping  equations  are  coupled  with  the 
rigid  body  equations  result  in  a  12th  order  set  of  differential  equations  (  an  8th  order  rigid 
body  model  and  two  coupled  2nd  order  models  representing  the  rotor  dynamics  ).  There  is 
also  coupling  with  the  rigid  body  linear  and  angular  velocities.  More  complex  rotor  models 
which  include  lag  and  torsional  dynamics  can  be  formulated  but  12th  order  models  of  the 
type  described  above  appear  to  be  adequate  for  control  law  synthesis  for  most  articulated 
rotor  helicopters. 

III.  Design  Model 


The  mathematical  model  used  for  design  represents  a  modern  attack  helicopter  in  the 
hover  condition.  The  state  space  representation  is  a  linear,  coupled,  eight  degree  of 
freedom  helicopter  model  which  includes  eight  rigid  body  states  (  u,v,w,p,q,r,0,0  ),  and 
two  second  order  modes  of  the  advancing  and  regressing  rotor  tip  path  plane  (  a^b]  ). 
Main  rotor  collective  pitch,  lateral  cyclic  pitch,  longitudinal  cyclic  pitch,  and  tail  rotor 
collective  pitch  are  the  control  inputs.  These  equations  are  written  in  the  standard  form 

x12  =  ^12  x12  +  ®12  u  (4.23) 

where  A!2,  B12,  are  dimensional  matrices  containing  the  analytical  data  for  the  linearized 
dynamics  of  a  helicopter  and  are  given  in  Appendix  A. 


-77- 


The  sensors  measure  vertical  velocity,  and  roll,  pitch  and  yaw  rates,  and  the  output 


equation  can  be  written  as 


y  =  Cx12 


(4.24) 
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Fig.  4.4  Singular  Value  Plot  of  12th  Order  Open  Loop  Helicopter  at  Hover. 


The  singular  values  for  the  state  space  equations  4.23  and  4.24  that  represent  the  12th  order 
helicopter  at  hover  are  plotted  in  Fig.  4.4.  The  twelve  eigenvalues  and  their  respective 
modal  identifications  are  given  in  Table  4.1. 


As  can  be  seen  from  Trble  4.1,  the  rotor  dynamics  result  in  a  set  of  high  and  a  set  of 
relatively  low  frequency  eigenvalues.  The  high  frequency  set  of  rotor  eigenvalues  have  a 
natural  frequency  of  73.4  rad/s  and  a  damping  factor  of  0.189.  This  is  sometimes  called 
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the  rotor  advancing  flapping  mode.  Examination  of  the  residues  of  this  mode  indicate  that 
it  does  not  have  much  affect  on  the  rigid  body  modes.  On  the  other  hand,  the  lower 
frequency  set  of  rotor  eigenvalues  has  a  natural  frequency  of  12.3  rad/s  -nd  a  damping 
factor  of  0.96.  The  mode  associated  with  this  set  of  eigenvalues  is  known  as  the 
regressing  flapping  mode.  Since  the  rigid  body  roll  mode  has  an  eigenvalue  of 
-4.3434  rad/s,  the  regressing  flapping  mode  will  affect  the  rigid  body  dynamics. 
Examination  of  the  residues  of  this  eigenvalues  verifies  that  it  does  have  a  significant 
impact  on  pitch  and  roll  rate  responses. 

Twelfth  Order  Model 


-4.3434  + 

.OOOOi 

roll  rate 

.2130  ± 

,5270i 

forward  velocity 

.0328  ± 

,7499i 

side  slip 

-.9421  + 

.OOOOi 

pitch  rate 

-.5654  + 

.OOOOi 

yaw  rate 

-.3234  + 

.OOOOi 

heave  velocity 

-13.6446  ± 

72. 1554i 

advancing  rotor  mode 

-11.7698  ± 

3.6620i 

regressing  rotor  mode 

Table  4.1  Open  Loop  Eigenvalues  and  Corresponding  Modes  of  12th  Order  Model. 

Thus  the  design  model  should  contain  the  effects  of  rotor  dynamics.  In  practice,  the 
measurement  of  the  rotor  states  appears  impossible  with  current  technology.  It  was 
decided  to  develop  a  reduced  eighth  order  design  model  which  explicitly  contained  only  the 
rigid  body  modes.  The  use  of  a  lower  order  model  for  control  law  design  simplifies  the 
design  and  implementation  of  the  feedback  compensator.  Two  methods  of  order  reduction; 
truncation  and  residualization,  were  applied  to  the  twelfth  order  model. 

Truncation  involved  eliminating  or  truncating  the  fast  modes  from  the  original  model. 
In  this  case,  the  tip  path  plane  modes  (  rotor  flapping  modes  )  w'ere  truncated  from  the 
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twelfth  order  model.  In  Eq.  4.25,  x8  represents  the  eight  rigid  body  states  and  x4 
represents  the  4  rotor  states.  Truncation  involves  crossing  out  the  rows  and  the  columns  of 
the  A  matrix  correspond  to  x4,  and  the  rows  of  the  B  matrix  correspond  to  x4,  leaving  just 
the  eight  rigid  body  equations  as  the  nominal  model. 


oo 

< 

oo 

00* 

< 

1 _ 

'  x8  ' 

■  Bg  - 

= 

+ 

< 

00 

< 
_ 1 

-  x4  - 

.  b4  . 

(4.25) 


For  residualization,  the  higher  frequency  motions  are  assumed  to  reach  steady  state 
conditions  so  rapidly  that  the  time  rate  of  change  of  these  states  (  x4  )  may  be  set  to  zero. 
Following  the  procedure  shown  in  Eq.  4.26,  this  assumption  allowed  unique  algebraic 
solution  for  the  four  tip  path  plane  states  in  terms  of  the  eight  rigid  body  states.  These  four 
state  equations  were  then  substituted  into  the  rigid  body  equations  and  placed  in  standard 
state  space  form  in  order  to  obtain  the  residualized  nominal  model. 


x8 

OO 

< 

00 

00 

< 

1 _ 

■  x8  " 

•  b8- 

= 

+ 

.  0  . 

1 

Tf 

< 

°0 

< 

_ 1 

-  x4  - 

-  B4_ 

x4  -  -  A4  4  1  A4  8  x8  -  A4  4  1  B4  u 
=*  xg  =  (  Ag  g  -  Ag  4  A4  4  1  A4  8  )  Xg  +  (  Bg  -  A8>4  A4  4  1  B4  )  u 


The  nominal  models  resulting  from  truncation  and  residualization  were  evaluated  by 
comparing  their  open  loop  eigenvalues  with  the  eight  rigid  body  eigenvalues  obtained  from 
the  original  twelfth  order  model.  As  can  be  seen  from  Table  4.2,  the  eigenvalues  of  the 
residualized  model  matched  the  eigenvalues  of  the  twelfth  order  model  much  more  closely 
than  did  those  of  the  truncated  model.  Thus  the  residualized  model  was  used  as  the  design 
model. 
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-4.3434  +  .OOOOi  -1.4032  +  .OOOOi 

.2130  ±  .527 Oi  .2705  ±  .83381 

.0328  ±  .7499i  .2076  ±  ,4825i 

-.9421  +  .OOOOi  -.6249  +  .OOOOi 

-.5654  +  .OOOOi  -.4476  +  .OOOOi 

-.3234  +  .OOOOi  -.1939  +  .OOOOi 

-13.6446  ±  72. 1554i 
-1 1.7698  ±  3.6620i 


-3.2377  +  .OOOOi  roll  rate 
.2110  ±  .5296i  forward  velocity 
.0353  ±  .743 li  sideslip 
-.9021  +  .OOOOi  pitch  rate 
-.5695  +  .OOOOi  yaw  rate 
-.3221  +  .OOOOi  heave  velocity 

advancing  rotor  mode 
regressing  rotor  mode 


Table  4.2  Comparison  of  Open  Loop  Eigenvalues  of  12th  Order  Model  with  8th  Order 

Truncated  and  Residualized  Models. 


The  design  model  is  written  in  standard  state  space  form  as  given  in  Eq.  4.27  where 
A8  and  Bg  are  obtained  from  Eq.  4.26.  The  measured  outputs  are  heave  velocity,  roll, 
pitch  and  yaw  rates.  The  dimensional  A8,  Bg  and  Cg  matrices  are  given  in  Appendix  A. 

x8  =  AgXg  +  B8u  (4.27) 

y  =  C8  x8  (4.28) 


The  singular  value  plot  of  the  8th  order  residualized  model  given  by  Eqs.  4.27  and 

4.28  is  shown  in  Fig.  4.5.  As  compared  to  the  12th  order  model  shown  in  Fig.  4.4,  the 
singular  values  for — ,  — ^ ^ - ,  and are  essentially  the  same 

^collective  ^lateral  ^longitudianl  ^directional 

for  both  models  for  frequencies  less  than  10  rad/s.  This  similarity  is  also  verified  by  the 
comparison  of  the  open  loop  eigenvalues  in  Table  4.2. 
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FREQUENCV  <rad/s) 

-  w/SeollectiYe 

. f/S  directional 

- p/5  lateral 

■  — q/S  longitudinal 

Fig.  4.5  Singular  Value  Plot  of  8th  Order  Residualized  Helicopter  Model. 

Due  to  the  fact  that  the  elements  of  a  given  dimensional  eigenvector  have  different 

units,  (linear  and  angular  units),  it  is  difficult  to  identify  the  dominant  modes  and  couplings 

associated  with  a  particular  dimensional  eigenvector.  Hence  it  is  desirable  to  non- 

dimensionalized  the  state  and  control  vectors  so  that  modal  coupling  can  be  more  easily 

identified.  The  maximum  expected  value  of  each  state  variable  applicable  to  the  model  was 

used  for  non-dimensionalization.  All  linear  velocities  were  divided  by  16.877  ft./s 

(  10  knots  ),  angular  velocities  by  0.349  rad/s  (  20  deg./s  ),  and  angular  displacement  by 

0.349  rad  (  20  degrees  ).  The  control  inputs  were  non-dimensionalized  by  dividing  by 

their  maximum  values.  The  collective  pitch  was  divided  by  0.157  rad  (  9  degrees ),  the 

lateral  cyclic  pitch  by  0.153  rad  (  8.75  degrees ),  the  longitudinal  cyclic  pitch  by  0.262  rad 

(  15  deerees  ).  and  the  tail  rotor  collective  pitch  by  0.323  rad  (18.5  degrees ).  The  open 
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loop  eigenvalues  and  non-dimensional  eigenvectors  for  the  residualized  hover  model  are 
tabulated  in  Table  4.3. 


ro II  rate 


side  slip 


forward  veloci ty 


X  = 

-3.2377  + 

. 0000 i 

X  =  .0353  ± 

. 743 1 i 

X  =  .2110 

±  . 5296 i 

u 

.0071  + 

.  0000 i 

. 1239  ± 

.  0691  i 

.  1078 

T 

. 4243 i 

o 

.0692  + 

.  0000 i 

-. 1740  ± 

.  5073 i 

-.1781 

T 

.  1 191  i 

w 

.0063  + 

,0000i 

-.0429  4 

.  0285 i 

-.0023 

± 

. 0273 i 

p 

.9407  + 

.  0000 i 

.0898  T 

.  4337 i 

.  1  138 

T 

.03011 

q 

-.0308  + 

.  0000 i 

.0936  ± 

.  0437 i 

- . 0952 

? 

. 1698 i 

r 

-.  1557  + 

,0000i 

.3095  ± 

.  1069  i 

-.4358 

T 

. 57 12 i 

$ 

-.2913  + 

.  0000 i 

-.5791  T 

.  14  19  i 

.0433 

? 

. 2 102 i 

e 

.0112  + 

.  0000 i 

.0703  T 

.  1368  i 

-.3800 

± 

. 0565 i 

p i tch  rate 

yaw  rate 

heave  ve 1 oc i ty 

X 

=  -.9021  + 

. 0000 i 

X  =  -.5695 

+  .OOOOi 

X  =  -.3221 

+  .0000 

u 

.3738  + 

. 0000 i 

-.  1589 

+ 

.  0000  i 

-.0858  + 

. 0000 i 

V 

.0364  + 

. 0000 i 

.0710 

+ 

.OOOOi 

-. 1097  + 

.OOOOi 

UJ 

.  1027  + 

. 0000 i 

- . 2290 

+ 

.OOOOi 

-.3863  + 

.OOOOi 

p 

.0031  + 

. 0000 i 

.0633 

+ 

.OOOOi 

-.0317  + 

.  0000  i 

q 

-.4160  + 

. 0000 i 

.0413 

+ 

.OOOOi 

.0419  + 

. 0000 i 

r 

-.6597  + 

. 0000 i 

.9421 

+ 

.  0000  i 

-.9080  + 

.OOOOi 

♦ 

-.0145  + 

. 0000 i 

-.0857 

+ 

.OOOOi 

.0550  + 

.OOOOi 

© 

.4860  + 

. 0000 i 

-. 1291 

+ 

.OOOOi 

-.0332  + 

.OOOOi 

Table  4.3  Open  Loop  Eigenvalues  and  Associated  Non-Dimensional  Eigenvectors. 


The  values  in  Table  4.3  show  that  cross  coupling  exists  between  all  modes.  Roll  rate 
is  coupled  to  yaw  rate.  Side  slip  is  coupled  to  yaw  rate  and  forward  velocity.  Forward 
velocity  is  coupled  to  side  slip,  roll  and  yaw  rate.  Pitch  rate  is  coupled  to  yaw  rate  and 
heave  velocity.  Yaw  rate  and  heave  aie  also  coupled  to  each  other.  It  is  also  evident  that 
there  are  two  pairs  of  unstable  complex  conjugate  roots  for  side  slip  and  forward  velocity 
with  a  time  to  double  amprtudc  of  19.5  seconds  a/.d  3.27  seconds  respectively.  Also  the 
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bandwidths  for  pitch  rate,  yaw  rate  and  heave  modes  are  far  lower  than  desired.  These 
factors  greatly  degrade  the  handling  qualities  of  the  helicopter. 

Besides  the  modal  coupling  described  above,  examination  of  the  control  distribution 
matrix  B  reveals  strong  coupling  of  yaw  angular  acceleration  in  all  inputs.  Collective  cyclic 
pitch  has  primary  effect  on  yaw  acceleration,  and  about  equal  effect  on  heave  velocity  and 
roll  angular  acceleration.  Lateral  cyclic  has  strong  effect  on  roll  acceleration,  and  pitch 
angular  acceleration  is  significantly  influenced  by  the  longitudinal  cyclic  pitch  only.  Lateral 
and  longitudinal  cyclic  pitch  have  about  equal  effect  on  yaw.  Roll  and  yaw  acceleration  are 
coupled  considerably  in  tail  rotor  collective  pitch. 

IV.  Error  Model 

An  error  model  which  gives  an  approximation  of  the  modeling  uncertainties  is 
developed  so  that  the  maximum  attainable  bandwidth  for  the  flight  control  system  that 
ensures  stability  robustness  can  be  estimated  and  the  performance  robustness  determined. 

Hying  quality  problems  of  modem  flight  control  systems  are  most  often  related  to  the 
introduction  of  time  delays  in  the  rotorcraft  flight  control  system.  These  time  delays  can 
produce  a  significant  degradation  in  the  flying  qualities  of  the  rotorcraft  during  demanding 
tasks.  Therefore  the  error  model  used  in  this  study  approximates  the  unmodeled  high 
frequency  dynamics  of  rotor  dynamics,  structural  modes,  actuators,  filters,  sensors, 
sampling  delays,  and  computational  delays  associated  with  digital  implementation  of  the 
flight  control  system  with  effective  time  delays.  The  total  effective  time  delay  for  each  of 
the  four  input  channels  is  assumed  to  be  the  same.  Additional  assumptions  are:  (1)  there  is 
no  cross-coupling  of  the  uncertainties  between  channels,  and  (2)  the  nominal  design  model 
has  a  10%  error  at  low  frequencies. 
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The  multiplicative  (  relative  )  error  model  [9],  given  by  Eq.  4.29  is  utilized  in  the 
calculation  of  the  unmodeled  dynamics  error  bound  for  the  stability  robustness  test.  The 
multiplicative  error  matrix  is  defined  as 

Gtrue(s)  =  G(s)(  I  +  E(s) )  (4.29) 

and  is  related  to  the  effective  time  delay,  At  by  the  following  equation  [10], 

I  +  E(s)  =  diag{  e'sAx  }  (4.30) 


Using  a  first  order  Pade  approximation  for  the  total  effective  time  delay,  the  error  can  then 
be  evaluated  as 


E(s)  =  ( 


-  2Ats 
Ats  +  2 


)I 


(4.31) 


V .  Simulation  Model 


This  section  formulates  a  helicopter  mathematical  model  suitable  for  evaluation  of 
control  laws.  The  basic  simulation  model  will  consist  of  12th  order  dynamics  model,  i.e. 
8th  order  rigid  body  dynamics  plus  two  second  order  rotor  flapping  dynamics,  first  order 
actuators  with  bandwidths  of  50  rad/s,  and  a  first  order  Pade  approximation  of  the 
computational  and  filter  delays  as  shown  in  the  block  diagram  below. 


Fig.  4.6  Block  Diagram  of  the  Simulation  Model. 


A  more  realistic  model  will  include  nonlinear  dynamics  and  kinematics  model, 
computational  and  filter  delays,  and  rate  and  position  limits  in  the  actuator  dynamics  as 
given  in  Fig.  4.7. 


actuator 


actuator  deflection 

rate  limiter  limiter 


Fig.  4.7  Actuator  Dynamics. 
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Chapter  5 


Nominal  Design 


I.  Introduction 

The  purpose  of  this  chapter  is  to  develop  a  simple  control  law  structure  for  low  speed 
and  hovering  control  of  an  advanced  combat  helicopter.  In  the  first  section  of  this  chapter, 
the  philosophy  for  the  feedback  system  design  is  presented  and  in  the  second  section,  a 
methodology  using  eigenstructure  assignment  to  design  the  inner  loop  control  law  is 
discussed.  In  the  third  section,  the  use  of  a  simple  feedback  or  feedforward  outer  loop  that 
wraps  around  the  inner  loops  to  achieve  selectable  command  response  types  is  discussed. 

II.  Design  Philosophy 

A  modem  four-bladed  attack  helicopter  operating  at  low  speed  or  hover  flight  regimes 
is  used  to  illustrate  the  design  of  the  flight  control  system.  The  dynamic  response 
characteristics  for  the  helicopter  modeled  in  this  study  are  typical  of  most  high  performance 
helicopters.  Simulations  and  flight  tests  have  shown  that  with  such  helicopter  dynamics, 
even  experienced  helicopter  pilots  are  unable  to  accomplish  divided  attention  operations  or 
relatively  simple  tasks  in  degraded  visual  environments.  Handling  qualities  that  minimize 
the  involvement  of  the  pilot  in  basic  stabilization  tasks  are  required  to  accomplish  such 
operations.  This  requires  a  high  bandwidth  (  high  gain  ),  multiply-redundant  flight  control 
system. 

As  discussed  by  Tischler  [1],  high  gain  flight  control  systems  require  high  bandwidth 
actuators.  Hence  the  rate  and  deflection  limit  characteristics  of  the  actuators  will  impose 
significant  limitations  on  the  feedback  gains.  Additional  factors  that  limit  the  maximum 
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feedback  gains  are  :  (1)  sensor  noise  amplification,  (2)  in-plane  ( lead-lag  )  rotor  coupling 
and  inflow  dynamics  [2],  (3)  phase-margin  requirements  and  high  frequency  modeling 
uncertainty  (  rotor  and  structural  flexure  modes )  and  (4)  total  effective  time  delays  [3]-[4] 
due  in  part  to  digital  implementation  of  the  flight  control  laws. 

This  chapter  presents  a  control  law  structure  that  will  provide  good  command  tracking 
(  decoupling  ),  gust  attenuation,  stability  robustness,  and  meet  the  proposed  handling 
qualities  specifications  in  the  face  of  errors  in  the  mathematical  model  used  for  design.  In 
this  Chapter,  the  performance  of  the  system  is  evaluated  using  the  design  model.  In 
Chapter  6,  the  performance  of  the  control  law  will  be  demonstrated  by:  (1)  time  response 
simulations  with  more  realistic  models,  (2)  unstructured  singular  value  analysis,  (3) 
structured  singular  values  (  p.  analysis  [5]-[8]  ).  The  control  law  structure  consists  of 
"inner"  and  "outer"  feedback  loops.  This  chapter  also  descnoes  an  eigenstructure 
methodology  used  to  design  inner  loop  control  laws  which  decouple  the  vertical  velocity, 
yaw,  longitudinal,  and  lateral  modes.  The  resulting  closed  inner  loop  system  has  an 

G)r 

integral  or-^p  [9]-[  1 3]  ( <oc  is  the  desired  crossover  frequency  )  characteristic  in  the  lateral 

and  longitudinal  attitude  loops  and  first  order  type  responses  in  the  vertical  velocity  and 
yaw  rate  loops.  If  altitude  and  directional  hold  response  is  desired,  the  inner  loop  control 
exhibits  integral  response  in  these  loops  as  well.  The  inner  loop  control  laws  require  all 
state"  be  fed  back.  Thus  measurements  of  linear  velocities,  body  angular  rates,  and  roll 
and  pitch  attitudes  are  all  required.  In  practice,  only  the  body  angular  rates  and  the  vertical 
velocity  relative  to  body  axes  are  measured,  therefore  an  estimator  is  required  in  the  inner 
loop.  Various  command  response  types  can  then  be  achieved  by  simple  feedforward  ana 
feedback  outer  loops  wrapped  around  the  inner  loops.  Since  the  inner  loop  feedbacks  are 
flight  critical,  redundancy  [14]  will  be  required  in  these  loops.  The  outer  loops  are  flight 
critical  only  under  severely  degraded  flight  conditions. 
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A .  Inner  Loops 


The  inner  loop  flight  control  system  involves  four  loops.  Vertical  velocity,  roll,  pitch 
and  yaw  rates  are  the  four  sensor  outputs  (  measurements  )  and  collective  pitch,  cyclic 
lateral,  cyclic  longitudinal,  and  tail  rotor  collective  pitch  are  the  four  actuator  inputs  for 
these  loops.  At  hover,  the  mathematical  model  of  the  open  loop  helicopter  exhibits  two 
undamped  unstable  modes  (  forward  velocity  and  side  slip  )  with  frequencies  of 
approximately  0.57  rad/s  and  0.74  rad/s,  and  damping  ratios  C,  of  -0.37  and  -0.048 
respectively.  The  doubling  times  for  forward  velocity  and  side  slip  are  3.27  seconds  and 
19.5  seconds  respectively.  Level  1  handling  qualities  specify  a  minimum  of  6.8  seconds 
for  allowable  time  to  double  amplitude  for  low  frequency  second  order  modes  [15]-fl6]. 


w/6  collective 
r/6  directional 
p/^lateral 
d/^longitudinal 


Fig.  5. 1  Singular  Value  Plot  of  Open  Loop  Helicopter  at  1  T.  n  or. 
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The  open  loop  frequency  response  of  the  helicopter  at  hover  is  shown  in  Fig.  5.i. 
From  this  figure  it  can  be  seen  that  there  are  two  channels,  namely  roll  and  pitch,  with 
singular  values  of  0.0006  and  0.006  at  frequency  about  0.01  rad/s.  This  means  that  closed 
loop  tracking  and  gust  attenuation  will  be  very  poor  in  the  roll  and  pitch  channels  since  low 
frequency  gains  in  these  channels  are  so  small.  The  difference  in  magnitude  between  the 
largest  and  the  smallest  singular  value  is  about  122  dB,  resulting  in  a  condition  number  of 
~106.  Therefore  the  model  is  close  to  being  singular  at  low  frequency.  This  implies  that 
any  attempt  to  provide  compensation  at  low  frequency  by  inverting  the  plant  will  result  in  a 
very  complicated  and  not  very  robust  controller. 


Regulated  variables  for  the  inner  loop  feedbacks  are  vertical  velocity  and  yaw  rates  for 
heave  and  directional  control  and  roll  and  pitch  attitude  rates  for  the  lateral/longitudinal 
control.  The  control  design  was  carried  out  to  obtain  desirable  performance  for  the  four 
regulated  variable  outputs  and  good  stability  robustness  properties  at  the  four  control 
inputs.  For  good  inner  loop  performance,  it  is  necessary  to  have  stable  closed  loop  poles, 
decoupled  closed  loop  responses  and  acceptable  gust  rejection  properties.  Gust  rejection 
and  command  tracking  are  all  improved  with  increasing  bandwidth  [17].  Eigenstructure 
assignment  is  used  to  design  the  inner  loop  control  laws  so  that  decoupled,  first  order  rate 
responses  are  achieved  about  each  axis.  For  a  first  order  system,  the  bandwidth  of  the 
system  is  directly  related  to  the  position  of  the  pole  location.  Therefore  the  pole  location 
associated  with  each  rate  response  is  selected  to  have  sufficiently  high  bandwidth  so  that 
crisp  response  to  command  is  assured.  The  desired  transfer  functions  between  commands 
and  inner  loop  regulated  variables  arc  given  by 


_w _ _ 

wc  (s  +  ^w) 


(5.1) 
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(5.2) 


j>  _  ^p 

Pc  s(s  +  Xp) 

9  _ 

qc  S(s  +  Xq) 

JL  _ 

rc  (s  +  Xr) 


(5.3) 

(5.4) 


The  transfer  functions  between  roll  and  pitch  attitudes  and  lateral  and  longitudinal 
commands  exhibit  1.0/s  characteristics  for  frequencies  below  crossover  and  Xpq/s^ 
characteristics  above  crossover  frequency  as  shown  in  Figure  5.2. 

Under  the  worst  usable  cue  environment.  Level  2  handling  qualities  are  acceptable  for 
the  collective  and  directional  loops.  Since  Level  1  handling  qualities  are  not  required  under 
such  a  situation,  Level  2  handling  qualities  will  suffice.  Level  2  handling  qualities  [17] 
require  the  response  type  [18]  to  be  either  RCDH  (  Rate  Command  Direction  Hold  )  or 
RCHH  (  Rate  Command  Height  Hold  ).  These  response  types  require  open  loop  transfer 
functions  between  collective  and  directional  commands  and  altitude  and  direction  as  shown 
in  Fig.  5.2. 

When  the  loop  between  the  roll  attitude  and  lateral  command  (  given  in  Eq.  5.2  )  is 
closed  with  an  external  pilot  command,  0C,  the  response  that  the  pilot  sees  will  have  integral 
characteristics  as  shown  in  Fig.  5.2.  The  control  law  and  the  kinematic  relationship 
(  between  roll  rate  and  roll  angle  )  enable  the  lateral  loop  to  have  high  gain  at  low 
frequency  and  low  gain  at  high  frequency.  This  satisfies  the  need  for  a  k/s  [9]-[13] 
attitude-to-stick  frequency  response  in  the  region  of  pilot  crossover.  The  high  inner  loop 
gain  at  low  frequency  provides  good  tracking  of  attitude  commands  and  gust  attenuation. 
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while  roll  off  at  high  frequency  helps  to  avoid  potential  instabilities  due  to  high  frequency 
rotor  and  structural  dynamics.  The  same  comments  apply  to  the  pitch  attitude  control  loop. 


FREQUENCY  <rad/s> 

Fig.  5.2  Integral  Characteristics  of  the  Desired  Transfer  Functions. 


The  desired  transfer  functions  between  forward  velocity  and  pitch  command  and  side 
velocity  and  roll  command  are 


u(s) 

qc(s) 

v(s) 

Pc(s) 


s(s  +  ?iu) 

_ 

s(s  +  >.v) 


(5.5) 

(5.6) 


The  eigenvalues  Xu  and  Xw  are  associated  with  the  linearized  drag  forces  in  the 
forward  and  side  directions.  Since  the  values  of  and  Xv  are  small,  the  lateral  and 
longitudinal  velocities  to  roll  and  pitch  commands  are  essentially  Xpq/s2.  Equations  5.1- 
5.6  are  then  used  to  define  the  desired  closed  inner  loop  eigenstructure  for  the  model. 


B .  Outer  Loop 


-93- 


Once  the  inner  loop  system  has  been  designed  to  approximate  the  desired  transfer 
functions  given  in  Eq.  5. 1-5.6,  various  command  response  types  can  be  achieved  by 
simple  single  loop  feedbacks  or  feedforwards. 

The  desired  pitch  and  roll  attitude  transfer  functions  for  ACAH  are  [2],  [  1 8]-[  19] 

<t>  w2a 

—  =  - ® -  (5  7) 

<!>c  (  s2  +  2Ccd({)s  +  co^  ) 

_0_  _  co2 0 _ 

ec  ( s2  +  2Cco0s  +  ci)2e )  (5'8) 

Similar  second  order  transfer  functions  are  desired  for  TRCPH. 

III.  Inner  Loops  via  Eigenstructure  Assignment 
A.  Full  State  Feedback 
1 .  Nominal  Design 

A  full  state  regulator  is  developed  using  eigenstructure  assignment  techniques.  Even 
though  good  decoupling  of  the  closed  loop  eigenvectors  is  achieved,  input  distribution 
matrix  coupling  of  the  helicopter  is  so  great  that  a  feedforward  gain  matrix  is  required. 
Both  the  feedback  control  law  and  the  feedforward  gain  matrix  are  designed  assuming  full 
state  feedback.  Since  not  all  the  states  can  be  measured,  the  feedback  control  law  cannot  be 
implemented  directly.  Thus  it  is  necessary  to  include  a  state  estimator  in  the  feedback  loop 
in  order  to  realize  the  control  law.  The  state  estimator  is  synthesized  using  eigenstructure 
assignment  techniques  so  as  to  achieve  loop  transfer  recovery  in  the  crossover  frequency 
range. 
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The  theory  of  eigenstructure  assignment  by  feedback  control  is  discussed  in  Chapter  3. 
This  section  will  be  limited  to  a  discussion  of  the  design  procedure  and  a  presentation  of  the 
results. 


The  feedback  control  law  is  a  linear  function  of  the  state  vector 


u  =  -Kx 


(5.9) 


The  feedback  gain  matrix,  K  is  designed  such  that  the  implementation  of  this  control  law 
results  in  desired  placement  of  the  closed  loop  eigenvalues  and  shaping  of  the 
corresponding  eigenvectors.  The  placement  of  the  closed  loop  eigenvalues  is  necessary  in 
order  to  obtain  stability  and  set  the  time  constants  and  natural  frequencies  for  the  rigid  body 
modes.  Eigenvector  shaping  is  used  to  obtain  modal  decoupling  needed  to  achieve 
handling  quality  specifications  and  reduce  pilot  workload.  The  values  of  the  closed  loop 
eigenvalues  and  eigenvectors  must  be  specified  by  the  designer.  The  ideal  closed  inner 
loop  state  equations  for  the  helicopter  model  are  derived  from  the  desired  transfer  functions 
(  Eqs.  5. 1-5.6  ).  The  equation  is  given  in  state  space  form  for  the  desired  response  as 


x  =  Ad  x  +  Bd  Xc 


(5.10) 


where 


Ad  = 


0 

o  -K 
o  o 
o  o 
o  o 
0  0 
0  0 
0  0 


0  0 

0  1 

■^•w  0 

0  '^-p 

0  0 

0  0 

0  1 

0  0 


1  0 

0  0 

0  0 

0  0 

'^■q  0 

0  -^r 

0  0 

1  0 


0  ^-q 

^•p  0 

0  0 

0  0 

0  0 

0  0 

0  0 

0  0 
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Bd  = 


0  0  0  0 
0  0  0  0 


Xw  0  0  0 

0  Xp  0  0 

0  0  Xq  0 

0  0  0  XT 

0  0  0  0 

0  0  0  0- 


It  can  be  seen  that  this  state  space  representation  gives  the  transfer  functions  in  Eqs. 
5. 1-5.6.  The  values  of  ^.W)P>qir  for  the  ideal  A  matrix  are  determined  by  the  handling 

quality  requirements.  For  a  first  order  system,  the  bandwidth  and  the  eigenvalue  are 
identical.  Bandwidths  of  4  rad/s  for  vertical  velocity,  roll,  pitch  and  yaw  rate  were  chosen 
to  assure  Level  1  handling  qualities.  The  values  of  ^.u  and  A.v  were  set  at  10%  of  their 
nominal  values  from  the  residualized  A  matrix  since  previous  experience  indicated  that 
these  values  resulted  in  good  steady  state  response  [19].  The  values  of  -.00199  and 
-.00526  for  A.u  and  Xw  are  small  enough  so  that  Xp.q/s2  characteristics  for  lateral  and 
longitudinal  responses  to  roll  and  pitch  commands  are  guaranteed. 


Ed  = 


-.0020 

-.0001 

-.0053 

-.0001 

-4.0000 

-4.0000 

-4.0000 

-4.0000 


Table  5.1a  Desired  Closed  Loop  Eigenvalues 
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With  the  closed  loop  eigenvalues  specified,  this  Ad  is  then  used  to  define  the  desired 
eigenvalue/eigenvector  configurations  which  are  given  in  Table  5.1a. 


Vd  = 


u 

1.0000 

.0000 

.0000 

1.0000 

.0000 

.0000 

.0000 

.0000 

V 

.0000 

1.0000 

1.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

w 

.0000 

.0000 

.0000 

.0000 

1.0000 

.0000 

.0000 

.0000 

p 

.0000 

.0000 

.0000 

.0000 

.0000 

.9701 

.0000 

.0000 

q 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.9701 

.0000 

r 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

r\r\r\r\ 

.Ouuu 

1.0000 

<t> 

.0000 

.0013 

.0000 

.0000 

.0000 

-.2425 

.0000 

.0000 

0 

.0000 

.0000 

.0000 

.0005 

.0000 

.0000 

-.2425 

.0000 

Table  5.1b  Desired  Closed  Loop  Eigenvectors 

The  desired  eigenvectors  associated  with  the  closed  loop  eigenvalues  are  given  in  Table 
5.1b.  Kinematic  constraints  necessitate  the  coupling  of  pitch  rate  to  pitch  angle  and  pitch 
attitude  to  forward  velocity  for  the  longitudinal  eigenvectors.  These  same  conditions  also 
apply  to  roll  rate  being  coupled  to  roll  angle  and  roll  attitude  to  side  velocity.  This  set  of 
eigenvectors  will  achieve  the  required  modal  decoupling  and  therefore,  are  used  as  the 
desired  eigenvectors  for  the  regulator  design  procedure.  The  achievable  non-dimensional 
eigenvectors  are  shown  in  Table  5.2a.  The  closed  inner  loop  system  exhibits  substantially 
reduced  modal  decoupling  when  compared  with  the  open  loop  system. 
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Va  = 


u 

.9992 

-.0007 

-.0006 

.9993 

.0260 

.0025 

-.1063 

0.0000 

V 

-.0006 

.9998 

.9998 

-.0007 

-.0339 

.0467 

-.0077 

-.0244 

w 

-.0010 

-.0002 

.0001 

-.0010 

.9991 

.0015 

.0025 

-.0008 

p 

-.0001 

0.0000 

0.0000 

0.0000 

.0015 

.9691 

.0006 

.0020 

q 

0.0000 

0.0000 

-.0001 

0.0000 

.0025 

.0006 

.9646 

-.0022 

r 

0.0000 

0.0000 

0.0000 

0.0000 

-.0008 

.0011 

-.0002 

.9997 

.0360 

.0138 

.0061 

.0360 

-.0004 

-.2423 

-.0003 

.0033 

0 

.0148 

.0168 

.0168 

.0120 

-.0006 

-.0002 

-.2410 

-.0080 

Table  5.2a  Achievable  Non-Dimensional  Regulator  Eigenvectors 


K  = 


.0248 

.0239 

-1.1628 

-.0092 

.0787 

-.1040 

.1210 

.0735 

-.0553 

-.0355 

.0976 

-.0054 

-.1162 

-.0365 

-.6114 

-.0598 

.0062 

-.0133 

-.0428 

.0209 

.0760 

.1779 

.0028 

.0036 

-.5350 

.0112 

.0096 

-.0073 

-.1664 

-.3627 

-.0205 

.0085 

Table  5.2b  Non-Dimensional  Regulator  Feedback  Gain  Matrix 

The  extensive  coupling  of  the  control  input  matrix,  B,  requires  a  feedforward  gain 
matrix  to  alleviate  the  problem  of  control  cross-coupling.  This  requires  the  desired  control 
distribution  matrix  be  determined.  The  structure  and  the  magnitude  of  the  desired  control 
matrix  should  be  the  same  as  that  of  the  ideal  B  matrix  in  order  to  obtain  the  desired  transfer 
functions.  As  can  be  seen  from  the  desired  control  distribution  matrix  B^  in  Eq.  5.8,  the 
desired  output  to  a  collective  input  is  a  pure  heave  command.  Similarly,  pure  roll  rate  and 
pitch  rate  are  the  desired  commands  to  lateral  and  longitudinal  stick  pitch  inputs.  Also,  a 
pure  yaw  rate  command  is  desired  for  the  pedal  input.  The  feedforward  gain  matrix,  H, 
was  selected  such  that 
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B  *  H  =  Bd 


(5.11) 


Since  the  control  distribution  matrix  B  is  non-square,  it  is  necessary  to  solve  for  H  using 
the  optimal  pseudo-inverse  method  [21],  which  is  optimal  in  the  sense  that 
trace[  (  B*BT)(  B*BT)  ]  is  minimized, 

H  =  Bt*(B*Bt)'1  *  Bd  (5.12) 


The  solution  of  this  relation  gives  the  feedforward  gain  matrix,  H  as  shown  in  Table  5.3a. 


H  = 


-1.2824 

-.0007 

-.0035 

.0011 

.1506 

.2130 

.1003 

.1757 

.0771 

.0090 

-.6223 

.0111 

-.6880 

-.0489 

-.1592 

-.3990 

Table  5.3a  Non-Dimensional  Feedforward  Gain  Matrix 


B  *  H  = 


.0949 

.0026 

-.2950 

.0044 

-.1577 

.0418 

-.0346 

-.0946 

3.9915 

.0016 

.0057 

-.0038 

.0016 

3.9996 

.0006 

.0010 

.0057 

.0006 

3.9778 

-.0005 

-.0038 

.0010 

-.0005 

3.9978 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

Table  5.3b  Achievable  Form  of  Non-Dimensional  Bd 

The  achieved  form  for  non-dimensional  Bd  is  shown  in  Table  5.3b.  This  result  achieves 

the  decoupling  desired  for  the  four  control  inputs.  It  also  includes  some  unwanted  input 
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dynamics  in  forward  and  side  velocities  for  heave  command,  forward  velocity  for 
longitudinal  command,  and  side  velocity  for  yaw  rate  command.  Although  the  coupling 
terms  in  the  achieved  form  of  Bd  are  undesirable,  they  are  tolerable  due  to  the  fact  that  the 
magnitude  of  these  terms  as  compared  to  4  of  the  command  input  is  insignificant. 

Full  state  step  responses  for  the  closed  inner  loop  system  were  evaluated  with  the 
inclusion  of  the  feedforward  gain  matrix.  Using  the  dimensional  regulator  and  feedforward 
gain  matrix,  the  results  of  a  simulation  of  the  response  of  the  full  state  model  to  a  5  deg/s 
step  roll  rate  command  are  shown  in  Figs.  5.3  and  5.4.  Both  the  open  loop  and  closed 
inner  loop  responses  are  plotted  on  the  same  scale.  Figs.  5.3  and  5.4  illustrate  the  unstable 
character  of  the  open  loop  responses.  The  responses  are  extensively  coupled.  In  contrast, 
the  heave,  the  longitudinal,  and  the  directional  modes  of  the  closed  inner  loop  system  were 
nicely  decoupled  from  the  lateral  response.  The  full  state  closed  inner  loop  responses  are 
very  similar  to  the  ideal  responses  given  in  Eqs.  5.1  to  5.6.  The  roll  rate  response  behaves 
like  a  first  order  system  with  a  time  constant  of  1/4  second.  The  roll  angle  increases 
monotonically  with  time  and  the  side  velocity  increases  quadratically  with  time.  The  lack  of 
cross  coupling  can  be  seen  in  the  transfer  function  matrix  in  the  appendix.  The  diagonal 
terms  in  the  transfer  matrix  relating  the  outputs  w(s),  p(s),  q(s),  and  r(s)  to  the  inputs 

4 

wc(s),  pc(s),  qc(s),  and  rc(s)  are  essentially  while  the  off-diagonal  terms  are 

essentially  zero.  The  non-zero  vertical  velocity,  longitudinal  and  directional  responses 
were  very  small  compared  with  the  lateral  response. 

The  bode  plot  for  the  roll  angle  to  lateral  command  transfer  function  can  be  seen  in 
Fig.  5.5.  When  the  loop  between  the  roll  angle  and  lateral  command  is  closed  with  an 
external  pilot  command.  Fig.  5.5  can  be  considered  as  the  loop  transfer  function  for  the 
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Fig.  5.4  Full  State  Open  and  Closed  Inner  Loop  Time  Responses  to  a  Step  Roll  Rate 

Command. 


FREQUENCY  <rad/s> 

Fig.  5.5  Bode  Magnitude  of  Roll  Angle  to  Lateral  Command. 

inner  loop  control  law.  This  plot  shows  high  gain  at  low  frequency  for  good  tracking  and 
low  gain  at  high  frequency  for  robustness.  This  also  assures  that  good  attenuation  of  gust 
can  be  achieved  by  the  inner  loop  control  law. 

In  addressing  the  handling  qualities  issues  related  to  the  disturbance  response  of  the 
rotorcraft,  the  significant  factors  are  the  magnitude  and  the  frequency  content  of  the 
helicopter  response  to  turbulence  for  which  the  pilot  must  compensate.  The  frequency 
range  of  pilot  stabilization  is  approximately  0.5  rad/s  to  5  rad/s  [22]. 

Figure  5.6  shows  the  singular  value  plot  of  identity  matrix  minus  the  transfer  function 
matrix  of  the  closed  inner  loop  system.  In  order  to  have  a  realistic  plot,  the  units  for  the 
transfer  functions  from  lateral,  longitudinal,  and  directional  commands  to  vertical  velocity 
have  to  be  properly  scaled.  Since  the  angular  rate  commands  in  the  other  channels  are 
given  in  deg/s  instead  of  rad/s,  the  units  of  the  command  had  to  be  changed  from  rad/s  to 
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deg/s.  This  was  achieved  by  dividing  the  transfer  functions  from  directional,  lateral,  and 
longitudinal  commands  to  vertical  velocity  by  a  factor  of  57.3. 


-  Col  I ec t i ue  Loop 

. .  Directional  Loop 

—————  Longitudinal  Loop 
— —  Lateral  Loop 

Fig.  5.6  Singular  Values  of  [  I  -  GCIosed  Inner  ^  ]. 


For  a  system  to  have  good  decoupling  characteristics,  the  a  [  I  -  Gclosed  Ioop  ]  must  be 
smaller  than  unity  in  the  frequency  range  of  interest.  As  shown  in  Fig.  5.6,  decoupling  of 
the  heave  mode  is  not  satisfied  at  low  frequencies  but  the  controller  has  been  able  to  achieve 
small  gains  over  the  frequency  range  of  interest,  0.1  to  10  rad/s,  in  each  of  the  four  control 
channels.  This  plot  indicates  that  response  in  each  channel  is  the  same  for  frequencies 
greater  than  1 .0  rad/s  as  the  singular  values  are  essentially  the  same  for  *ach  channel.  Fig. 
5.6  shows  improvement,  in  terms  of  the  decoupling  and  bandwidth  criteria  for  each 
controlled  loop  compared  with  the  open  loop  system.  The  responses  in  Fig.  5.6  exhibit 

decoupling.  This  is  also  shown  by  the  transfer  matrix  given  in  the  Appendix  B.  The 
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singular  value  plot  illustrates  essentially  first  order  response  with  a  bandwidth  of  4  rad/s  in 
each  channel.  At  low  frequencies,  some  coupling  exists  because  it  is  impossible  to  achieve 
full  pole/  zero  cancellation.  Since  the  coupling  occurs  below  the  frequency  range  of 
interest,  it  is  tolerable. 


FREQUENCY  <rad/s> 

-  Col lective  Loop 

■  . .  Directional  Loop 

— — *  — Longitudinal  Loop 
<  —  Lateral  Loop 

Fig.  5.7  Singular  Values  of  o  [  G closed  inner  Loop  !• 

Another  mean  to  verify  the  inner  loop  control  law  is  to  look  at  the  singular  values  of 
closed  inner  loop  system,  as  given  by  Gcioscd  loop  in  Eq.  5.13.  Fig.  5.7  shows  the 
singular  values  plot  of  Gdoscd  ,oop.  The  singular  values  are  First  order  with  a  bandwidth  of 

about  4  rad/s,  except  at  frequencies  below  0.1  rad/s.  This  plot  indiactes  that  the  achievable 
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transfer  functions  between  w  and  wc,  p  and  pc,  q  and  qc,  and  r  and  rc  are  very  close  to  the 
desired  transfer  functions.  The  singular  values  at  low  frequencies  are  not  flat  is  because  the 
achievable  transfer  functions  do  not  have  total  pole  and  zero  cancellation.  These  low 
frequency  dynamics  will  not  become  a  concern  because  they  can  be  corrected  by  the  pilot. 
Fig.  5.7  also  shows  that  response  in  each  channel  is  essentially  the  same  for  frequencies 
greater  than  1.0  rad/s. 


FREQUENCY  (RAD/S) 

-  u 

.  v 

-  vv 

Fig.  5.8  Full  State  Closed  Inner  Loop  Response  of  Accelerations  to  Vertical  Gust. 
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5.9a.  Inner  Loop  Feedback  system  with  Vertical  Gust  to 
Acceleration  Response. 
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5.9b.  Estimator  in  Feedback  Loop. 


Fig.  5.9  Block  Diagram  for  Inner  and  Outer  Loops  Feedback  Configuration. 


A  plot  of  the  transfer  function  response  of  the  three  components  of  the  linear 
acceleration  of  the  center  of  gravity  to  a  vertical  gust  is  shown  in  Fig.  5.8.  The  block 
diagram  of  this  closed  loop  configuration  can  be  seen  in  Fig.  5.9a.  As  would  be  expected, 
longitudinal  and  lateral  accelerations  are  insensitive  to  a  vertical  gust.  In  the  case  of  the 
vertical  acceleration,  the  small  singular  value  at  low  frequencies  [22]  shows  good 
attenuation  of  gust  response.  This  plot  also  shows  that  the  vertical  acceleration  has  a 
bandwidth  of  4  rad/s.  Even  at  frequencies  above  the  bandwidth,  the  gain  between  vertical 
acceleration  and  vertical  gust  is  only  0.4.  Although  no  specifications  exist  for  gust 
response.  Key  [  17]  indicates  that  the  bandwidth  of  the  gust  should  be  about  the  same  as  the 
bandwidth  of  command  response  to  ensure  adequate  disturbance  rejection.  Fig.  5.8  shows 
that  the  closed  inner  loop  system  has  very  good  gust  attenuation  for  all  frequencies. 

B .  Full  State  Estimator 

The  full  state  regulator  in  the  inner  loop  requires  that  all  states  be  fed  back.  Since  only 
vertical  velocity,  pitch  rate,  roll  rate  and  yaw  rate  are  measured  in  the  inner  loop,  it  is 
necessary  to  include  a  state  estimator  in  the  feedback  loop  in  order  to  implement  the  control 
law.  The  use  of  eigenstructure  assignment  techniques  to  achieve  loop  transfer  recovery  is 
discussed  in  Chapter  3.  This  section  presents  a  brief  discussion  of  the  design  procedures 
and  the  results  of  implementing  the  control  law. 

The  control  law  with  the  estimator  in  the  inner  loop  is  given  by 

u  = -  K  x  (5.15) 

where  x  is  the  estimate  of  the  state  vector  given  by  the  state  estimator 

$  =  Ax  +  Bu  +  L(y-Cx)  (5.16) 
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The  block  diagram  of  the  system  with  the  estimator  in  the  feedback  loop  is  given  in 
Fig.  5.9b.  The  estimator  gain  matrix  L,  is  to  be  selected  such  that  A  -  LC  is  stable  and  the 
frequency  response  of  the  system  with  the  regulator/estimator  in  the  feedback  loop 
approximates  the  full  state  loop  transfer  matrix.  That  is 

K8(s)G8(s)  «  [  K  (  si  -  A  )-l  B  ](s)  (5.17) 

where  K8(s)  =  [  K  (  si  -  A  +  BK  +  LC  f1  L  ](s)  and 
G8(s)  =  [  C  (  si  -  A  )-l  B  ](s) 

This  is  called  loop  transfer  recovery  (  LTR  )  [23]-[24], 

The  four  finite  open  loop  transmission  zeroes  of  the  transfer  matrix  relating  the  outputs 
w(s),  p(s),  q(s),  and  r(s)  to  the  inputs  wc(s),  pc(s),  qc(s),  and  rc(s)  of  the  helicopter  are 

Zj2  =  0.0 

z34  =  -.000734  ±. 01 398i 

The  four  desired  eigenvalues  and  eigenvectors  of  the  estimator  were  placed  at  the  four 
transmission  zero  positions  and  their  respective  left  zero  directions.  The  values  of  zj  and 
Z2  were  approximated  by  -.001  and  -.0015  because  of  computational  difficulties  arising 
from  multiple  roots.  The  remaining  four  closed  loop  eigenvalues  were  placed  at  -10.0  rad/s 
(  Table  5.4  )  and  their  corresponding  eigenvectors  were  arbitrarily  chosen.  The  desired  and 
achievable  estimator  eigenvectors  are  given  in  Table  5.5.  Figure  5.10  shows  the  full  state 
and  the  combined  regulator/estimator  inverse  closed  loop  transfer  functions.  This 
configuration  achieves  full  LTR  from  0.1  rad/s  up  to  5  rad/s.  This  frequency  range  is 
desirable  because  it  includes  the  inverse  closed  loop  transfer  function  characteristics  of  the 
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bandwidth  of  the  regulator.  Below  and  above  this  frequency  range,  the  presence  of  the 
estimator  in  the  feedback  loop  results  in  an  additional  20  dB/decade  of  roll  off.  The  high 
frequency  roll  off  enhances  robustness  compared  with  the  full  state  regulator. 
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Full  State  Inner  Loop 
Estimator  in  Inner  Loop 


Fig.  5.10  5  [  I  +  (  K(s)G(s)  )_1  ]  Full  State  Inner  Loop  and  Estimator  in  the  Feedback 

Loop. 
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ED  = 


-.0007  ± 

.0140i 

-.0010  + 

.0000i 

-.0015  + 

.OOOOi 

-10.0000  + 

,0000i 

-10.0000  + 

.OOOOi 

-10.0000  + 

.OOOOi 

-10.0000  + 

.OOOOi 

Table  5.4  Desired  Estimator  Eigenvalues 


VD  =  (  Desired  estimator  eigenvectors ) 


u 

.0083  ± 

.0 1 52i 

0.0000 

0.0000 

.0000 

.0000 

.0000 

.0000 

V 

.01 19  ± 

,0006i 

0.0000 

0.0000 

.0000 

0000 

.0000 

.0000 

w 

.0003  ± 

.0003i 

0.0000 

0.0000 

1.0000 

.0000 

.0000 

.0000 

p 

-.0001  ± 

O.OOOOi 

0.0000 

0.0000 

.0000 

1.0000 

.0000 

.0000 

q 

.0007  ± 

.001  li 

0.0000 

0.0000 

.0000 

.0000 

1.0000 

.0000 

r 

.0003  ± 

O.OOOOi 

0.0000 

0.0000 

.0000 

.0000 

.0000 

1.0000 

0 

-.0006  ± 

.5642i 

.7690 

.7761 

.0000 

.0000 

.0000 

.0000 

e 

-.7032  ± 

.432 1  i 

-.6393 

-.6307 

.0000 

.0000 

.0000 

.0000 

VA  =  ( Attainable  estimator  eigenvectors  ) 


u 

.0083  ± 

.01 52i 

-.0010 

-.0014 

.0078 

-.0446 

-.0366 

-.1094 

V 

.0119  ± 

.0006i 

-.0012 

-.0018 

.0070 

.2510 

-.0182 

-.0693 

w 

.0003  ± 

.0003i 

0.0000 

0.0000 

.9999 

-.0016 

.0004 

.0014 

p 

-.0001  ± 

O.OOOOi 

0.0000 

0.0000 

-.0015 

.9667 

.0031 

.0133 

q 

.0007  ± 

.001  li 

0.0000 

0.0000 

.0004 

.0032 

.9991 

-.0054 

r 

.0003  ± 

O.OOOOi 

-.0001 

-.0001 

.0014 

.0136 

-.0054 

.9914 

0 

-.0006  ± 

.5642i 

.7690 

.7761 

-.0027 

-.0167 

.0012 

.0046 

e 

-.7032  ± 

.432 1  i 

-.6393 

-.6307 

-.0005 

-.0030 

-.0024 

-.0073 

Table  5.5  Desired  Estimator  Eigenstructure 
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L  = 


.0139 

-.0307 

-2.0074 

.0055 

-.0212 

-.0276 

2.0934 

-2.0215 

.5700 

.0018 

-.0133 

.1357 

-.0045 

20.0522 

-2.1698 

1.7284 

-.0103 

.2049 

26.9162 

-.0361 

-.0021 

1.1677 

1.0953 

27.0865 

0.0000 

2.8616 

.0368 

-.0648 

0.0000 

-.0043 

2.9044 

.0744 

Table  5.6  Non-Dimensional  Estimator  Gain  Matrix 

The  non-dimensional  estimator  gain  matrix  is  given  in  Table  5.6.  Improved  LTR  is 
obtained  as  the  finite  eigenvalues  are  moved  closer  to  the  transmission  zeroes  and  as  the 
infinite  eigenvalues  are  moved  to  the  left.  It  is  important  to  place  the  infinite  eigenvalues  no 
further  to  the  left  than  is  required  to  achieve  LTR  necessary  for  stability  robustness,  since 
moving  the  infinite  eigenvalues  further  to  the  left  increases  the  estimator  bandwidth  and 
sensor  noise  may  become  a  problem. 

The  frequency  response  of  the  inner  loop  compensator  (  regulator/  estimator  ),  [  K  (  si 
-  A  +  BK  +  LC  )_1  L  ](s)  is  shown  in  Fig.  5.1 1.  It  can  be  seen  that  the  controller  has  a 
condition  number  which  is  of  the  same  order  as  the  plant.  This  seems  reasonable  because 
the  open  loop  model  has  a  large  separation  in  the  singular  values  at  low  frequency  with  the 
roll  and  pitch  channels  being  very  small  and  the  vertical  velocity  and  yaw  rate  channels  very 
large.  In  order  to  achieve  good  tracking  and  gust  rejection  characteristics,  it  is  necessary  to 
increase  the  gain  in  each  loop  proportionally  so  that  the  loop  gain  of  K(s)G(s)  as  defined  in 
Eq.  5.17  has  an  integral  characteristic  at  frequencies  below  crossover  and  1/s2  characteristic 
above  crossover  frequency.  Therefore  the  increases  in  tne  gains  in  the  roll  and  pitch  loops, 
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5  longitudinal 

.  5  lateral 

“ “  ^directional 

”  ^collective 

Fig.  5.1 1  Singular  Value  Plot  of  Compensator  (  Regulator/  Estimator ). 

which  have  a  very  low  open  loop  DC  gains,  will  have  proportionally  greater  increases  in 
gains  than  in  the  vertical  velocity  and  yaw  rate  loops,  which  have  higher  open  loop  DC 
gains.  This  requirement  is  indicated  by  the  singular  value  plot  of  the  compensator,  K(s)  in 
Fig.  5.11. 

IV.  Outer  Loop 

Once  the  inner  loop  system  has  been  design  to  approximate  the  desired  transfer 
functions  given  in  Equations  5. 1-5.6,  various  command  response  types  can  be  achieved  by 
simple  single  loop  feedbacks  or  feedforwards.  The  inner  loop  system  already  exhibits  RC 
response,  so  only  ACAH  and  TRCPH  are  discussed. 


A.  ACAH 
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For  ACAH,  the  lateral  and  longitudinal  angles,  0  and  0,  can  be  fed  back  by  an  outer 
loop  with  a  proportional  gain  of  and  Kq  to  achieve  Attitude  Command/  Attitude  Hold 
type  response  about  the  roll  and  pitch  axes.  Thus  only  the  roll  and  pitch  angle 
measurements  are  needed  for  the  outer  loop  design.  If  the  outer  loop  sensors  fail,  the  inner 
loop  control  controller  will  maintain  stability  of  the  helicopter,  therefore  multiple  redundant 
sensors  are  not  required  in  the  outer  loop. 


The  transfer  functions  between  commanded  and  actual  pitch  and  roll  angle  are  given  in 
Eq.  5.18-5.19.  Both  and  Kq  were  selected  as  2.  The  transfer  functions  are  second 
order  with  a  natural  frequency,  ci)nA0  =  2.83  rad/s  and  a  damping  ratio,  £  =  0.707 


_0_  _  _ _ 

0c  (  s2  +  \pS  +  2A.p  ) 

9  _  _ KgXq _ 

0C  (  s2  +  XqS  +  2Xq  ) 


(5.18) 

(5.19) 


Fig.  5.12  Definition  of  Phase  Margin  Limited  Bandwidth. 
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Using  the  handling  quality  specification  for  a  second  order  system,  the  desired 
bandwidth  is  defined  in  terms  of  the  natural  frequency  and  damping  ratio  of  the  system  as 
in  Eq.  5.20.  This  type  of  bandwidth  is  phase  margin  limited  and  is  defined  as  the 
frequency  where  the  phase  angle  becomes  -135  degrees  as  shown  in  Fig.  5.12.  The  roll 
and  pitch  attitude  command  response  give  a  Level  1  handling  quality  with  a  bandwidth  of 
5.46  rad/s. 

®bw  =  con(C  +  Vc2  +  1  )  (5.20) 

Since  direct  force  control  is  available  only  in  the  heave  direction,  the  lateral  and 
longitudinal  velocities  of  the  helicopter  in  hover  must  be  controlled  using  roll  and  pitch 
moments  about  the  aircraft  center  of  gravity.  The  relationship  between  pilot  commands  a 
change  in  roll  or  pitch  attitude  and  lateral  or  longitudinal  velocities  is  shown  by  the 
following  transfer  functions 


u(s) 

Kg^q  (s  +  Xq) 

(5.21) 

0C(S) 

(s  +  A.u)  (s^  +  A,qS  +  2)tq) 

v(s) 

(s  +  X.p) 

(5.22) 

(S  +  Xy)  (s^  4-  XpS  +  2Xp) 

The  resulting  closed  roll  and  pitch  attitude  outer  loops  eigenvalues  are  obtained  as 

-.0006  ±  ,0140i 

-2.0107  ±  1.9866i 

-1.9925  ±  2.0038i 

-4.0000  +  .0000i 

-4.0000  +  .0000i 

As  can  be  seen  from  the  closed  outer  loop  poles,  the  eigenvalues  corresponding  to  vertical 
velocity  mode,  roll  and  pitch  attitudes,  and  yaw  rate  were  very  much  as  specified.  There  is 
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a  low  frequency,  oscillatory  mode  with  a  frequency  of  con  =  0.014  rad/s  and  a  damping 
ratio  of  C,  =  0.043.  This  phugoid  period  mode  might  become  unstable  due  to  uncertainties 
in  the  aerodynamic  coefficients.  The  amount  of  instability  which  can  be  tolerated  can  be 
estimated  by  examining  the  minimum  allowable  time  to  double  amplitude  for  a  particular 
mode.  The  equation  for  time  to  double  amplitude  is  given  as 

tdouble  =  if  X  is  real  (5.23) 

or 

^double  =  if  Cco,,  is  the  real  part  of  the  complex  root  (5.24) 

W 

For  ACAH  response  type,  the  specified  minimum  tdoub]c  is  6.8  seconds.  Under  the  worst 
case  perturbations  in  the  A  and  B  matrices,  it  can  be  shown  [25]  that  this  phugoid  mode  has 
a  time  to  double  amplitude  of  tdou},ie  =  22.5  seconds.  This  time  is  long  enough  that  the 
pilot  can  instinctively  correct  for  any  instability  with  his  normal  control  motions.  There  is 
some  indication  that  many  pilots  prefer  a  slightly  unstable  phugoid  mode. 

The  results  of  a  simulation  of  the  response  of  the  helicopter  to  a  step  roll  rate  command 
of  5  deg/s  with  only  the  inner  loop  closed  and  to  a  step  roll  attitude  command  of  5  degrees 
with  both  the  inner  and  outer  loops  closed  are  shown  in  Figs.  5.13  and  5.14.  The  actual 
helicopter  responses  shown  in  Figures  5.13  and  5.14  are  very  similar  to  the  ideal  responses 
given  in  Eqs.  5. 1-5. 6,  Eqs.  5.18-5.19  and  Eqs.  5.21-5.22.  With  only  the  inner  loop 
closed,  the  roll  rate  response  behaved  like  a  first  order  system,  and  with  both  the  inner  and 
outer  loops  closed,  the  roll  angle  exhibited  a  second  order  response  as  desired.  It  can  be 
seen  from  Figs.  5.13-5.14  that  the  inner  loop  gives  rate  response  to  pilot  input  and  the 
outer  loop  gives  attitude  response.  Also,  the  decoupling  of  the  states  is  further  improved 
when  the  outer  loop  is  closed.  An  effect,  called  "rate  cross-coupling",  which  holds  for  a 
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Closed  Inner  Loop 
Closed  Outer  Loop 


Fig.  5.14  Full  State  Closed  Inner  and  Outer  Loop  Time  Responses  to  a  Step  Roll 

Command. 


8 


Fig.  5. 15  Full  State  Closed  Outer  Loop  Response  to  a  5  degrees  Step  Roll  Command  of 

Duration  4  seconds. 
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right  (  positive  )  roll,  i.e.  the  helicopter  will  pitch-up  when  rolling  right,  is  shown  in  the 
inner  and  outer  loop  responses.  The  reason  why  the  outer  loop  pitch  rate  returns  to  zero  is 
that  once  a  steady  roll  attitude  is  established,  there  is  no  longer  any  rolling  angular  velocity 
and  this  particular  effect  disappears.  Also,  the  side  velocity  for  the  outer  loop  is  essentially 
linear  after  1  second  but  the  side  velocity  response  for  the  inner  loop  crosses  the  side 
velocity  response  for  the  outer  loop  at  about  2.5  seconds  and  increases  quadratically.  This 
implies  that  in  order  to  keep  a  constant  roll  rate,  higher  side  thrust  force  is  required. 

The  decoupled  character  of  the  response  is  further  ilias^ated  in  Fig.  5.15  which  shows 
the  response  to  a  5  deg  roll  angle  command  which  is  maintained  for  four  seconds  and  then 
returned  to  zero.  The  results  were  plotted  such  that  all  linear  velocities  and  angular 
displacements  are  on  the  same  graphs  to  better  illustrate  the  closed  loop  decoupling.  The 
response  of  the  side  velocity  to  roll  command  results  in  an  acceleration  command/  velocity 
hold  type  of  response.  Once  the  attitude  achieves  its  commanded  value,  the  control  is 
returned  to  the  zero  stick  position.  The  attitude  drops  back  to  its  original  value  within  two 
seconds  after  the  control  is  released,  whereas,  the  side  velocity  remains  almost  constant. 
Fig.  5.15  also  shows  that  after  the  stick  has  been  released,  there  is  a  slight  increase  in  the 
forward  velocity.  Pitch  attitude  remains  unchanged  at  zero  through  out  the  10  seconds 
indicating  that  the  decoupling  is  good. 

Control  actuator  deflections  were  shown  in  Fig.  5.16  for  the  maneuver  illustrated  in 
Figs.  5.13-5.14.  It  can  be  seen  that  the  maneuver  requires  less  than  one  degree  deflection 
from  collective,  lateral,  longitudinal,  and  tail  rotor  pitch.  These  values  were  very  small 
compared  to  the  maximum  control  input  values  of  9  degrees  for  collective,  8.75  degrees  for 
lateral,  15  degrees  for  longitudinal,  and  18.5  degrees  for  collective  tail  rotor.  Hence  the 
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inputs  will  not  saturate  the  system  with  a  aggressive  maneuver  that  requires  a  ±30  deg.  roll 
angle  command.  This  would  be  a  very  large  roll  angle  command  for  a  helicopter  in  hover. 


6  collective 

.  5  lateral 

^longitudinal 
—  Sdirectional 

Fig.  5.16  Control  Actuator  Deflections. 

Closed  outer  loop  time  responses  were  also  calculated  for  a  5  ft/s  step  command  in 
vertical  velocity  (  Figs.  5.17  and  5.18  )  and  a  5  deg/s  command  in  yaw  rate  (  Figs.  5.19 
and  5.20 ).  Pilot  collective  and  pedal  commands  are  of  primary  importance  in  maintaining 
hover.  Considerable  training  and  coordination  are  required  for  pilots  to  maintain  a 
helicopter  at  a  constant  altitude  during  hover.  Changes  in  collective  result  in  changes  in 
torque  produced  by  the  engine,  this  in  turn  produces  yawing  moments  which  must  be 
compensated  for  by  the  changes  in  tail  rotor  collective  pitch.  These  changes  in  tail  rotor 
collective  pitch  effect  the  power  supplied  to  the  main  rotor  which  requires  changes  in  main 
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Fig.  5.17  Full  State  Closed  Inner  and  Outer  Loop  Time  Responses  to  a  5  ft./s  Heave 

Command. 
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Fig.  5.19  Full  State  Closed  Inner  and  Outer  Loop  Time  Responses  to  a  5  deg./s  Yaw  Rate 

Command. 
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Closed  Inner  Loop 


Closed  Outer  Loop 


Fig.  5.20  Full  State  Closed  Inner  and  Outer  Loop  Time  Responses  to  a  5  deg./s  Yaw  Rate 

Command. 


rotor  collective  pitch  and  so  on.  Decoupling  response  to  collective  and  pedal  inputs  results 

in  improved  handling  qualities.  Figures  5.17  and  5.18  illustrate  the  response  to  a  collective 

input  from  the  pilot.  A  collective  pitch  change  from  the  pilot  input  gives  a  change  in  the 

mean  blade  angle  of  attack,  which  produces  a  change  in  the  thrust  magnitude.  The  change 

in  thrust  magnitude  increases  the  main  rotor  torque,  which  in  term  requires  the  increase  in 

the  tail  rotor  thrust  power  to  provide  the  anti-torque  forces  to  balance  the  helicopter.  There 

is  a  non  negligible  side  velocity  which  is  caused  by  the  change  in  tail  rotor  pitch  required  to 

balance  the  main  rotor  torque.  Fig.  5.18  shows  that  the  side  velocity  increases  linearly  in 

order  to  keep  a  constant  yaw  rate.  Correction  of  this  side  velocity  would  require  a  small 

lateral  stick  input.  It  can  be  seen  that  the  yaw  rate,  roll  rate,  pitch  rate,  roll  angle  and  pitch 

angle  are  all  negligible.  The  qualitative  requirement  given  for  collective  control  is  to 

achieve  vertical  velocity  response  and  shall  not  cause  objectional  response  in  yaw  or  more 

other  axes  [18].  Limits  on  this  form  of  cross-coupling  are  specified  in  terms  of  yaw  rate 

r,  r, 

oscillation  parameters  per  rate  of  climb.  The  values  for  —  and  —  are  0.00111  and 

h3  h3 

0.001 14  respectively.  These  values  are  much  smaller  than  the  values  of  0.65  and  0.15  as 
required  by  Level  1  handling  qualities. 

Figures  5.19  and  5.20  show  the  effect  of  a  5  deg./s  yaw  rate  command.  The  change 
in  vertical  velocity  is  negligible  however  the  lateral  and  longitudinal  responses  are  small  but 
not  negligible  and  would  require  both  lateral  and  longitudinal  stick  inputs  for  stabilization. 
During  yawing  maneuvers,  the  drag  force  of  a  tail  rotor  must  be  countered  by  the  main 
rotor.  The  decoupling  of  all  states  is  improved  when  the  outer  loop  is  closed. 

Since  pilots  are  sensitive  to  the  shape  of  the  phase  curve  at  frequencies  above  the 
bandwidth  frequency,  phase  delay  parameter,  xp  as  defined  in  Eq.  5.25  is  used  for  the 
evaluation  of  the  design. 
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Fig.  5.21  Phase  Delay  of  ACAH  System  plus  Effective  Time  Delay. 
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Fig.  5.22  Level  1  Small-Amplitude  Requirements  for  Minimum  and  Maximum  Gain 

Tasks. 
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Fig.  5.23  Closed  Outer  Loop  Full  State  Roll  Attitude  Bode  Magnitude  rind  Phase  Plots, 


AC>2colsn 
57.3  (  2(0]go  ) 


(5.25) 


Hoh  [19]  indicates  that  for  large  values  of  phase  delay,  the  phase  curve  drops  off  more 
rapidly  than  for  small  values.  The  evaluation  of  xp  is  based  on  the  assumption  that  the 
ACAH  system  has  a  total  effective  time  delay  of  0.15  seconds  in  the  roll  axis.  The 
resulting  bode  phase  plot  about  the  roll  axis  of  the  ACAH  system  plus  total  effective  time 
delay  is  shown  in  Fig.  5.21.  From  figure  5.21,  A02(0igo  is  determined  to  be  57.57 
degrees  and  2coigo  *s  determined  to  be  10.98  rad/s.  Using  Eq.  5.25,  the  phase  delay  of  the 
overall  system  is  calculated  to  be  xp  =  0.092  seconds.  As  shown  in  Fig.  5.22,  this  value 
satisfies  Level  1  handling  qualities  requirements  for  ultra  high  gain  tasks. 

The  bode  magnitude  plot  in  Fig.  5.23  shows  that  the  roll  attitude  response  to  lateral 
input  is  flat  until  about  2.83  rad/s,  it  then  rolls  off  at  40  dB/decade.  The  bandwidth  is 
determined  to  be  5.46  rad/s  from  the  phase  plot  in  Fig.  5.23  using  the  phase  margin  limited 
criteria.  Similar  results  are  obtained  for  the  pitch  transfer  function.  The  vertical  velocity 
and  yaw  rate  transfer  functions  are  characterized  by  simple  decoupled  first  order  responses 
and  are  similar  to  each  other. 


B.  TRCPH 


With  the  inner  rate  loops  closed,  the  transfer  function  between  the  forward  velocity 
and  the  rate  command  is 


(5.26) 


The  desired  transfer  function  for  forward  velocity  TRCPH  is 

u  _  _ _ 

uc  s2  +  2£cous  +  Cl)2u 


(5.27) 
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This  transfer  function  can  be  obtained  as  shown  below 


Fig.  5.24  Block  Diagram  for  Forward  Velocity  TRCPH. 


with  the  equivalent  transfer  function  being 


S“  +  ^  ^q^u 


(5.28) 


From  Eqs.  5.27  and  5.28  with  Xq  fixed  from  the  inner  loop  design,  the  gain  Ku  and  time 
constant  xu  can  be  selected  to  give  the  desired  damping  ratio  and  natural  frequency.  This 


system  will  give  a  constant  forward  velocity  output  for  a  constant  non-zero  pilot  input  and 
zero  velocity  for  zero  input.  If  more  precise  position  control  is  desired,  it  can  be  achieved 
by  the  addition  of  a  positional  feedback  loop  as  shown 


Fig.  5.25  Block  Diagram  for  Position  Control. 


Similar  result  can  be  obtained  for  control  of  TRCPH  of  lateral  velocity  and  position. 
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Chapter  6 


Robustness  Analysis 

I.  Introduction 

In  control  system  design,  the  physical  system  to  be  controlled  is  represented  by  a 
nominal  mathematical  model  derived  from  assumptions  concerning  the  physics  governing 
the  system.  The  control  system  is  then  designed  and  evaluated  using  this  model.  Due  to 
our  inability  to  mathematically  model  any  system  exactly,  the  physical  system  upon  which 
the  control  system  actually  operates  responds  differently  to  control  inputs  than  its  nominal 
mathematical  model.  Thus,  no  matter  how  well  a  control  law  may  work  when  applied  to 
the  nominal  model,  it  will  work  differently  when  applied  to  the  actual  system  because  of  the 
presence  of  model  uncertainties.  The  effects  of  these  modeling  errors  must  be  evaluated 
before  the  control  law  can  be  considered  for  implementation.  This  evaluation  is  called 
robustness  analysis.  The  robustness  of  SISO  systems  has  been  traditionally  determined 
using  the  concepts  of  gain  and  phase  margins.  Robustness  of  MIMO  systems  requires 
more  sophisticated  concepts,  and  has  been  the  area  of  intensive  research  for  the  last  decade 
[1]. 

The  purpose  of  this  Chapter  is  to  use  several  tools  to  analyze  the  robustness  of  the 
control  laws  designed  in  Chapter  5.  The  methods  used  are  unstructured  and  structured 
singular  value  analyses  and  simulauons  using  mathematical  models  which  are  more  realistic 
than  the  nominal  design  model.  The  properties  which  are  of  most  interest  are  :  (1)  nominal 
performance,  that  is  how  well  the  control  system  meets  performance  specifications  when 
applied  to  the  model  used  for  control  law  design,  (2)  stability  robustness,  the  ability  of  the 
control  system  to  maintain  stability  despite  the  presence  of  model  uncertainties,  and  (3) 
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robust  performance,  the  ability  of  the  control  system  to  meet  the  performance  specifications 
and  maintain  stability  in  the  presence  of  model  uncertainties  [2].  Nominal  performance  was 
evaluated  in  Chapter  5  so  this  Chapter  will  concentrate  primarily  on  the  evaluation  of 
stability  and  performance  robustness. 

II.  Unstructured  Singular  Value  Analysis 

A.  Stability  Robustness 

Since  the  inner  loop  system  is  flight  critical,  the  stability  robustness  of  this  stabilization 
loop  with  respect  to  modeling  errors  and  parameter  variations  is  very  important.  A 
sufficient  condition  [3]  for  the  closed  inner  loop  system  to  be  robustly  stable  with  respect  to 
the  multiplicative  error  E(s)  [4],  is 

2  ( I  +  [  K(s)G(s)  l*1 )  >  a  (E(s))  foralls=jco  (6.1) 

where  0  <  co  <  °°. 

As  discussed  in  Chapter  4,  the  multiplicative  error  is  modeled  as  an  effective  time 
delay.  The  stability  boundary  is  determined  by  varying  the  effective  time  delay  in  the 
closed  inner  loop  system.  Many  components  of  the  effective  time  delay  such  as  delays  due 
to  rotor  response  and  structural  dynamics  are  essentially  fixed;  however,  others  such  as 
computational  delays  and  sampling  times  may  be  influenced  by  the  design  of  the  flight 
control  system  hardware  and  software.  For  example,  a  faster  computer,  more  effective 
algorithms  or  a  faster  sampling  rate  may  be  used  to  reduce  actuator  delays. 

Based  on  the  data  in  Reference  5,  a  break  down  of  the  various  delays  is  given  in  Table 

6.1. 
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Element 

Delay  ( ms ) 

Rotor 

66 

Actuators 

31 

Zero-order  hold 

17 

Computations 

22 

Notch  filter 

11 

Total  delay 

147 

Table  6.1  Summary  of  Time  Delay  Contributors. 


The  total  effective  time  delay  is  chosen  to  be  0.15  seconds.  Since  large  relative 
modeling  errors  cannot  be  avoided  as  frequency  increases,  the  modeling  error  is  estimated 
to  lie  within  the  region  shown  in  Fig.  6.1.  If  the  modeling  error  is  indeed  within  this 
region,  stability  of  the  closed  loop  system  is  assured  if  £  (  I  +  [  K(s)G(s)  ]_1  )  does  not 
enter  this  region.  The  time  delay  is  modeled  by  a  first  order  Pade  approximation.  The 
error  is  then 


E(s)  =  ( 


-  2Ats 
Ats  +  2 


)I 


(6.2) 


The  maximum  bandwidth  that  the  closed  loop  system  can  achieve  results  when  the 

-  2  A  T  S 

error  model  has  a  magnitude  of  one,  i.e.  1 - 1  =  1.  Letting  s  =  jo),  the  maximum 

Axs  +  2 

bandwidth  of  the  first  order  closed  loop  system  ( i.e.  a  loop  transfer  function  of  L  =  ) 


can  be  related  to  total  effective  time  delay  of  the  uncertainty  model  as  follows 


wmax  bw 


2 

V" 3  At 


(6.3) 


The  above  equation  provides  a  simple  quantitative  relationship  useful  in  preliminary 
design  studies  as  it  allows  the  designer  to  trade  off  various  hardware  and/or  software 
options  versus  bandwidth.  As  shown  in  Eq.  6.3,  the  achievable  closed  loop  crossover 
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8 


frequency  depends  inversely  on  the  effective  time  delay  Ax  in  the  stabilization  loop.  Using 
this  equation,  the  maximum  achievable  crossover  frequency  is  estimated  to  be  7.7  rad/s, 
which  is  substantially  higher  than  the  design  bandwidth  of  4  rad/s.  However,  it  would  not 
be  prudent  to  try  and  achieve  the  maximum  achievable  bandwidth,  because  the  error  model 
is  itself  only  an  estimate  and  approximation  of  the  actual  unknown  error. 


FREQUENCY  (RAD/S) 


-  a  (  E  )  (  Effective  Time  Delays,  x  =  0. 1 5  sec.  ) 

.  Full  State  Inner  Loop 

Fig.  6.1  <2  [  I  +  (  K(s)G(s)  l*1  ]  Full  State  Inner  Loop  and  Estimator  in  the  Feedback 

Loop. 

As  can  be  seen  in  Fig.  6.1,  the  stability  condition  given  in  Eq.  6.1  is  violated  below 
0.01  rad/s.  Condition  6.1  is  conservative  so  that  the  helicopter  may  not  be  unstable; 
however,  this  is  expected  because  of  the  high  condition  number  at  low  frequency.  The 
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frequency  of  0.01  rad/s  is  much  lower  than  the  closed  loop  bandwidth  of  4  rad/s  so  even  if 
the  helicopter  exhibits  low  frequency  instabilities,  it  should  not  have  any  significant  effect 
on  the  handling  qualities. 


FREQUENCY  (RAD/S) 

-  a(E)  (  Effective  Time  Delays,  i  =  0. 1 5  sec.  ) 

.  Full  State  Closed  Outer  Loop 

Fig.  6.2  Minimum  Singular  Value  for  Full  State  Closed  Outer  Loop  and  Maximum 
Singular  Value  of  Error  Matrix  vs.  Frequency. 

Figure  6.2  illustrates  that  the  stability  robustness  of  the  helicopter  is  further  improved 
when  the  outer  loop  is  closed.  With  the  outer  loop  closed,  the  roll  and  pitch  attitude 
responses  become  a  second  order.  Closing  the  outer  loop  improves  the  stability  robustness 
of  the  system  at  frequencies  below  0.01  rad/s.  A  measure  of  the  robustness  of  a  MIMO 
system  is  the  minimum  difference  between  the  minimum  singular  value  of 
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[  I  +  (  K(s)G(s)  )_1  ]  and  the  maximum  singular  value  of  the  multiplicative  error,  E(s). 
This  can  be  regarded  as  a  MIMO  gain  margin.  As  shown  in  Fig.  6.2,  the  minimum 
difference  between  the  two  curves  is  6.7  dB  at  about  7.39  rad/s. 

B .  Nominal  Performance 


Fig.  6.3  Performance  Test  for  Closed  Inner  Loop  System. 


The  plot  of  the  maximum  singular  value  of  a  closed  loop  frequency  response  matrix 

can  be  used  to  test  performance  requirements  for  disturbance  rejection  and  command 

following.  For  the  closed  inner  loop  system  shown  in  Fig.  6.3,  it  is  desired  to  achieve 

-  w  - 

to  the  output  q  .  The  performance 
.  r  . 

4 

specification  can  then  be  represented  by  a  Bode  magnitude  plot  of  r(s)  =  g  If  the 

closed  inner  loop  gain  a[  G(  I  +  K  G  )'*  ]  from  the  inputs  to  the  outputs  remains  near  the 
curve  r(s)  at  all  frequencies,  the  system  meets  its  specification.  A  sufficient  test  [3]  is  to 
require  that 

— - —  O  [  G(  jco  )  ( I  +  K(  jco )  G(  jco  )  r1  ]  <  1  for  all  co  e  SK 

I  r(jco)  I 

(6.4) 

A  plot  of  the  maximum  singular  value  of  the  closed  inner  loop  system  multiplied  by  the 
inverse  of  the  performance  specification  is  shown  in  Fig.  6.4.  The  plot  shows  that  the 


+  4  type  of  response  from  the  input 
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system  does  not  achieve  the  desired  -  —  ^  transfer  function  requirement  at  frequencies 


below  1  rad/s,  however,  the  deviation  from  1  is  not  large  for  frequencies  greater  than 
0.1  rad/s.  This  further  verifies  the  results  obtained  in  Chap.  5  that  the  closed  inner  loop 
system  does  not  have  total  poles/  zero  cancellation  at  low  frequencies.  However,  the  time 
histories  given  in  Chapter  5  show  that  the  response  of  the  system  to  commands  is  excellent. 


Fig.  6.4  — - -  d[G(jCo)(I  +  K(jco)G(jco))1]. 

I  r(jco)  I 


III.  Structured  Singular  Value  Analysis 


A.  Theory 


Unstructured  singular  value  analyses  are  sometimes  overly  conservative  because  they 
do  not  take  into  account  any  information  which  may  be  available  on  the  structure  of  the 
uncertainties.  Structured  singular  value  analyses  methods  have  been  developed  to  provide 
less  conservative  estimates  of  robustness  properties  than  can  be  obtained  from  unstructured 


singular  value  analyses. 
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This  section  gives  a  fairly  brief  introduction  to  structured  singular  value,  or  |i-based 
methods.  The  p.-based  methods  discussed  here  have  proven  to  be  useful  for  analyzing  the 
performance  and  robustness  properties  of  linear  feedback  systems.  The  general  robust 
performance  framework  used  was  developed  by  Doyle  [7]-[  12]  in  the  early  1980's.  Doyle 
uses  the  general  interconnection  structure  illustrated  in  Fig.  6.5 


Fig.  6.5  General  Interconnection  Structure. 


Any  linear  interconnection  of  inputs,  outputs,  commands,  perturbations  and  a 
controller  can  be  rearranged  in  the  form  of  Fig.  6.5.  The  nominal  model  P(s),  a  3x3 
partitioned  matrix,  the  controller  K(s)  and  the  perturbations  A(s)  are  assumed  throughout  to 
be  linear  and  time  invariant.  Four  of  the  signals  are  familiar  from  standard  control  system 
literature.  These  are  external  inputs  (  disturbances  ),  d,  errors,  e,  measured  plant  outputs 

y,  and  control  commands  u.  These  four  signals  form  four  pairs  of  terminals,  one  being  the 

V  .  c  . 

jj  pair  across  which  the  control  law  is  connected,  one  being  j  representing  errors  due  to 

inputs  or  disturbances,  one  being  representing  sensor  outputs  due  to  inputs,  and  the 

e 

other  is  the  -  pair  which  represents  the  errors  due  to  command  inputs.  The  remaining 

terminals  are  connected  to  a  bounded  perturbation  model  (  A  )  representing  various 
uncertainties  in  the  system. 
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The  robust  performance  problem  for  helicopter  control  is  defined  in  the  block  diagram 
shown  in  Fig.  6.6.  The  eighth  order  rigid  body  dynamic  equations  are  in  the  helicopter 
equations  block,  which  is  equivalent  to  the  nominal  model  P(s)  in  Fig.  6.5.  Sensors 
measure  linear  velocities,  body  angular  rates,  and  roll  and  pitch  attitudes.  Model 
uncertainty  is  parameterized  with  a  multiplicative  perturbation  at  the  control  inputs.  The 
controller  K(s)  consists  of  full  state  feedback,  feed  forward,  and  outer  loop  feedback  gain 
matrices. 


vertical  velocity  error 
roll  attitude  error 
pilch  attitude  error 
yaw  rate  error 


Fig.  6.6  Robust  Performance  Block  Diagram  Description. 

The  interconnection  structure  shown  in  Fig.  6.6  has  four  inputs  and  four  outputs.  The 
four  inputs  correspond  to  the  pilot's  collective,  lateral  cyclic,  longitudinal  cyclic  and  tail 
rotor  collective  inputs.  The  four  outputs  correspond  to  vertical  velocity  error,  roll  and  pitch 
attitude  errors,  and  yaw  rate  error. 

For  the  purpose  of  analysis,  the  diagram  of  Fig.  6.6  reduces  to  that  of  Fig.  6.7  which 
is  analogous  to  closing  the  lower  loop  of  P  with  K  in  Fig.  6.5.  The  feedback  control  law 
can  be  thought  of  as  another  system  component  and  which  is  incorporated  into  the  model  to 
yield  a  stable  closed  loop  system,  M.  This  forms  the  closed  loop  response  interconnection 
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structure  suitable  for  (i  analysis.  The  closed  loop  system  has  a  frequency  response  M(jco) 
that  depends  on  P(jco)  and  K(jco). 


A  mult. 


coll,  stele 
lat.  cyclic 
long,  cyclic 
TR  pedals 


vertical  vel.  error 
roll  attitude  error 
pitch  attitude  error 
yaw  rate  error 


Fig.  6.7  Perturbed  Pilot  Stick  Inputs  to  Errors. 


B .  Mathematical  Formulation 


This  section  begins  with  an  explanation  of  the  matrix  and  block  diagram  notation  that 
will  be  used  throughout  this  Chapter.  9TX*  is  a  real  n  x  k  matrix.  Let  M  e  'Rnxk,  a  real 
rational  matrix  function  of  the  Laplace  operator,  s,  whose  elements  can  be  represented  by 
rational  transfer  functions,  u  and  v  are  real  vectors,  with  u  e  9t\  v  e  91",  and  v  =  M  u  as 
given  below 


Fig.  6.8  Matrix -Vector  Multiplication. 


The  block  diagram  for  the  standard  feedback  loop  is  illustrated  in  Fig.  6.9  where  G  is 
the  plant  and  K  is  the  controller.  It  is  conventional  to  absorb  any  scalings,  weights,  or 
coloring  filters  into  G  to  form  P  as  given  in  Fig.  6.5.  Since  G  is  unweighted  for  the 
analysis  purpose,  it  is  the  same  as  P  in  Fig.  6.5.  If  the  eighth  order  dynamic  equations  of 
the  helicopter  in  Fig.  6.6  are  defined  to  be  G,  and  the  full  state  feedback,  feed  forward,  and 
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outer  loop  feedback  gain  matrices  are  defined  to  be  the  controller  K,  then  Fig.  6.6  is 
equivalent  to  the  standard  feedback  loop. 


Fig.  6.9  Block  Diagram  for  the  Standard  Feedback  Loop. 


In  order  to  introduce  a  class  of  general  linear  feedback  loops  called  the  Linear 
Fractional  Transformation,  the  matrix  M  is  partitioned  in  2  x  2  fashion  as 


Fig.  6.10  Partitioned  Matrix-Vector  Multiplication. 


Now,  consider  the  analysis  model  which  isolates  the  nominal  part  of  the  system  M 
from  the  perturbations  A  as  shown  in  Fig.  6.10.  Suppose  there  is  a  defined  block  structure 
A  which  is  compatible  in  size  with  Mn.  Let  A  =  diag  (  Aj, ... ,  A„  )  be  a  block  diagonal 
matrix.  Let  A  be  the  set  of  all  matrices  with  a  given,  fixed  block  diagonal  structure.  The 
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individual  i  blocks  may  be  euher  scalars  or  matrices.  Assume  that  A,  e  that  is, 

n 

each  A,-  is  a  square  &,•  x  ki  complex  matrix.  The  total  matrix  dimension  of  A  is  m  =  k  ,• . 

i  =  1 

It  is  known  that  for  a  block  diagonal  A  e  A,  a  (  A  )  =  max,  a  (  A,  )  [12].  That  is,  the 
overall  gain  of  A  is  equal  to  the  gain  of  its  largest  diagonal  block. 

For  A  e  A,  consider  the  following  loop  equations 

z  =  Mn  v  +  M12  d 

e  =  M2j  v  +  M22  d  (6.6) 

v  =  A  z 

The  elements  of  M  can  be  defined  for  the  standard  feedback  loop  shown  in  Fig.  6.9  with 
performance  and  plant  perturbations  specified  at  output  and  input  respectively,  as 


M|j  =  KG(I  +  KG)'1 

(6.7) 

M12  =  Kd  +  GK)1 

(6.8) 

M21  =  G(I  +  KG)‘ 

(6.9) 

M22  =  (I  +  GK)'1 

(6.10) 

The  set  of  equations  (6.6)  is  well  posed,  if  for  any  vector  d,  there  exist  unique  vectors 
v,  z,  and  e  satisfying  the  loop  equations.  Substituting  v  into  the  equation  for  e,  it  is 
obvious  that  the  set  of  equations  is  well  posed  if  and  only  if  the  inverse  of  [  I  -  Mn  A  ] 
exists.  If  not,  then  depending  on  d  and  M,  there  is  either  no  solution  to  the  loop  equations, 
or  there  are  an  infinite  number  of  solutions.  If  the  inverse  of  ( I  -  Mu  A  )  does  exist,  the 
transfer  function  from  d  to  e  can  be  expressed  as  e  =  Fu  ( M,A )  d  where 

Fu  ( M,A )  =  M22  +  M21  A(I-MU  A)'1  M12 
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(6.11) 


Fu  ( M,A  )  is  called  a  Linear  Fractional  Transformation  on  M  by  A,  and  in  a  feedback 
diagram  appears  as 


Fig.  6.1 1  Linear  Fractional  Transformation. 


As  can  be  seen  from  Fig.  6.1 1,  M22  is  the  nominal  map  between  the  disturbance  and 
error,  and  A  represents  perturbations,  which  affect  the  map  through  Mn,  M12,  and  M21. 
The  subscript  «onf  means  that  the  "upper"  loop  of  M  is  closed  by  A.  An  analogous 
"lower"  linear  fractional  transformation  F\  ( P,K  ),  describes  the  resulting  matrix  M 
obtained  by  closing  the  "lower"  loop  of  P  in  Fig.  6.5. 

The  structured  singular  value  of  the  matrix  M  e  C myjn  denoted  |i  [  M  ]  is  defined  by 

MM]  :=  - - - - -  (6.12) 

min  {  c  (A)  :  det  (  I  +  M  A  )  =  0  } 

Ae  A 

If  no  Ae  A  makes  I  +  M  A  singular,  and  then  pi  [  M  ]  =  0.  In  addition,  the  structured 
singular  value  has  the  following  properties 

p  [  M  ]  <  (i  [  M  ]  <  d[M]  (6.13) 

where  p  [  M  ]  is  the  spectral  radius  (largest  eigenvalue  )  of  M  [8], 

C.  Robust  Stability 
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The  robust  stability  question  is  cast  in  terms  of  the  closed  loop  response 
M  =  F[(  P,K  )  shown  in  Fig.  6.6  and  it  involves  determining  whether  the  error  remains 
in  a  desired  set  for  the  sets  of  inputs  and  perturbauons  A  of  size  1  or  less.  The  robust 
stability  quesuon  is  answered  by  computing  the  structured  singular  value  of  Mn  for  block 
diagonal  A.  Since  A  has  only  one  block  in  this  case,  the  structured  and  unstructured 
singular  values  are  the  same.  In  general,  this  is  not  true.  If  this  structured  singular  value  is 
less  than  one  for  all  frequencies,  then  stability  is  guaranteed  for  all  possible  uncertainties. 
That  is,  stability  is  robust  if  and  only  if 

a[Mn(jco)]  <1  V  to  (6.14) 

or 

p  [  Mu  (jeo)  ]  <1  V  co  for  block  diagonal  A  (6.15) 

D.  Nominal  Performance 

In  the  helicopter  control  problem,  the  nominal  performance  is  acceptable  if  the 
responses  meet  the  handling  qualities  specifications.  The  Level  1  handling  qualities  for  a-' 
ACAH  system  require  the  following  transfer  functions.  Bandwidth  must  be  2  rad/s  or 
greater 

_w _  ^w 

WC  (s  +  Xw) 

±  _  _ CO^ _ 

<l)c  (  s2  +  2Cco4)s  +  co2^  ) 

Q_  _  _ 0)^9 _ 

0C  (  S2  +  2£cOqS  +  (O20  ) 


T_ _ V 

rc  (S  +  Xr) 

Nominal  performance  means  that  with  no  model  uncertainty  as  illustrated  in  Fig.  6.12, 
the  system  satisfies  the  above  transfer  functions  from  inputs  d,  to  outputs  e.  If  M22 
achieves  nominal  performance,  the  closed  outer  and  inner  loop  system  will  satisfy  the 
handling  qualities  specifications.  The  nominal  performance  question  is  answered  by 
computing  the  maximum  singular  value  of  M22.  If  thij  maximum  singular  value  is  less 
than  one  for  all  frequencies  then  nominal  performance  has  been  achieved. 


Fig.  6.12  Linear  Fractional  Transformation  for  Nominal  Performance. 

The  nominal  performance  specification  is  satisfied  if  and  only  if 
a[M22(jco)]  <1  V  co  (6.16) 

E.  Robust  Performance 

The  performance  is  robust  if  the  desired  response  meets  the  handling  qualities 
specifications  (  performance  )  in  the  presence  of  nonzero  model  uncertainty.  The  robust 
performance  question  is  answered  by  computing  the  structured  singular  value  of  M  for  a 
block  diagonal  perturbation  structure.  If  the  structured  singular  value  of  M  is  less  than  one 
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for  all  frequencies,  then  robust  performance  has  been  achieved.  That  is  performance  is 
robust  if  and  only  if 

(i  [  M  (jco)  ]  <  1  V  co  (6.17) 

F.  Weighting  Functions 

The  structured  singular  value  analysis  will  be  discussed  with  reference  to  the 
equivalent  block  diagrams  shown  in  Figs.  6.6  and  6.7.  The  frequency  response  of  the 
perturbation  is  used  to  model  the  effective  time  delay,  e_sAT.  Using  a  first  order  Pade 
approximation  for  the  total  effective  time  delay,  the  error  can  be  expressed  as 

A(s)  =  [  ~  2S2  ]I  (6.18) 

s  +  — 

At 

The  uncertainty  is  formulated  in  such  a  way  that  the  perturbations  enters  the  system  in  a 
feedback  form,  or  a  linear  fractional  way.  For  the  uncertainty  model  to  be  more  realistic,  a 
10%  error  at  low  frequencies  is  included  in  the  error  model.  In  the  case  of  the 
multiplicative  model  uncertainty  parameterized  at  the  control  inputs,  the  formalism  for 
robust  . stability  analysis  requires  I  Amull  ( jco  )  I  <  1  for  all  co.  For  Ax  =  0. 15  seconds.,  the 
unweighted  transfer  function  of  Aj,^  is 

Amuk(s)  =  (  s  +  0.001  )  (  ~  2  )  (6.19) 

s  +  0.15 

This  transfer  function  will  have  magnitude  greater  than  one  at  frequencies  greater  than  7.7 
rad/s.  Therefore,  a  weighting  function  for  the  multiplicative  uncertainty  is  required  to  make 
I  Amuit  ( jco  )  I  <  1  for  all  co.  This  weighting  transfer  function,  wA  (s),  is  selected  such  that 

Amuit  ( jco  )  =  A'mull  (jco  )  wA(jco)  (6.70) 
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where 


IA'mult(jco)!<l  Vcd  (6.21) 

Hence  the  weighting  transfer  function  for  the  multiplicative  uncertainty  is  given  by 

WaOO  =  (  sVo°001  >  < - ^2 — >  = 

s  +  0.15 

The  weighting  function  for  the  inputs  was  chosen  such  that  maximum  pilot  input  to 

each  channel  satisfies  the  ACAH  handling  qualities  specifications  [  1 3]-[  1 4]  for  small 

amplitude  inputs.  The  mathematical  model  [  15)-[  17]  of  the  human  pilot's  central  nervous 

and  neuromuscular  systems  is  very  complicated.  Therefore,  a  simple  pilot  model  of  a  pure 

gain  with  a  representative  transport  delay  is  used  (  high  frequency  pilot  neuromuscular 

dynamics  are  included  in  the  time  delay  )  [  1 8]-[  19],  The  addition  of  pilot  lead  is  also 

studied  by  Aponso  [18J.  References  14  to  18  show  that  a  first  order  low  frequency  phase 

lag  characteristics  are  expected  when  pure  gain  dynamics  are  being  controlled.  The 

nominal  design  in  Chapter  5  can  be  considered  as  a  unity  gain  system.  For  simplicity,  the 
pilot  is  modeled  as  a  unity  gain  low  pass  filter  of  the  form  — — — .  Based  on  experimental 

S  T  s 

results  [20],  a  bandwidth  of  0.1  rad/s  was  chosen  for  this  filter. 

From  Reference  13,  a  maximum  input  of  15  ft./s  in  vertical  velocity,  ±  10°  in  roll 
attitude,  ±  5”  in  pitch  attitude  and  ±  10  deg./s  in  yaw  rate  are  required  for  small  amplitude 
maneuvers.  The  resulting  input  weighting  function  is 
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r  i5(o.i) 

s  +  0.1 

0 

0 

0 

"  ft./s  ” 

0 

10(0.1) 
s  +  0.1 

0 

0 

deg. 

0 

0 

5(0.1) 
s  +  0.1 

0 

deg. 

0 

0 

0 

10(0.1) 
s  +  0.1  -1 

L-  deg./s  -* 

According  to  Reference  14,  the  desired  performance  for  precision  hover  is  to  maintain 
altitude  within  ±  1  meter  (  3.28  ft. )  and  heading  within  ±  2  degrees.  Similar  requirements 
are  apply  to  roll  attitude  and  pitch  attitude.  We  will  assume  that  errors  must  be  reduced  to 
zero  in  one  second  or  less.  Therefore,  the  output  weighting  function  is  selected  so  that 
vertical  velocity  is  within  ±  3.28  ft./s,  yaw  rate  is  within  ±  2  deg./s,  and  roll  and  pitch 
attitude  angles  are  within  ±  2  degrees.  The  resulting  output  weighting  matrix  is  given  as 


3.28  0  0  0 

-  s/ft  - 

0  i  0  0 

1/deg 

O 

—  ICN 

O 

o 

1/dcg 

1 — 

o 

o 

o 

Nil  •— * 

1 _ 

L  s/deg  -1 

When  the  multiplicative  uncertainty,  input  and  output  weighting  matrices  are  included 
in  M  -  f-  (  P,K  ),  the  interconnection  structure  M'  is 

wA  KG  (  I  +  KG  )_1  wA  K  (  I  +  GK  )-1  wd''  i 


M’  = 


we  G  (  I  +  KG  r1  we(I+GK  r1  wd-> 


(6.22) 


Therefore,  the  robust  stability  test  is 


o[M'n  ]  <  1 


V  to 


(6.23) 


or 
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a[wAKG(I  +  KG)1]  <  1  Va> 

The  nominal  performance  test  is 

o  [  M'22  ]  <1  Vco 


(6.24) 


or 


a  [  we  ( I  +  GK  )A  wd4  ]  <  1  V  co 

and  the  robust  performance  test  is 

(i  [  M’  ]  <  1  Vco  (6.25) 

G.  Results 

A  plot  of  the  maximum  singular  value  of  M'n  is  shown  in  Figure  6.13.  This  indicates 
whether  stability  is  guaranteed  for  model  uncertainty  discussed  in  the  previous  section. 
The  maximum  singular  value  has  a  value  of  greater  than  one  at  frequencies  lower  than  0. 1 
rad/s.  Since  this  is  in  the  frequency  range  in  which  instabilities  will  still  result  in  Level  1  or 
2  handling  qualities,  the  robust  stability  test  is  considered  satisfactory  even  though  there  is 
potential  fc"  low  frequency  instabilities. 
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log  magnitude 


M  (P,  K)  Mull  .max 


log  frequency  ( rad/s  ) 

Fig.  6.13  o  [  Mjj  ],  Maximum  Singular  Value  for  Robust  Stability. 

The  maximum  singular  value  plot  of  M'22  is  shown  in  Fig.  6.14.  Since  it  is  less  than 
one  for  all  frequencies,  the  vertical  velocity,  roll  and  pitch  attitudes,  and  yaw  rate  deviations 
from  desired  are  within  the  errors  admitted  by  the  weighting,  we,  discussed  in  Part  E  of 
this  section.  Thus  nominal  performance  is  guaranteed.  This  result  also  agrees  with  the 
time  history  responses  obtained  in  Chapter  5. 
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log  magnitude 


M(P,K)Mu22.max 


log  frequency  ( rad/s  ) 

Fig.  6. 14  a  [  M 22  ],  Maximum  Singular  Value  for  Nominal  Performance. 

The  robust  performance  question  is  answered  by  computing  the  structured  singular 
value  of  M'  for  a  two  block  structure.  One  block  of  this  two  block  structures  is  Amu][  and 
one  block  is  "unit  sized"  A0  (  nominal  performance  )  [11].  A  plot  of  this  structured 
singular  value  is  shown  in  Figure  6.15.  The  structured  singular  value  is  greater  than  one  at 
frequencies  lower  than  0.15  rad/s.  The  robust  performance  is  considered  acceptable 
because  frequencies  lower  than  0.15  rad/s  are  associated  to  the  low  frequency  "phugoid" 
modes.  This  will  be  discussed  in  detail  the  following  section.  The  vertical  velocity,  roll 
and  pitch  attitude,  and  yaw  rate  responses  deviations  from  desired  are  within  the  errors, 
wc,  in  the  presence  of  model  uncertainty. 
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log  magnitude 


log  frequency  ( rad/s  ) 


Fig.  6.15  |i  [  M'  ],  Structured  Singular  Value  for  Robust  Performance. 

IV.  Simulations 

A .  Analysis  of  Worst  Case  Perturbation 

A  different  stability  robustness  test  was  used  to  examine  the  effect  of  variations  of 
parameters  in  the  nominal  model.  This  is  due  mainly  to  the  difficulty  of  determining  the 
aerocoefficients  of  the  rotorcraft  accurately  and  results  in  low  frequency  uncertainties. 
Reference  21  used  the  stochastic  root  locus  method  to  determine  the  influence  of  the 
simultaneous  variations  of  between  ±  20%  the  nominal  values  of  ail  the  elements  of  the  A 
anu  B  matrices  on  the  closed  inner  loop  eigenvalues.  It  was  determined  that  varying  the 
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element  A4  2  of  the  state  matrix  by  -  20%  and  the  element  B4  2  of  the  control  matrix  by 
+  20%  resulted  in  the  largest  positive  eigenvalue  of  any  combination  of  variations. 
Therefore,  the  goal  of  this  section  is  to  determine  whether  the  controller  is  able  to  maintain 
acceptable  Level  1  handling  qualities  for  the  closed  inner  and  outer  loops  system  with  this  a 
worst  case  variation  of  the  elements  of  the  state  and  control  matrices. 


No  Variation 

Varying  B4  2 

Varying  A4  2  & 
B4.2 

-.0053 

.2801 

.2481 

.3635 

-.0020 

-.3107 

-.2716 

-.4093 

-.0001 

-.0001 

-.0001 

-.0001 

-.0001 

-.0020 

-.0020 

-.0020 

-4.0001 

-3.9749 

-4.2882 

-4.2659 

-4.0000 

-4.0001 

-4.0001 

-4.0001 

-4.0001 

-4.0000 

-4.0000 

-4.0000 

-4.0000 

-4.0000 

-4.0000 

-4.0000 

Table  6.2  Comparison  of  Closed  Inner  Loop  Eigenvalues  of  the  Nominal  Design  with 

Parameter  Variations. 

Table  6.2  shows  the  eigenvalues  obtained  from  three  different  case  studies  and 
compares  them  to  the  eigenvalues  of  the  nominal  design.  The  maximum  eigenvalue  on  the 
right  half  plane  that  will  satisfy  Level  2  mid-term  handling  qualities  is  0.139,  therefore, 
with  the  worst  case  parameter  variation  of  the  state  and  control  matrices,  the  closed  inner 
loop  system  mid-term  response  will  be  degraded  to  Level  3. 

Table  6.3  shows  the  variations  in  eigenvalues  when  the  outer  loops  are  closed  with 
gains  in  the  roll  and  pitch  attitude  loops.  The  phugoid  mode,  though  unstable  is  low 
frequency  and  can  be  stabilized  with  little  effort  by  the  pilot.  These  eigenvalues  do  not 
satisfy  Level  1  limit  on  roll  (  pitch  )  oscillation.  As  given  in  Chapter  4,  Level  1  limits  for 
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tasks  which  do  not  require  significant  division  of  attention  requires  a  minimum  time  to 
double  amplitude  of  6.76  seconds  and  a  damping  ratio  £  of  -0.2.  The  above  three  cases 
have  tdouble  of  38.3  seconds  (  varying  A42  ),  55.6  seconds  (  varying  B4  2  )  and  25.2 
seconds  (  varying  A4  2  &  B4  2  simultaneously  )  which  satisfy  the  doubling  time 
requirement,  but  their  damping  ratios  £  of  -0.8,  -0.51  and  -0.96  are  beyond  the  limit  of 
Level  1.  However,  these  values  are  adequate  for  Level  2  handling  qualities. 


No  Variation 

Varying  A4  2 

Varying  B4  2 

Varying  A4  2  & 
B4.2 

-.0006  ±  .01 40i 

.0180  ±  .01211 

.0124  ±  .021  li 

.0274  ±  ,0084i 

-2.0107  ±  1.9866i 

-2.0118  ±  2.019H 

-2.1602  ±  2.2786i 

-2.1754  ±  2.2958i 

-1.9925  ±2.0038i 

-2.0101  ±  1.9933i 

-2.0094  ±  1.9932i 

-2.0091  ±  1.9932i 

-4.0000 

-4.0000 

-4.0000 

-4.0000 

-4.0000 

-4.0000 

-3.9997 

-3.9996 

Table  6.3  Comparison  of  Closed  Outer  Loop  Eigenvalues  of  the  Nominal  Design  with 

Parameter  Variations. 

The  results  of  a  simulation  of  the  response  of  the  helicopter  to  a  step  roll  attitude 
command  of  5  degrees  with  both  the  inner  and  outer  loops  closed,  which  includes  a  -  20% 
variation  of  A4  2  in  the  state  matrix  are  shown  in  Figs.  6.16  and  6.17.  The  responses  of 
the  system  that  includes  -  20%  variation  of  A42  are  very  similar  to  the  responses  with  no 
variation  in  A4  2  for  the  first  two  seconds.  After  two  seconds,  the  responses  in  vertical 
velocity,  roll  and  pitch  attitudes,  and  yaw  rate  diverge,  however,  cross  coupling  is  still 
small  and  performance  acceptable.  As  given  by  the  doubling  time  and  also  shown  in  Figs. 
6.16  and  6.17,  the  rate  of  divergence  is  slow  and  should  be  able  to  be  stabilized  by  the  pilot 
with  little  effort.  Similar  responses  are  shown  in  Figs.  6.18  to  6.21  for  the  case  with  20% 
variation  of  B4>2  of  the  control  matrix  and  varying  A4  2  and  B42  simultaneously.  Again, 
performance  is  satisfactory  even  with  these  variations. 
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Closed  Outer  Loop 

Closed  Outer  Loop  with  Variation  i  n  A(  4,2)  &  B(  4,2) 


Fig.  6.21  Full  State  Closed  Inner  and  Outer  Loop  Time  Responses  to  a  5  Degrees  Step 
Roll  Command  with  Parameter  Variation  of  -  20%  in  i  and  +  20%  in  B* 


B .  Time  Histories  with  High  Order  Model 


Mode 

l  V  l 

-.0006  ± 

.01 39i 

-.0006  ± 

.0 1 40i 

Pitch  Attitude 

-2.4448  ± 

2.423 1  i 

-2.0107  ± 

1.9866i 

Roll  Attitude 

-2.2628  ± 

2.3004i 

-1.9925  ± 

2.0038i 

Yaw  Rate 

-3.8038  + 

•0000i 

-4.0000  + 

,0000i 

Vertical  Velocity 

-3.9995  + 

,0000i 

-4.0000  + 

,0000i 

-13.7551  ± 

71.7906i 

-10.2323  ± 

3.6497i 

Table  6.4  Comparison  of  Closed  Inner  and  Outer  Loop  Eigenvalues  of  12th  Order  and  8th 
Order  Models  Using  the  Same  Full  State  Regulator. 

The  full  state  regulator  designed  for  the  eighth  order  nominal  system,  was  evaluated 
using  the  12th  order  model  which  includes  rotor  dynamics.  The  closed  inner  and  outer 
loop  eigenvalues  are  shown  in  Table  6.4.  The  eigenvalues  corresponding  to  the  tip  path 
plane  flapping  dynamics  are  essentially  the  same  as  the  open  loop  system  since  there  is  no 
compensation  for  these  two  modes  from  the  full  state  regulator.  The  remaining  eigenvalues 
for  the  12th  order  system  are  very  close  to  the  8th  order  system,  therefore,  similar  time 
responses  are  expected  for  the  12th  and  8th  order  models. 

The  12th  order  closed  inner  and  outer  loop  time  responses  were  also  calculated  for  a 
5  ft./s  step  command  in  vertical  velocity  (  Figures  6.22  to  6.23  ),  a  5  degree  step  command 
in  roll  and  pitch  attitude  (  Figures  6.25  and  6.26,  and  Figures  6.28  and  6.29  ),  and  a  5 
deg./s  step  command  in  yaw  rate  (  Figures  6.31  to  6.32  ).  The  longitudinal  and  lateral 
rotor  flapping  angles  and  angular  velocities  for  the  four  maneuvers  are  shown  in  Figs. 
6.24,  6.27,  6.30,  and  6.33. 
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1  2th  Order  Closed  Outer  Loop 


Fig.  6.22  12th  Order  and  8th  Order  Closed  Loop  Responses  to  a  5  ft./s  Step  Vertical 

Velocity  Command. 
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8th  Order  Closed  Outer  Loop 
1  2th  Order  Closed  Outer  Loop 


Fig.  6.23  12th  Order  and  8th  Order  Closed  Loop  Responses  to  a  5  ft./s  Step  Vertical 

Velocity  Command. 
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8th  Order  Closed  Outer  Loop 
1 2th  Order  Closed  Outer  Loop 


Fig.  6.25  12th  Order  and  8th  Order  Closed  Loop  Responses  to  a  5  Degrees  Step  Roll 

Attitude  Command. 
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Fig.  6.27  12th  Order  and  8th  Order  Closed  Loop  Responses  to  a  5  Degrees  Step  Roll 

Attitude  Command. 


K 


Figures  6.22  and  6.23  illustrate  the  response  to  a  collective  command  of  5  ft./s  from 
the  pilot.  It  can  be  seen  that  the  vertical  velocity  response  is  first  order  with  little  coupling 
in  the  other  seven  states.  In  general,  the  responses  of  the  12th  order  model  in  Figures  6.22 
and  6.23  are  very  close  to  the  8th  order  full  state  responses  with  the  exception  of  the 
overshoot  for  the  first  two  seconds  in  pitch  angle,  pitch  rate,  roll  angle,  roll  rate  and  yaw 
rate.  This  is  expected  because  one  of  the  flapping  modes  has  very  small  damping  effects 
and  the  full  state  regulator  does  not  have  any  effect  on  these  modes.  Fig.  6.23  shows  the 
responses  of  the  lateral  and  longitudinal  tip  path  plane  tilting  angles  and  angular  velocities. 
These  include  the  motions  of  the  blades  which  have  very  low  damping.  The  angular 
velocities  settled  down  to  zero  after  one  second,  the  tilting  angles  are  small  and  should  not 
affect  the  performance  of  the  helicopter. 

Figures  6.25  to  6.26  show  the  roll  attitude  response  resulting  from  pilot  lateral 
command  of  5  degrees.  The  response  of  roll  angle  to  roll  attitude  command  was  second 
order  and  settled  out  to  approximately  5  degrees  in  three  seconds.  The  responses  of  the 
12th  order  model  in  Figures  6.25  and  6.26  are  close  to  the  8th  order  responses  with 
overshoot  in  pitch  angle,  pitch  rate  and  yaw  rate  for  the  first  two  seconds.  The  coupling 
between  the  pitch  attitude  and  roll  attitude  command  is  smaller  than  the  full  state  model. 
The  responses  of  the  three  linear  velocities  are  identical  up  to  5  seconds  for  both  full  state 
and  12th  order  models.  Figure  6.27  shows  that  the  flapping  angular  velocities  settle  down 
to  approximately  zero  after  one  second.  The  longitudinal  tilting  angle  is  essentially  linear 
after  two  seconds  and  is  increasing  at  a  rate  of  0.47  deg./s.  Since  the  12th  order  model 
with  the  full  state  regulator  is  a  stable  system,  this  linear  response  will  eventually  reach  a 
steady  state  value. 


-  171  - 


O  UJ  O  ^  tfl  -0-  O  UJ  o 


1.000 


r 


.5000  -f 


-  .500  4—* - 1 

0  5 


TIME  (SEC) 


P 

D 

E 

G 

/ 

S 


TIME  (SEC) 


0  5 

TIME  (SEC) 


8th  Order  Closed  Outer  Loop 
1  2th  Order  Closed  Outer  Loop 


Fig.  6.28  12th  Order  and  8th  Order  Closed  Loop  Responses  to  a  5  Degrees  Step  Pitch 

Attitude  Command. 
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Fig.  6.29  12th  Order  and  8th  Order  Closed  Loop  Responses  to  a  5  Degrees  Step  Pitch 

Attitude  Command. 
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Fig.  6.31  12th  Order  and  8th  Order  Closed  Loop  Responses  to  a  5  deg./s  Step  Yaw  Rate 

Command. 
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Fig.  6.33  12th  Order  and  8th  Order  Closed  Loop  Responses  to  a  5  deg./s  Step  Yaw  Rate 

Command. 


Figures  6.28  and  6.29  show  the  pitch  attitude  response  resulting  from  the  pilot 
longitudinal  command  of  5  degrees.  The  decoupling  in  all  axes  are  essentially  the  same  as 
the  8th  order  model.  Similar  overshoots  in  roll  attitude,  roll  rate  and  yaw  rate  are  seen  in 
the  plots.  There  is  small  coupling  between  the  pitch  attitude  command  and  the  lateral  tilting 
angle  as  illustrated  in  Figure  6.30.  Yaw  rate  response  to  a  pilot  command  of  5  deg./s  is 
shown  in  Figures  6.31  and  6.32.  The  small  lag  in  the  vertical  velocity  response  is  expected 
due  to  higher  order  dynamics  of  the  12th  order  model.  Good  decoupling  in  the  longitudinal 
and  lateral  flapping  modes  is  shown  in  Figure  6.33. 
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APPENDIX  A 

DIMENSIONAL  12TH  &  8TH  ORDER  A  AND  B  MATRICES 
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Dimensional  12th  Order  Helicopter  A  Matrix 


A12  - 


Columns  1  thru  8 


-.0232 

.0055 

-.0059 

-.0521 

.7174 

.0211 

.0000 

-32.1662 

-.0569 

-.0557 

-.0061 

-.8480 

-.0519 

.7144 

32.1473 

-.0169 

-.0772 

-.0748 

-.3803 

-.0413 

-.3013 

2.0310 

1.1024 

.4938 

-.0219 

-.0339 

-.0039 

-.5980 

-.0431 

.4003 

.0000 

.0000 

.0039 

-.0012 

-.0036 

.0077 

-.1631 

.0029 

.0000 

.0000 

.0293 

.0161 

-.0007 

.1352 

.0423 

-.4930 

.0000 

.0000 

.0000 

.0000 

.0000 

1.0000 

.0005 

-.0154 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.9994 

.0343 

.0000 

.0000 

1.1291 

-.3616 

.0036 

-70.2077 

-34.8969 

.3893 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

-.3414 

-1.0981 

.0058 

-34.5620 

70.2831 

.4193 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

Columns  9  thru  12 


.5483 

-5.0970 

-.0860 

9.7470 

-.0881 

9.7240 

-.5523 

5.0490 

-.0078 

-.9641 

.0288 

-1.0320 

-.0687 

7.9190 

-.4014 

32.8146 

-.0703 

5.8210 

.0111 

-1.2540 

.0377 

-5.0972 

.0573 

-3.4304 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

-27.5297 

-511.8210 

-67.8411 

-941.9460 

1.0000 

.0000 

.0000 

.0000 

67.9287 

935.7810 

-27.2686 

-540.1146 

.0000 

.0000 

1.0000 

.0000 
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Dimensional  12th  Order  Helicopter  B  Matrix 


Bl2  “ 


-4.9320 

-8.8350 

25.4900 

.0004 

-3.9280 

25.5100 

8.7880 

18.7500 

-334.6000 

.6213 

-.2959 

.0133 

-5.3705 

18.5905 

6.4560 

10.0067 

-.6877 

1.1350 

-3.2690 

-.1231 

12.8487 

-2.5707 

-1.0964 

-11.9561 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

-1.7003 

855.9650 

-395.4310 

.1293 

.0000 

.0000 

.0000 

.0000 

4.1285 

380.5095 

849.5440 

-9.9949 

.0000 

.0000 

.0000 

.0000 

Dimensional  12th  Order  Helicopter  C  Matrix 

C12  = 


0. 

0. 

1. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0 

0. 

0. 

0. 

1. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0 

0. 

0. 

0. 

0. 

1. 

0. 

0. 

0. 

0. 

0. 

0. 

0 

0. 

0. 

0. 

0. 

0. 

1. 

0. 

0. 

0. 

0. 

0. 

0 
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Dimensional  8th  Order  Design  Model  A  Matrix 


A8  - 


-.0199 

-.0058 

-.0058 

-.7304 

1.1197 

.0268 

.0000 

-32.1662 

-.0452 

-.0526 

-.0061 

-1.2567 

-.7517 

.7154 

32.1473 

-.0169 

-.0788 

-.0747 

-.3803 

.0375 

-.2334 

2.0306 

1.1024 

.4938 

.0094 

-.0536 

-.0037 

-2.9979 

-.5308 

.4155 

.0000 

.0000 

.0076 

.0040 

-.0036 

.0710 

-.5943 

.0013 

.0000 

.0000 

.0226 

.0151 

-.0007 

.4058 

,4069 

-.4940 

.0000 

.0000 

.0000 

.0000 

.0000 

1.0000 

.0005 

-.0154 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.9994 

.0343 

.0000 

.0000 

Dimensional  8th  Order  Design  Model  B  Matrix 

B8  = 


-4.9064 

-.9103 

30.4980 

-.0834 

-3.9661 

30.7245 

.5566 

18.8079 

34.5965 

-.3527 

.4905 

.0099 

-5.3885 

47.5204 

1.3582 

9.9301 

-.7124 

.5789 

-8.4358 

-.0700 

12.8683 

-5.9781 

3.1746 

-11.9831 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

Dimensional  8th  Order  Design  Model  C  Matrix 

C8  = 


0. 

0. 

1. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

1. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

1. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

1. 

0. 
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APPENDIX  B 

FULL  STATE  CLOSED  INNER  LOOP  TRANSFER  FUNCTION 

MATRIX 
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Full  State  Closed  Inner  Loop  Transfer  Function  Matrix  for 
the  8th  Order  Design  Model 


Column  1 

Column  2 

UJ 

3.992 

Ul 

0.077 

wc 

(-4) 

Pc 

(-4X-.00526X-.0001) 

p 

3.29E-5(-.l)(- .0000171) 

P 

4 

Wc 

(-4X-.00015X-.002) 

Pc 

(-4) 

q 

1 .1 74E-4(-.00784)(. 00497) 

q 

5.554E-4(-5.1  76E-6) 

wc 

(-4X-.002X- .00526) 

Pc 

(-4)(-.00526)(.0001 ) 

r 

-7.956E-5(.007;.1  1) 

r 

9.896E-4(-.001 43) 

wc 

(-4X-.002X-.00015) 

Pc 

(-4X-.00526X-.0001) 

Column  3 

Column  4 

uu 

2745(-.0357) 

UJ 

-.0186(1 .122X-.0214) 

^c 

(-4)(-.00015)(-,002) 

f'c 

(-4)(- .0001  5)(- .002) 

p 

5.554E-4(- .00001 542) 

p 

9 ,896E-4(3 .5456X-9 .504E-6) 

qc 

(-4X-.002X-.0001 5) 

rc 

(-4)(-.00015)(-.002) 

q 

3.9778(.0077) 

q 

-4.91  E-4(-.0042)(4.265E-5) 

% 

(-4X-.002) 

rc 

(-.0001 5X-.002X- .00526) 

r 

-4 .909E~4(  .44 ;  .07) 

r 

3. 9978(. 000226) 

(-4X-.002X-.00015) 

r'c 

✓“N 

1 

A 

W 

/“\ 

1 

o 

o 

o 

Ul 

V _ / 

Mote 

:  (  X  )  for  real  eigenvalues  ; 

)  for  complex  eigenvalues 
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